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PREFACE 


Cumberland  Island  National  Seashore  was  established  in  1972  to 
preserve  the  scenic,  scientific  and  historical  values  of  the 
largest  and  most  southerly  island  off  the  coast  of  Georgia.  It  is 
well  known  for  its  marine  turtles,  abundant  shorebirds,  dune 
fields,  maritime  forests,  fishing,  marshes  and  tidal  creeks  and 
flats,  and  historic  structures.  The  St.  Marys  Inlet,  at  the 
border  between  Florida  and  Georgia,  is  a  Federally  maintained 
entrance  channel  to  the  Intracoastal  Waterway,  ports  at 
Fernandina,  Florida,  and  St.  Marys,  Georgia,  and  the  US  Naval 
Submarine  Base  at  Kings  Bay,  Georgia.  Construction  of  coastal 
engineering  works  and  channel  dredging  over  the  past  100  years 
have  had  noticeable  effects  on  the  St.  Marys  Entrance,  Cumberland 
Island,  Georgia  and  Amelia  Island,  Florida.  In  the  early  1960 's, 
Kings  Bay  was  selected  as  the  Navy's  home  port  for  Poseidon-class 
submarines.  In  the  mid-1970' s  Kings  Bay  was  selected  to  homeport 
the  Navy's  new  class  of  Trident  submarines.  In  upgrading  the 
Kings  Bay  base  from  the  smaller  Poseidon  submarines,  it  was 
necessary  to  deepen,  widen,  and  lengthen  the  entrance  channel  to 
Kings  Bay.  The  five  year,  Kings  Bay  Environmental  Research 
Program  was  conceived  in  1986  by  the  U.S.  Departments  of  Interior 
and  Navy.  This  Department  of  the  Navy  funded  Program  focuses  on 
evaluating  the  potential  effects  on  the  natural  resources  of 
Cumberland  Island  and  vicinity  of  the  deepening  of  the  Kings  Bay 
Trident  Submarine  ship  channel  from  42  ft.  (12.7  m)  to  51  ft. 
(15.5  m)  .  The  channel  is  almost  22  miles  (35.2  km)  long  and 
required  the  removal  of  approximately  3  5  million  cubic  yards 
(26.8  million  cu  m)  of  dredged  material.  The  potential 
biophysical  effects  of  dredging  are  being  evaluated  by  the 
National  Park  Service  through  a  series  of  biological  and 
geological  research  projects.  The  Department  of  the  Navy,  through 
the  US  Army  Engineers,  is  monitoring  the  physical  aspects  of  the 
ocean  shoreline  of  Cumberland  and  Amelia  Islands  and  the 
Cumberland  Sound  estuary.  Technical  direction  and  guidance 
during  the  study  were  provided  by  Dr.  Albert  Greene,  Jr., 
National  Park  Service  (NPS) ;  Messrs.  Thomas  J.  Peeling,  Naval 
Facilities  Engineering  Command  (NAVFAC) ;  John  Headland,  NAVFAC; 
Darryll  Molzan,  NAVFAC;  Dr.  Robert  Dean,  University  of  Florida, 
Gainesville  (NPS) ;  Dr.  Stephen  Cofer-Shabica  (NPS) ;  and  the  late 
Dr.  William  Odum,  University  of  Virginia,  Charlottesville  (NPS) . 

The  ultimate  goal  of  this  research  is  to  document  the  potential 
for  short-  and  long-term  changes  on  the  resources  of  Cumberland 
Island  and  Cumberland  Sound  estuary  related  to  channel  dredging. 
The  work  described  in  this  report  is  one  of  a  series  of  National 
Park  Service  studies  directed  towards  this  goal. 
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ABSTRACT 

A  hydrogeologic  investigation  of  southern  Cumberland  Island  was  conducted  to 
determine  the  effect  of  channel  dredging  on  the  groundwater  quality  in  the  Pliocene- 
Miocene  aquifer.  As  a  part  of  the  study  10  wells  were  installed  on  the  southern  portion  of 
the  island.  Three  wells  were  installed  at  two  cluster  sites  and  four  wells  were  installed  at  a 
third  cluster  site. 

The  aquifers  of  interest  included  the  Miocene  sand  aquifer,  the  Pliocene-Miocene 
aquifer,  and  the  surficial  aquifer.  Water-level  measurements  during  high  tide  indicate  that 
the  horizontal  gradient  of  the  surficial  aquifer  was  toward  the  southwest  and  the  Pliocene- 
Miocene  aquifer  was  toward  the  west.  The  direction  of  groundwater  flow  within  the 
Pliocene-Miocene  aquifer  may  change  in  response  to  fluctuating  tides  and  storm  surges. 

An  aquifer  test  performed  in  the  Pliocene-Miocene  aquifer  indicated  that  the 
transmissivity  of  the  aquifer  ranged  from  235  to  650  ft2/d,  the  hydraulic  conductivity 
ranged  from  34  to  94  ft/d,  the  storage  coefficient  ranged  from  1.05  x  10-5  to  5.6  x  10-5, 
and  the  diffusivity  ranged  from  5.9  x  103  to  3.5  x  107  ft2/d.  Test  data  indicated  that  lateral 
and  vertical  heterogeneities  existed  in  the  aquifer.  The  seepage  velocity  of  the  Pliocene- 
Miocene  aquifer  ranged  from  0.049  to  0.085  ft/d 

The  surficial  aquifer  was  in  a  steady-state  condition  with  constant  long-term 
recharge,  equivalent  discharge  and  no  appreciable  withdrawals  due  to  pumping. 
Approximately  2.85  x  105  gallons  per  square  mile  recharged  the  aquifer  each  day  and 
nearly  the  same  amount  was  discharged  to  the  low  lying  wetlands  (1.92  x  105  gal/d),  the 
underlying  Pliocene-Miocene  aquifer  (8.5  x  104  gal/d),  and  along  the  coast  at  the  seepage 
face  (7.58  x  103  gal/d). 

Analysis  of  hydraulic  data  suggest  that  intrusion  of  seawater  from  the  channel  into 
the  Pliocene-Miocene  aquifer  was  insignificant.  The  amount  of  freshwater  discharged  from 
the  aquifer  during  flushing  was  an  order  of  magnitude  greater  than  the  volume  of  seawater 
that  intruded  the  aquifer.  The  net  volume  of  freshwater  discharged  from  the  Pliocene- 
Miocene  aquifer  per  unit  area  was  0.019  ft3/ft2  d  (0.142  gal/ft2  d). 
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INTRODUCTION 

During  the  mid-seventies,  the  U.S.  Navy  initiated  a  channel-dredging  program  to 
facilitate  Trident  submarine  traffic  to  and  from  the  Kings  Bay  Naval  base,  located  on  the 
mainland  near  St.  Marys,  Ga.  (Figure  1).  The  U.S.  Army  Corps  of  Engineers  drilled  over 
300  boreholes  in  the  Kings  Bay  ship  channel  to  determine  the  shallow  sediment  and  rock 
composition  and  the  depth  to  the  top  of  rock.  Based  on  their  log  descriptions  of  1984,  the 
amount  of  overburden  covering  the  Miocene  limestone  ranged  from  none  in  some  areas 
near  the  center  of  the  channel  to  as  much  as  54  feet*  along  the  banks.  The  average 
thickness  of  the  overburden  was  1 8  feet  At  that  time  the  depth  to  the  bottom  of  the  channel 
ranged  from  17  to  48  feet  with  an  average  depth  of  30  feet.  The  Navy  has  since  dredged 
the  channel  to  a  depth  of  51  feet  below  mean  sea  level  (Figure  2).  This  depth  of  dredging 
removed  clay  and  silt  deposits  that  protected  the  freshwater  aquifers  from  seawater 
intrusion  and  in  many  areas  breached  the  Miocene  limestone. 

The  U.S.  Departments  of  the  Interior  and  Defense,  in  cooperation  with  Georgia 
State  University,  developed  a  5-year  cooperative  program  to  evaluate  the  effects  of  channel 
dredging  on  the  hydrogeological  system.  This  thesis  describes  the  results  of  the 
hydrogeologic  investigation  of  Cumberland  Island  and  determines  the  effect  of  channel 
dredging  on  the  Pliocene-Miocene  aquifer. 

Purpose  and  Scope 
The  purpose  of  this  investigation  was  to  determine  the  effect  of  channel  dredging  on 
the  quality  of  water  in  the  Pliocene-Miocene  aquifer.  As  part  of  the  study,  the  geologic  and 
hydrogeologic  framework  of  southern  Cumberland  Island  was  characterized  by: 
completing  a  literature  and  data  search;  installing  new  wells;  describing  core  and 
*1  foot  =  0.3048  meters 
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Figure  1.  Cumberland  Island,  Ga.,  and  surrounding  geographical  features. 
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Figure  2.  Location  and  depth  of  channel  dredging  by  the  US  Navy  and  well  site  locations. 


cutting  samples;  conducting  surface  and  borehole  geophysics;  and  constructing  structure 
and  potentiometric  maps.  The  hydraulic  properties  of  the  Pliocene-Miocene  aquifer  were 
determined  through  aquifer  testing.  Water  level  data  and  aquifer  test  results  were  used  to 
determine  aquifer  interaction  and  the  potential  for  seawater  intrusion  into  the  Pliocene- 
Miocene  aquifer  from  Cumberland  Sound  Channel. 

Description  of  Study  Area 
Cumberland  Island  is  located  in  the  Barrier  Island  sequence  of  the  Sea  Island 
section  of  the  Coastal  Plain  physiographic  province.  It  is  a  compound  barrier  island  on  the 
Georgia  coast  just  north  of  the  Georgia- Florida  State  line  (Figure  1).  The  island  is  19  miles 
long  and  3  miles  wide  at  the  broadest  point  having  an  area  of  about  60  square  miles 
(McLemore  and  others,  1981).  It  is  separated  from  Jekyll  Island,  Ga.,  by  the  St.  Andrews 
Sound  to  the  north  and  from  Ameila  Island,  Fla.,  by  the  Cumberland  Sound  to  the  south. 
The  study  area  extends  from  the  south  end  of  Cumberland  Island  to  the  area  surrounding 
Dungeness  and  encompasses  an  area  of  9  square  miles  (Figure  1).  Cumberland  Island  is 
administered  as  a  National  Seashore  by  the  National  Park  Service. 

Climate 
The  island's  climate  is  characterized  by  long,  humid  summers  and  mild,  dry 
winters.  The  mean  annual  rainfall  is  about  52  inches  of  which  40  inches  are  lost  to 
evapotranspiration,  6  inches  to  storm  runoff  and  the  remaining  6  inches  are  available  for 
recharge  (Brown,  1984).  From  the  months  of  June  through  September  the  island  receives 
60  percent  of  the  annual  rainfall.  The  dry  season  occurs  from  October  to  May  (Brown, 
1984). 


PREVIOUS  WORK 

Martinez  (1980)  described  the  Neogene  stratigraphy  of  Cumberland  Island 
including  lithologic  descriptions,  paleontology,  and  thicknesses  of  each  zone.  Data  were 
collected  from  core  and  cuttings  from  18  boreholes  and  seismic  profiles  taken  from  the 
periphery  of  the  island.  McLemore  and  others  (198 1)  provided  a  brief  description  of  the 
stratigraphy  from  late  Eocene  to  Holocene  age  and  discussed  the  geomorphology  of  the 
island  The  report  also  evaluated  the  hydrogeology  of  the  island  with  emphasis  on  the 
Floridan  aquifer  system  but  did  not  define  the  hydraulic  properties  of  these  aquifers.  The 
water  quality  of  the  Floridan  and  surficial  aquifers  was  included. 

Brown  (1984)  defined  the  hydrogeologic  framework  of  northeastern  Florida  and 
southeastern  Georgia.  The  majority  of  the  data  were  collected  from  Fernandina  Beach, 
Fla.,  and  Sl  Marys,  Ga.,  areas.  Cumberland  Island  was  included  in  the  study  area,  but 
the  hydrogeology  of  the  island  was  not  investigated.  The  report  contained  a  detailed 
description  of  the  hydrogeologic  characteristics  of  the  Floridan  aquifer  system.  Also,  a 
range  of  estimated  transmissivity  values  were  given  for  the  surficial  aquifer  (100  to  1000 
square  feet  per  day  (ft2/d)  based  on  lithology)  and  the  Pliocene-Miocene  aquifer  (750  ft2/d 
based  upon  aquifer  testing)  at  Kingsland,  Ga.,  15  miles  from  the  study  area.  The  recharge 
and  discharge  areas  for  the  upper  and  lower  Brunswick,  Pliocene-Miocene,  and  surficial 
aquifers  were  discussed  for  the  region. 

Miller  (1986)  described  the  hydrogeologic  framework  of  the  Floridan  aquifer 
system,  and  Bush  and  Johnston  (1988)  determined  the  ground- water  hydraulics  and 
regional  flow  in  its  entire  aerial  extent  Krause  and  Randolph  (1989)  defined  the 
hydrogeology  of  the  Floridan  aquifer  system  and  evaluated  its  development  potential  using 
digital  modeling.  Although  these  investigations  included  Cumberland  Island,  they  were 


part  of  a  regional  aquifer  evaluation  and  were  not  specific  to  the  hydrogeology  of  the 
island. 

Clarke  and  others  (1990)  provided  a  detailed  description  of  the  late  Eocene  through 
late  Miocene  age  deposits  along  the  coast  of  Georgia.  In  their  report,  core,  cuttings,  and 
geophysical  logs  were  described  from  a  well  on  the  island.  The  report  also  described  the 
Floridan  aquifer  system,  the  upper  and  lower  Brunswick  aquifers,  and  the  surficial  aquifer 
along  the  coast  of  Georgia.  A  map  of  the  potentiometric  surface  of  the  Floridan  aquifer 
(May  1985)  at  Cumberland  Island  was  included. 

In  summary,  much  hydrogeologic  information  is  available  for  the  Floridan  aquifer 
system  in  the  southeast  Georgia,  northeast  Florida  area.  However,  information  is  limited 
for  the  aquifers  above  the  Floridan  aquifer  system  in  the  Cumberland  Island  area. 
Information  pertaining  to  the  characteristics  of  those  aquifer  systems  overlying  the  Floridan 
on  Cumberland  Island  generally  does  not  exist. 


GEOLOGIC  FRAMEWORK 

Geomorphology 

The  geomorphology  of  Cumberland  Island  is  typical  of  barrier  islands  along  the 
Georgia  coast.  From  seaward  to  landward  margin,  the  geomorphic  provinces  are  the 
beach,  foredunes,  interdune  meadow,  back  dunes,  interior  forest,  and  marshes. 
McLemore  and  others  (1981)  gave  a  thorough  description  of  the  Cumberland  Island 
geomorphology.  The  following  discussion  is  a  brief  synopsis  of  their  description. 

The  beaches  of  Cumberland  Island  are  composed  of  fine  quartz  sand.  The  beaches 
average  125  feet  in  width  and  are  as  much  as  300  feet  wide  on  the  southeast  coast  due  to 
construction  of  a  jetty  on  the  south  end  of  the  island  in  1876.  The  landward  margin  is  not 
as  expansive  as  the  seaward  margin.  The  landward  margin  is  composed  of  Pleistocene 
sands  eroded  from  the  island  and  clay  derived  from  Cumberland  Sound.  Some  areas  along 
the  seaward  margin  consist  of  the  oyster  Crassostrea  virginica.  However,  these  areas  are 
not  extensive. 

The  foredunes  lie  behind  the  eastern  beach  and  extend  the  length  of  the  island.  The 
dunes  vary  in  height  from  1.5  feet  to  6  feet  and  consist  of  cross-bedded  fine  quartz  sand 
and  some  heavy  minerals.  In  general  these  dunes  erode  during  the  winter  and  rebuild  in  the 
summer. 

Interdune  meadows  are  broad  basin  areas  adjacent  to  the  foredunes.  The  meadows 
receive  sediment  from  wind-blown  sand  off  the  beach  and  foredunes  and  periodically  by 
washover  as  a  result  of  storm  surges.    The  meadows  are  composed  of  well-sorted,  fine 
quartz  sand.  Sediment  remains  in  the  interdune  meadow  only  if  stabilized  by  vegetation,  as 
is  the  case  on  the  south  end  of  the  island. 

Back  dunes  lie  behind  the  interdune  meadow  over  most  of  the  island's  length  and 


extend  to  a  maximum  height  of  45  feet.  They  are  composed  of  well-sorted,  cross-bedded, 
fine  quartz  sand  and  heavy  minerals.  The  back  dunes  migrate  to  the  west  because  of  the 
prevailing  easterly  winds.  Migration  may  be  rapid  where  there  is  a  lack  of  vegetation. 

Live  oaks  and  palmettos  of  the  maritime  forest  prevent  rapid  migration  of  dunes 
along  the  Raccoon  Keys  salt  marsh.  The  undergrowth  of  the  interior  maritime  forest 
consists  almost  entirely  of  saw  palmetto.  The  forest  lies  on  flat  rolling  Pleistocene  sand 
and  is  the  prominent  feature  on  the  island. 

The  salt  marshes  dominate  the  west  side  of  Cumberland  Island  which  is  typical  of 
Georgia  barrier  islands.  These  marsh  deposits  consist  of  clay  sediment  in  low-lying  areas 
and  silty  sand  in  elevated  areas  such  as  on  hammocks.  The  freshwater  inland  marshes 
form  in  areas  of  low  topography  such  as  Johnson  Pond,  or  where  drainage  is  impeded  by 
dunes  along  the  east  coast  such  as  Lake  Whitney  and  the  Sweetwater  Complex  (see  Figure 
1).  These  marshes  are  not  influenced  by  tidal  currents  (McLemore  and  others,  1981). 

Regional  Island  Geology 
Cumberland  Island  is  underlain  by  a  thick  sequence  of  consolidated  and 
seimconsolidated  coastal  plain  sediments.  Lithological  and  geophysical  data  exists  for  a 
few  of  the  deep  wells  (approximate  total  depth  of  600  feet)  on  the  island  (see  Figure  3  for 
well  locations).  Stratigraphy  below  this  depth  are  not  addressed  in  this  thesis;  only  the  late 
Eocene  through  Holocene  age  deposits  are  discussed  herein. 

The  late  Eocene  age  unit  consists  of  the  Ocala  Limestone,  the  top  of  which  is 
located  548  feet  below  land  surface  (see  Figure  4).  The  Ocala  is  a  massive  fossiliferous 
limestone  that  contains  bryozoan  remains,  foraminifera  tests,  and  mollusk  shells.  This 
limestone  formed  in  a  warm,  shallow- water,  carbonate-bank  environment  (Clarke  and 
others,  1990). 
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Bate  from  McLemore  and  others  (1981). 


Figure  3.  Location  of  wells,  Cumberland  Island. 
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The  Oligocene  series  is  located  at  a  depth  of  53 1  feet  below  land  surface  and  consists  of  a 
17-foot  section  of  the  Suwannee  Limestone  (see  Figure  4).  The  Suwannee  is  a  phosphatic 
limestone  and  is  distinguished  by  an  abundance  of  miliolid  foraminifera.  The  limestone 
was  deposited  in  a  carbonate  bank  environment  (Clarke  and  others,  1990). 

Deposits  of  Miocene-age  consist  of  the  Charlton  and  Hawthorn  Formations  and  are 
composed  of  sand,  clay,  limestone,  and  dolomite  (Clarke  and  others,  1990).  The  Miocene 
strata  are  about  458  feet  thick. 

Clarke  and  others  (1990)  divided  the  Miocene  strata  into  three  units  based  on 
paleontologic  evidence  and  geophysical  characteristics.  Each  unit  consists  of  a  basal 
carbonate  layer,  overlain  by  clay,  and  an  upper  sand  layer.  In  descending  order,  these 
units  are  Miocene  units  A,  B,  C  (see  Figure  4).  Miocene  unit  C  is  about  91  feet  thick  in  the 
study  area.  The  basal  carbonate  layer  of  unit  C  is  composed  of  sandy,  fossiliferous, 
phosphatic  limestone,  and  dolomite.  It  grades  upward  into  a  silty  clay  and  clayey  silt  that 
contains  molds  of  diatoms.  The  clay  layer  grades  upward  into  a  poorly  sorted,  very  fine  to 
granule-size  quartz  sand  containing  sparse  phosphate  grains  and  dolomite  rhombohedrons. 
Miocene  unit  B  is  about  66  feet  thick  and  Miocene  unit  A  is  about  215  feet  thick  in  the 
study  area.  The  composition  of  Miocene  units  A  and  B  are  similar  to  unit  C.  However,  the 
sand  unit  of  Miocene  unit  A  contains  alternating  beds  of  argillaceous  dolostone  (Clarke  and 
others,  1990).  The  section  of  the  Charlton  Formation  examined  during  this  study  was 
equivalent  to  Miocene  unit  A.  Figures  6  and  7  show  geologic  sections  from  the  top  of  the 
Charlton  Formation  to  land  surface  (see  Figure  5  for  cross  section  location).  (Figure  8 
shows  the  altitude  of  the  top  of  the  Charlton  Formation).  According  to  Henry  and  Kellam 
(1988),  the  middle  Miocene  age  sediments  were  deposited  in  an  open  marine  environment . 
The  core  taken  from  Site  1  (see  Figure  3  for  site  location)  contained  abundant  marine 
fossils  and  evidence  of  bioturbation  indicating  that  a  shallow  open  marine  or  marginal  sea 
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EXPLANATION 


B'      Alignment  of  Geologic  Section 


Data  point  identification  number 

6 

•     Borehole  drilled  by  the  Georgia 

Geologic  Survey  (McLemore  and 
others,  1981) 
SITE  3 

KB  MP  10  Number  on  bottom  is  the  well  drilled 
0  for  the  Kings  Bay  Monitoring  Project. 

Number  on  top  is  the  site  location. 
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KBMP 


Base  Irom  McLemore  and  othe.s  (1981) 


Figure  5.  Location  of  boreholes  and  geologic  sections,  Cumberland  Island. 
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EXPLANATION 


-60  —  STRUCTURE  CONTOUR-- 

Shows  altitude  of  top  of  Charlton 
Formation,  in  feet.  Contour  interval, 
10  feet.  Datum  is  mean  sea  level. 
DATA  POINT- 

Borehole  drilled  by  the  Georgia  Geologic 
Survey  (McLemore  and  others,  1981) 

Well  drilled  for  the  Kings  Bay 
Monitoring  Project 

Number  is  altitude  of  top  of 
Charlton  Formation,  in  feet. 
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Bas*  from  McL»more  and  oth«rs  (1981). 


Figure    8 .  Altitude  of  the  top  of  the  Charlton  Formation. 
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environment  existed  during  deposition. 

In  the  study  area  the  middle  Pliocene  series  is  the  equivalent  of  the  Duplin 
Formation  of  southeastern  Georgia  and  consists  of  well- sorted,  calcareous  sand  which  is 
locally  gravelly  (Henry  and  Kellam,  1988).  The  thickness  of  the  formation  ranges  from  2 
to  42  feet  thick  (McLemore  and  others,  1981).  The  early  Pliocene  sediments  were 
deposited  in  a  restricted  basin  and  estuarine  environment.  These  deposits  were  subaerially 
eroded  and  channel  cut  during  middle  Pliocene  (Henry  and  Kellam,  1988).  (Figures  9  and 
10  show  the  altitude  and  thickness  of  the  Duplin  Formation.) 

The  Pleistocene  age  deposits  consist  of  22  to  62  feet  of  sand  and  clayey  sand 
overlying  5  to  40  feet  of  clay  and  sandy  clay.  The  Pleistocene  sand  is  the  Satilla  Formation 
equivalent  and  consists  of  very  fine  quartz  sand.  This  sediment  represents  barrier-island 
facies  deposits  and  ranges  in  color  from  brown  in  the  soil  horizon  to  light  brown  in  deeper 
zones.  The  fossil  burrow  Ophiomorpha  nodosa  of  the  shrimp  Callianassa  major  is  present 
which  indicates  that  deposition  occurred  in  a  near- shore  environment  (McLemore  and 
others,  1981).  Below  the  sand  deposits  are  deposits  which  represent  back-barrier 
sedimentation  and  consist  of  silt,  clay,  and  sandy  clay,  that  range  from  gray  to  tan  in  color. 
(See  Figures  1 1  and  12  for  the  altitude  and  thickness  of  the  clay  sediments  of  Pleistocene 
age.)  Commonly,  the  clay  is  interbeded  with  sand  ranging  from  less  than  1  feet  to  about  10 
feet  in  thickness.  The  lithology  here  indicates  an  estuarine  back-barrier  depositional 
environment  interrupted  by  periods  of  high  energy  (Martinez,  1980). 

The  Holocene  deposits  occur  as  fine  quartz  sand  in  the  beach  and  dune  areas  along 
the  east  coast.  Along  the  west  coast,  the  deposits  consist  of  clays  at  the  base  of  the 
freshwater  ponds  and  marshes.  On  the  island's  interior  the  deposits  occur  as  floodplain 
and  stream  deposits  (McLemore  and  others,  1981). 
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EXPLANATION 


STRUCTURE  CONTOUR-- 
Shows  altitude  of  top  of  Duplin 
Formation,  in  feet.    Contour 
interval,  10  feet.   Datum  is 
mean  sea  level. 

DATA  POINT- 

Borehole  drilled  by  the  Georgia  Geologic 
Survey  (McLemore  and  others,  1981) 

Well  drilled  for  the  Kings  Bay 
Monitoring  Project 

Number  is  altitude  of  top 
of  Duplin  Formation,  in  feet. 
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Base  from  McL»mor»  and  oth»r»  (1981). 


Figure  9.  Altitude  of  the  top  of  the  Duplin  Formation. 
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EXPLANATION 


— 10 STRUCTURE  CONTOUR-- 

Shows  thickness  of  Duplin 
Formation,  in  feet    Contour 
interval,  10  feet    Dotted 
where  approximately  located. 
Datum  is  mean  sea  level. 

DATA  POIMT-- 

Borehole  drilled  by  the  Georgia  Geologic  Survey 
(McLemore  and  others,  1981) 

Well  drilled  for  the  Kings  Bay  Monitoring 
Project 

Number  is  thickness 
of  Duplin  Formation,  in  feet. 
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Ba»»  from  Mclsmore  and  othert  (19B1). 


Figure  10.  Thickness  of  Duplin  Formation. 
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EXPLANATION 


•STRUCTURE  CONTOUR-- 

Shows  altitude  of  top  of  clay  sediments 
of  Pleistocene  age,  in  feet.  Contour 
interval,  10  feet  Datum  is  mean  sea  level. 
DATA  POINT- 

Borehole  drilled  by  the  Georgia  Geologic  Survey 
(McLemore  and  others,  1981) 

Well  drilled  for  the  Kings  Bay  Monitoring 
Project 

Number  is  altitude  of  top  of  clay 
sediments  of  Pleistocene  age,  in  feet. 
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Baa*  from  McLemore  and  others  (1981). 


Figure  1 1.  Altitude  of  the  top  of  the  clay  sediments  of  Pleistocene  age. 
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EXPLANATION 


—20  —  STRUCTURE  CONTOUR-- 

Shows  thickness  of  the  clay  sediments 
of  Pleistocene  age,  in  feet  Contour 
interval,  10  feet.  Datum  is  mean  sea  level 
16  DATAPOINT- 

•      Borehole  drilled  by  the  Georgia  Geologic  Survey 

(McLemore  and  others,  1981) 
32 
o      Well  drilled  for  the  Kings  Bay  Monitoring  Project 

Number  is  thickness  of  clay  sediments 
of  Pleistocene  age,  in  feet. 
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Base  from  McL»mor«  and  oth»r«  (1981). 


Figure  12.  Thickness  of  clay  sediments  of  Pleistocene  age. 
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Local  Geology 
Geology  from  the  upper  Charlton  Formation  through  the  Holocene  deposits  were 
investigated  in  the  study  area.  The  following  sections  include  a  site  specific  discussion  of 
the  Charlton  Formation,  the  Duplin  Formation,  and  Holocene  to  Pleistocene  deposits 
undifferentiated.  Core  samples  were  examined  at  Site  1;  cutting  samples  were  examined  at 
Sites  1  and  2.  A  detailed  description  of  the  core  and  cuttings  are  provided  in  Appendix  B. 

Miocene  Series 

The  Charlton  Formation  consists  of  alternating  layers  of  sand,  clay,  limestone,  and 
dolomite  locally  characterized  by  high  concentrations  of  phosphate.  At  Site  1,  the  lower 
part  of  the  Charlton  Formation  (151  to  141  feet  below  sea  level)  consists  of  argillaceous 
dolomitic  limestone.  The  light  gray  section  is  soft,  bioturbated  with  shell  molds  (20 
percent).  The  dolomitic  limestone  is  overlain  by  10  feet  of  light  gray,  poorly  indurated,  to 
loose  quartz  sand.  At  Site  1,  the  SPR  log  was  useful  for  determining  the  contacts  between 
the  quartz  sand  and  the  dolomitic  limestone  (see  Figure  13).  The  sand,  subrounded  to 
subangular,  varies  from  fine  to  medium  (25  percent)  to  very  coarse  (60  percent).  Light 
gray  dolomitic  limestone  (31  ft  thick)  overlies  the  quartz  sand.  This  section  consists  of  fine 
to  medium  sand  (35  percent),  clay  (35  percent),  and  very  coarse,  phosphatized  calcite 
grains  as  large  as  0.25  inches  in  diameter.  The  section  is  highly  bioturbated  and  is  locally 
high  in  phosphate  content.  Three  feet  of  very  light  gray,  calcareous  clay  separates  the 
dolomitic  limestone  from  8  feet  of  dolomite.  The  dolomite  is  argillaceous  (30  percent  clay) 
with  very  fine  to  fine  sand  (35  percent)  and  very  coarse  quartz  sand  and  phosphatized 
calcite  (5  percent).  The  overlying  zone  (89  to  74  feet  below  sea  level)  consists  of 
alternating  layers  of  dolomite,  sand,  dolomitic  limestone  and  limestone  (Figure  13).  This 
section  consists  of  fine  to  medium  sand  (45  percent),  very  coarse  quartz  grains  (7  percent) 
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Figure  13.  Lithology  and  geophysical  logs  at  well  KBMP  1,  site  1. 
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and  weathered  shell  fragments.  It  is  highly  pitted,  weathered,  fractured,  and  has  good 
secondary  porosity.  Dolomite  overlying  the  89  to  74  foot  section  increases  in  clay  content 
(25  percent)  and  is  not  as  fractured  and  weathered  as  the  underlying  unit.  The  dolomite 
consists  of  fine  sand  (45  percent)  and  becomes  very  hard  and  dense  at  the  top  of  the 
Charlton  Formation,  61  feet  below  sea  level. 

At  Site  2,  the  upper  109  feet  of  the  Charlton  Formation  were  examined  as  cutting 
samples.  Variations  in  the  lithology  do  not  exist  as  are  present  at  Site  3.  The  section 
consists  of  olive  gray  dolomitic  limestone  (Figure  14).  The  dolomitic  limestone  is 
argillaceous  (30  percent  clay)  with  fine  to  medium  sand  (40  percent)  subangular,  broken 
shell  fragments  (5  percent)  and  locally  contains  very  coarse  sand  to  gravel  (15  percent). 
The  gamma  log  indicates  that  this  zone  is  highly  phosphatic  (see  Figure  15). 

At  Site  3,  an  upper  53  foot  section  of  the  Charlton  Formation  was  examined. 
Figure  16  is  a  combination  of  data  from  Site  3  and  lithologic  descriptions  of  boreholes  by 
Martinez  (1980)  that  were  drilled  to  the  west  of  Site  3  .  The  basal  45  feet  of  this  section 
consist  of  alternating  layers  of  argillaceous  dolomite,  clayey  sand,  and  sandy  clay.  This 
zone  contains  phosphatic  calcite  and  abundant  phosphate  at  depths  of  100  and  125  feet,  as 
indicated  by  the  sharp  increase  in  natural  gamma  emission  (see  Figure  17).  The  light  gray 
dolomite  layers  are  very  hard  and  consist  of  clay  (30  percent)  and  fine  to  coarse  sand  (20 
percent),  subangular  to  subrounded.  The  sandy  clay  layers  are  about  40  percent  clay  and 
30  percent  fine  to  coarse  sand,  subangular  to  subrounded,  with  some  shell  fragments.  The 
upper  8  feet  of  the  Charlton  Formation  consists  of  weathered,  fractured  limestone  with  fine 
to  coarse  sand  (15  percent)  and  shell  fragments  (15  percent).  This  layer  is  continuous  on 
the  island  north  of  Dungeness  based  on  lithologic  descriptions  by  Martinez  (1980). 
Core  taken  from  the  Charlton  Formation  contained  abundant  marine  fossils  such  as 
Mercenaria prodroma  cf.,  Chesapecten  cf.,  and  the  brachiopod  Lingula.  Callianassa 
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'Originally  well  KBMP  6  was  drilled  to  a  total  depth  of  1 70  ft-  below  land  surface.  Because  of  partial  borehole  collapse,  the 
borehole  was  fj Ued,  abandoned,  and  a  new  well  was  drilled  to  a  total  depth  of  7 1  ft  below  land  surface. 


Figure  14.  Geologic  section  D  to  D'. 
cross  section  location. 


Refer  to  Figure  5  for 
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Figure  15.  Lithology  and  geophysical  logs  at  well  KBMP  6,  site  2. 
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Figure  17.  Lithology  and  geophysical  logs  at  well  KBMP  10,  site  3. 
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burrows  cf.  were  found  in  the  core  as  well.  These  fauna  indicate  that  a  shallow  open 
marine  or  marginal  sea  environment  existed  during  deposition. 

Pliocene  Series 

The  Pliocene  deposits  on  Cumberland  Island  are  equivalent  to  the  Duplin 
Formation.  In  general  the  formation  consists  of  marine  sand  that  fines  upward  (McLemore 
and  others,  1981).  The  Pliocene  deposits  at  sites  1  and  2  are  much  different  than  the 
deposits  at  Site  3  and  those  described  by  McLemore  and  others  (1981)  northward  on  the 
island.  At  these  sites,  the  sand  contains  a  significant  amount  of  clay  and  the  sediment  is 
phosphatic  as  indicated  by  the  active  gamma  emission  (see  Figures  13  and  15). 

At  Site  1  the  deposits  are  26  feet  thick  and  consist  of  brown  gray  to  olive  brown 
clayey  sand,  very  fine  to  fine  (Figure  14).  These  deposits  are  slightly  calcareous  and 
consist  of  shell  fragments  ranging  from  a  trace  to  20  percent.  At  Site  2,  the  deposits  are  12 
feet  thick  and  consist  of  clayey  sand,  very  fine  to  fine  with  shell  fragments  (10  percent). 

At  Site  3,  an  unconformable  contact  exists  between  the  weathered  limestone  of  the 
Charlton  Formation  and  sand  of  the  Duplin  Formation.  The  basal  8  feet  of  the  Duplin 
Formation  consists  of  sandy  clay  having  very  coarse  quartz  sand  and  abundant  broken 
pectin  fragments  (20  percent),  which  is  overlain  by  10  feet  of  alternating  layers  of  sand  and 
clayey  sand  with  20  percent  shell  fragments  (see  Figure  16).  The  8-foot  sandy  clay  layer 
was  not  described  by  Martinez  (1980),  thus  it  likely  is  a  localized  lens  and  probably 
represents  stream  deposits.  The  stratigraphy  and  broken  pectin  shells  found  in  the  Duplin 
Formation  indicate  a  near  shore  environment 
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Holocene  and  Pleistocene  Deposits 

There  are  no  distinguishing  characteristics  between  the  Holocene  and  Pleistocene 
sand  on  the  south  end  of  the  island.  Therefore,  these  deposits  are  referred  to  as  Holocene 
to  Pleistocene  deposits  undifferentiated.  At  Site  1,  the  basal  clayey  zone  (17  feet  thick) 
consists  of  alternating  layers  of  sandy  clay  and  clay.  The  medium  gray,  quartz  sand  is  very 
fine,  and  the  dark  gray  clay  is  clean.  The  overlying  yellowish  gray  sand  (23  feet  thick)  is 
very  fine,  subangular,  and  well  sorted.  At  Site  2,  the  lower  part  of  the  Holocene  to 
Pleistocene  deposits  consist  of  light  gray  sandy  clay  with  very  fine  to  fine  sand,  (35 
percent)  and  heavy  minerals  (3  percent).  The  overlying  olive  gray,  quartz  sand  (45  feet 
thick)  is  very  fine  to  fine,  and  well  sorted. 

At  Site  3,  the  lower  part  of  the  Holocene-Pleistocene  deposits  (32  feet  thick)  consist 
of  alternating  layers  of  light  gray  sandy  clay,  clean  clay,  and  yellowish  gray  clayey  sand, 
very  fine.  These  deposits  represent  back  barrier-island  deposits.  The  clayey  zone  is 
overlain  by  30  feet  of  well  sorted,  very  fine,  quartz  sand  which  represent  barrier-island 
deposits.  Clay  is  completely  absent  in  the  upper  section  of  the  Holocene-Pleistocene 
deposits  at  each  site  (see  Figures  14  and  16). 
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METHODS  OF  INVESTIGATION 

Various  methods  were  employed  to  collect  data  needed  to  characterize  the 
geology/hydrogeology  of  the  study  area.  Surface  geophysical  data  provided  preliminary 
geologic  and  water  quality  data  that  were  used  to  select  the  most  appropriate  locations  for 
the  well  sites.  A  total  of  10  wells  were  drilled  at  the  well  sites.  Three  wells  were  emplaced 
in  a  cluster  at  well  Site  1;  three  wells  at  Site  2;  and  four  wells  at  Site  3  (see  Figure  3  for 
location  of  well  sites).  Core,  split-spoon  samples,  and  cuttings  were  collected  and 
analyzed  during  well  drilling.  Borehole  geophysical  logs  enhanced  core  and  cuttings  data 
by  providing  information  such  as  contacts  between  formations,  porosity  and  lithology  of 
the  aquifers,  and  marker  beds  for  the  tops  of  aquifers.  Surface  geophysical  surveys 
provided  specific  conductivity  data  of  the  shallow  sediments  and  pore  water.  The  hydraulic 
characteristics  of  the  aquifer  were  determined  by  performing  an  aquifer  test. 

Surface  Geophysics 
An  electromagnetic  (EM)  terrain  conductivity  meter  (the  Geonics  EM  34-3)  was 
used  to  determine  the  specific  conductivity  of  the  shallow  sediments  and  formation  water  in 
the  study  area.  The  instrument  measures  conductivity  in  millimhos  per  meter  at  theoretical 
depths  of  25  feet  to  197  feet  below  land  surface  depending  on  the  coil  configuration  and 
spacing  (McNeil,  1980).  Three  transects  were  made  across  the  southern  end  of  the  island. 
At  each  location  of  EM  conductivity  measurement,  a  reading  for  the  33,  66,  and  132  feet 
intercoil  spacing  for  each  dipole  configuration  (horizontal  and  vertical)  was  taken  so  that 
data  acquisition  ranged  from  theoretical  depths  of  25  to  197  feet  below  land  surface.  The 
conductivity  meter  was  useful  for  detecting  cavities  in  limestone,  mapping  lithology,  and 
locating  the  zone  of  diffusion  between  the  freshwater  and  seawater. 
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The  conductivity  meter's  transmitter  coil  receives  an  alternating  audio  frequency 
current  which  produces  a  magnetic  field  in  the  coil.  The  primary  field  generated  by  the  coil 
propagates  through  the  geologic  medium  creating  electrical  currents  in  the  strata.  The 
amount  of  current  flow  depends  on  the  conductivity  of  the  layers  (Zohdy  and  others, 
1974).  A  secondary  magnetic  field  is  produced  that  has  the  same  frequency  but  a  different 
phase  than  that  produced  by  the  transmitter  coil.  The  secondary  field  is  detected  by  the 
receiver  coil,  which  cancels  signals  from  the  primary  field.  The  strength  of  the  secondary 
field  is  directly  proportional  to  terrain  conductivity  which  is  a  function  of  intercoil  spacing, 
operating  frequency,  and  ground  conductivity  (Duran,  1984). 

Most  of  the  current  flow  in  a  saturated  zone  is  along  the  surface  of  the  grains  and 
through  the  pore  water.  The  bulk  conductivities  of  saturated  sedimentary  rocks  are  based 
largely  on  the  effective  porosity  and  formation- water  conductivity.  As  these  increase  so 
does  the  bulk  conductivity.  For  strata  saturated  with  saline  water,  fluid  properties  dominate 
the  EM  conductivity  response.  In  sediment  saturated  with  freshwater,  the  porosity  controls 
the  bulk  conductivity  (Stewart,  1982). 

Drilling  Methods 

Three  well  sites  were  selected  in  the  study  area  to  provide  hydrogeological, 
stratigraphic  and  geophysical  data  in  areas  representing  diverse  hydrogeologic  conditions 
and  to  monitor  the  water  levels  and  water  chemistry  (see  Figure  3  for  well  site  locations). 
At  each  site,  a  cluster  of  three  to  four  wells,  spaced  about  10  feet  apart,  were  drilled  to 
monitor  the  water-bearing  zones  of  interest. 

Drilling  of  well  KBMP  1  (Kings  Bay  Monitoring  Project)  was  initiated  without  the 
use  of  a  drilling  additive  or  hole  stabilizers.  After  drilling  12  feet,  the  hole  began  to 
collapse  and  a  5  1/2  inch  diameter  fishtail  bit  and  a  drilling  additive  was  used  to  finish  the 
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hole  to  the  top  of  rock.  A  biodegradable  drilling  mud  was  used  which  naturally 
decomposes  after  72  hours.  Split  spoon  samples  were  taken  at  5  foot  intervals  to  the  top  of 
the  Charlton  Formation  (66  feet  below  land  surface).  The  remainder  of  the  well  was  cored 
using  a  steel  rock  bit  core  barrel  to  a  depth  of  146  feet  with  an  average  core  recovery  of 
94.8  percent.  At  the  remaining  wells,  cuttings  were  taken  and  analyzed  at  5  to  10  foot 
intervals  or  where  a  significant  change  in  lithology  occurred.  The  core,  split-  spoon 
samples,  and  cuttings  were  analyzed  for  lithological  and  textural  composition  under  a 
binocular  microscope.  Ubiquitous  phosphate  in  samples  was  detected  by  directly  applying 
ammonium  molybdate.  If  phosphate  was  present,  the  sample  turned  bright  yellow.  The 
relative  permeability  was  estimated  based  on  sediment  grain  size  and  sorting,  and  for  rock 
zones,  the  amount  of  weathering  and  secondary  porosity.  Dolomite  was  distinguished 
from  dolomitic  limestone  by  noting  the  abundance  of  dolomite  rhombohedrons.  Dolomite 
was  very  dense,  hard,  and  reacted  very  slowly  to  hydrochloric  acid.  Dolomitic  limestone 
was  medium  hard,  a  lighter  color,  and  not  as  dense  as  the  dolomite.  Limestone  reacted 
violendy  to  hydrochloric  acid. 

The  remaining  wells  were  drilled  using  the  hydraulic  rotary  method  and  usually  a  7 
5/8  inch  diameter  bit.  Where  temporary  casing  was  used  to  prevent  caving,  a  9  3/4  inch  bit 
was  used  to  the  top  of  the  Charlton  Formation  and  finished  with  a  7  5/8  inch  bit. 

Well  Design  and  Construction 

All  wells  were  cased  and  screened  using  4  inch  polyvinyl  chloride  (pvc),  and  were 
capped  at  the  base  (see  Figures  18  through  20).  The  annulus  of  each  well  was  packed  with 
sand  to  a  distance  of  about  2  to  5  feet  above  the  screen.  A  2-foot  bentonite  plug  was  placed 
above  the  sand  pack  to  seal  the  annular  space,  except  in  shallow,  sandy  zones  where  a  plug 
was  not  necessary,  or  where  the  annular  space  caved  in  before  a  plug  could  be  installed, 
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Figure  18.  Well  construction  at  well  Site  1.  The  altitude  is  4.59  feet  above  mean  low 
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Originally  well  KBMP  6  was  drilled  to  a  total  depth  of  170  ft.  below  land  surface.  Because  of 
partial  borehole  collapse,  the  borehole  was  filled,  abandoned,  and  a  new  well  was  drilled  to  a 
total  depth  of  71  ft.  below  land  surface. 


Figure  19.  Well  construction  at  well  Site  2.  The  altitude  is  6.86  feet  above  mean  low  water. 
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Figure  20.  Well  construction  at  well  Site  3.  The  altitude  is  15.9  feet  above  mean  low  water. 
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such  as  in  well  KB  MP  4.  Above  the  bentonite  plug,  the  remaining  annular  space  was  filled 
with  cement  grout  to  seal  the  annular  space  and  prevent  leakage  from  overlying  or 
underlying  zones.  Table  1  summarizes  the  construction  of  each  well. 

A  grain-size  analysis  was  conducted  for  the  sandy  aquifers  to  determine  suitable 
screen  slot  sizes.  The  screen  slot  size  chosen  for  each  well  was  selected  to  retain  40 
percent  of  the  sediment  and  permit  60  percent  to  pass  through  the  screen  (see  Figure  21). 
This  will  promote  grading  around  the  screened  area  and  increase  the  porosity  and 
permeability  of  the  near  bore  formation  during  well  development.  As  a  result,  the 
drawdown  around  the  well  will  be  reduced  during  pumping  and  the  well  efficiency  will  be 
increased  (Driscoll,  1984). 

Appropriate  slot  sizes  for  each  well  were  selected  by  permitting  60  percent  of  the 
aquifer  medium  to  pass  through  the  screen  during  development.  However,  an  exception  to 
this  procedure  was  made  for  well  KBMP  9.  According  to  the  grain-size  analysis,  a  50-slot 
screen  (width  of  the  slot  is  in  thousandths  of  an  inch)  would  have  been  appropriate  for  this 
zone.  Because  a  screen  size  this  large  was  not  available  at  the  time  of  drilling,  the  largest 
available  screen  slot  size,  12  slot,  was  installed. 

The  slot  size  selected  for  monitor  wells  tapping  limestone  aquifers  was  based  on  the 
percentage  and  size  of  sand  present  in  the  aquifer.  The  saturated  zones  tapped  by  well 
KBMP  3  and  KBMP  7  contained  a  significant  amount  of  fine  sand,  therefore  12-slot 
screens  were  installed.  Wells  KBMP  4  and  KBMP  10  were  installed  in  an  argillaceous 
limestone  and  10  and  12  slot  screens  were  installed. 

Well  Development 

All  wells  installed  on  Cumberland  Island  as  a  part  of  this  investigation  were 
developed  to  reduce  compaction  of  the  formations  caused  during  drilling,  remove  drilling 
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fluid  from  the  formation,  and  increase  the  permeability  and  effective  porosity  around  the 
well,  which  in  turn  increases  the  yield  and  specific  capacity.  The  wells  were  developed  by 
overpumping  with  a  4-inch  submersible,  0.5  horsepower  pump  for  varying  times  ranging 
from  45  minutes  to  2.5  hours  at  approximately  20  gal/min.  Well  KBMP  10  could  only  be 
pumped  at  5  minute  intervals  because  the  zone  tapped  by  this  well  produced  litde  water. 

During  overpumping,  the  pump  was  periodically  shut  off  and  restarted  to  create  a 
surging  effect.  One  drawback  of  overpumping  is  that  water  moves  only  toward  the  well, 
so  backwashing  is  not  possible.  The  unidirectional  movement  of  water  may  cause 
sediment  to  bridge  outside  the  screen.  If  the  formation  is  agitated  during  normal  pumping, 
sediment  may  collapse  and  enter  the  well.  However,  if  unstable  bridging  occurred,  the 
sand  pack  placed  around  the  screen  will  prevent  most  of  the  sediment  from  entering  the 
well  (Driscoll,  1986). 

Borehole  Geophysics 
At  each  drilling  site,  borehole  geophysical  logs  were  derived  using  an  EG&G 
Mount  Sorpris  logger  from  the  deepest  well  of  the  cluster.  Natural  gamma,  single-point 
resistance,  and  spontaneous  potential  logs  were  run  at  each  site.  Figures  13,  15,  and  17  in 
the  geology  section  contain  the  geophysical  data  collected  at  each  site. 

Natural  Gamma 

The  natural  gamma  log  measures  the  natural  radioactive  emissions  from  the 
geologic  material.  The  gamma  probe  consists  of  a  scintillation  receiver  and  counter  circuit. 
Gamma  kick  intensity  in  a  geologic  medium  is  proportional  to  the  number  of  pulses 
detected  by  the  probe  per  unit  time.  The  intensity  is  expressed  as  the  average  number  of 
counts  per  second  (Driscoll,  1986).  In  sedimentary  rocks,  clay  and  shale  usually  emit 
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higher  counts  of  gamma  rays,  due  to  their  high  potassium  content,  than  do  limestone  and 
sand,  which  are  characterized  by  low  levels  of  potassium,  thorium,  and  uranium. 
However,  micaceous  sand  emits  high  levels  of  radiation  as  do  dolomite  and  phosphorite 
(Hilchie,  1982).  Much  of  the  radiation  emitted  from  carbonates  and  phosphates  originate 
from  thorium  and  uranium  particles.  Some  phosphate  pebbles  may  contain  up  to  150  parts 
per  million  (ppm)  of  uranium,  and  therefore  are  gamma  active  zones  (Kwader,  1982).  At 
Site  1,  the  high  intensity  kicks  recorded  by  the  gamma  probe  did  not  match  zones  in  the 
core  that  were  high  in  phosphate.  The  difference  between  the  depth  read  by  the  probe  and 
the  actual  depth  of  the  phosphatic  zone  appeared  to  be  5  feet  (a  kick  at  75  feet  below  land 
surface  on  the  log  compared  to  the  lithology  at  70  feet). 

Electric  Logs 

The  single-point  resistance  (SPR)  and  spontaneous  potential  (SP)  logs  were  run 
simultaneously.  Instead  of  using  the  conventional  system  in  which  the  probe  contains  a 
down-hole  lead  electrode  and  a  ground  electrode  placed  at  the  surface,  the  differential  SPR 
and  SP  log  system  was  used.  In  this  system  both  electrodes  are  included  in  the  probe  and 
are  separated  by  an  insulator.  The  current  flows  around  the  insulator  from  the  lead 
electrode  to  the  probe  shell.  This  system  provides  higher  resolution  than  the  conventional 
system  (Keys,  1988). 

Spontaneous  Potential 

The  SP  log  measures  changes  in  natural  electrochemical  potential  between  the 
borehole  fluid,  formation  fluid,  and  geologic  material  (Kwader,  1982).  To  interpret  the  SP 
log,  a  clean  clay  or  shale  should  be  present  in  the  strata  so  that  a  clay  baseline  can  be 
established.  Without  a  defined  clay  baseline,  interpretation  of  the  log  is  almost  impossible. 
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The  SP  deflections  are  to  the  left  of  the  baseline  where  the  formation  water  is  more  saline 
than  the  drilling  mud  and  are  to  the  right  if  the  formation  water  is  fresher  than  borehole 
water  (Hilchie,  1982).  If  the  formation  and  borehole  water  are  both  fresh,  insignificant 
deflection  occurs  and  the  log  will  be  practically  a  straight  line. 

The  SP  logs  run  in  the  wells  on  Cumberland  Island  tended  to  drift  with  deflection 
that  yielded  little  information.  The  geologic  sequence  logged  was  mostly  calcareous  sand, 
sandy  clay,  and  argillaceous  limestone.  Clean  clays  were  not  present  in  these  sections. 
Although  a  comparatively  pure  clay  layer  did  exist  at  two  locations  logged,  temporary  metal 
casing  was  in  place  during  logging  and  extended  below  these  zones.  Therefore,  these 
cased  zones  could  not  be  logged  with  the  electric  log,  and  a  clay  baseline  could  not  be 
established. 

SP  drifting  also  occurred  because  the  formation  water  alternated  from  saline  to 
fresh.  The  direction  of  deflection  is  dependent  on  the  relation  of  formation  and  borehole 
water  chemistry,  and  an  informative  curve  could  not  be  generated  due  to  the  changing 
formation  water  chemistry.  As  a  result,  the  SP  logs  taken  on  Cumberland  Island  yielded 
little  hydrogeological  information. 

Single-Point  Resistance 

The  SPR  log  measures  the  electrical  resistance  of  the  geologic  formation  and 
borehole  water  between  two  electrodes.  The  ability  of  a  rock  to  conduct  a  current  depends 
on  the  effective  porosity  and  conductivity  of  the  water  in  the  pore  spaces.  Resistivity 
increases  as  porosity,  water  content,  and  hydraulic  conductivity  decrease  (Driscoll,  1986). 
The  SPR  log  was  useful  in  providing  information  such  as  location  of  contacts  between 
formations,  relative  water  quality,  and  the  distinction  of  sandy  zones  and  sandy  limestones 
from  argillaceous  limestones  and  dolostones. 
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Grain  Size  Analysis 

A  grain  size  analysis  of  samples  representing  each  sandy  zone  (the  Pliocene  sand, 
Miocene  sand  aquifer,  and  the  surficial  aquifer  at  well  Sites  1,  2,  and  3)  was  performed  to 
determine  the  appropriate  screen  sizes  for  the  new  wells,  and  estimate  the  hydraulic 
conductivity  of  these  zones  (Figure  21).  The  sandy  zones  were  well  sorted  quartz  sand  and 
contained  little  to  no  clay  deposits.  Table  Al  in  Appendix  A  shows  the  interval  sampled  in 
each  zone  and  the  raw  data  from  the  analysis. 

Each  sample  (35-45  grams)  was  oven  dried  for  2  hours  and  weighed  to  the  nearest 
0.01  gram.  The  'nest'  of  sieves  ranged  from  -2  to  4  0  (0  =  -log2  d,  where  d  is  the  particle 
diameter  in  millimeters).  The  samples  were  vibrated  for  10  minutes  on  80  percent  of  the 
power  setting.  The  fraction  retained  on  each  mesh  was  weighed  and  recorded.  Finally  the 
measurement  error  in  recovery  percent  of  the  initial  sample  weight  was  determined  (Table 
Al).  The  grain  size  distribution  from  wells  KBMP  1,  2,  5,  8,  and  9  are  given  in  Figure 
21. 

Aquifer  Test 
A  48-hour  aquifer  test  was  performed  at  Site  3  in  the  Pliocene-Miocene  aquifer  near 
Dungeness.  This  aquifer  was  chosen  for  testing  because  it  may  be  susceptible  to  seawater 
intrusion  due  to  channel  dredging.  Prior  to  the  aquifer  test,  a  step-drawdown  test  was 
performed  to  determine  the  well  efficiency  and  the  optimum  pumping  rate  for  the  aquifer 
test. 
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Step  Drawdown  Test 

The  step-drawdown  test  was  conducted  by  selecting  five  pumping  rates  or  steps 
each  lasting  1.5  hours,  and  measuring  drawdown  in  the  pumped  well.  The  data  provided 
from  the  test  can  be  used  to  determine  the  specific  capacity  of  the  well  and  the  optimum 
pumping  rate  (Jacob,  1947).  During  the  step-drawdown  test,  a  centrifugal  pump  having  25 
feet  of  intake  hose  and  200  feet  of  discharge  hose  was  used.  The  discharge  hose  was 
connected  to  an  orifice  bucket  that  measured  the  discharge  rate.  The  orifice  bucket  was 
frequently  monitored  to  ensure  a  constant  pumping  rate.  Each  successively  higher 
pumping  rate  was  implemented  by  increasing  the  throttle  20  percent. 


Multi-well  Aquifer  Test 

The  aquifer  test  was  conducted  December  9-11,  1989,  at  Site  3  on  Cumberland  Island 
(Figure  3).  The  pumping  well  (KB MP  7)  penetrated  the  limestone  part  of  the  Pliocene- 
Miocene  aquifer  so  that  the  transmissivity  of  this  zone  could  be  determined.  During  the 
aquifer  test,  six  wells  were  monitored,  three  of  which  penetrated  the  limestone  part  of  the 
aquifer.  The  pumping  well  was  pumped  for  48  hours  using  a  4-inch  submersible  pump 
powered  by  a  generator.  The  pump  discharge  hose  was  5/8  of  an  inch  in  diameter  and  200 
feet  long.  This  distance  for  discharge  prevented  direct  local  recharge  of  the  surficial  aquifer 
during  the  test  The  discharge  was  monitored  for  the  pumping  duration  by  an  orifice  bucket. 

The  pumping  and  observation  wells  were  monitored  for  the  duration  of  the  test  by 
continuous  water-level  recorders,  pressure  transducers,  and/or  the  wetted- tape  method. 
The  water-level  recorders  were  installed  on  4  of  the  monitoring  wells.  During  the  test,  the 
recorders  were  set  to  collect  data  for  12  hour  time  periods.  The  cylindrical  drum  of  the 
recorder  is  connected  to  a  pulley  with  a  steel  tape  threaded  through  it.  The  steel  tape  has  a 
float  attached  to  one  end  and  a  counter  weight  on  the  other.  As  the  water  level  fluctuates  in 
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the  well,  the  tape  causes  the  drum  to  rotate.  A  pen,  controlled  by  a  clock,  moves  across  the 
rotating  drum  with  time,  thus  producing  a  hydrograph. 

The  transducer  is  a  variable  reluctance  device  operated  in  conjunction  with  a  data 
logger- 120.  Deflection  of  the  transducer  diaphragm  changed  the  reluctance  of  the 
transducer  and  produced  a  voltage  proportional  to  the  pressure  difference,  or  the  water  level 
difference,  across  the  transducer.  The  voltage  was  digitally  displayed  on  the  data  logger  in 
feet  Hollow  tubing  from  the  port  of  the  transducers  to  the  data  logger  allowed  the 
diaphragm  to  be  exposed  to  atmospheric  pressure  changes  and  hence  negated  these  effects. 
A  transducer  was  used  in  the  pumping  well  and  three  of  the  monitoring  wells.  Periodic 
water-level  measurements  were  made  using  the  wetted-tape  method  for  verification. 

During  the  aquifer  test,  two  monitoring  wells  were  located  where  the  pressure 
transducers  could  not  reach.  These  wells  were  measured  solely  by  the  wetted-tape  method. 
Water  levels  measured  for  the  potentiometric  surface  maps  on  June  28,  1990  also  were 
measured  by  this  method. 

Barometric  pressure  data  were  collected  from  the  National  Oceanic  and  Atmospheric 
Administration  at  Jacksonville,  Fla.,  and  St.  Simons  Island,  Ga.,  to  evaluate  the  effects  of 
changes  in  barometric  pressure  on  water  levels  during  the  aquifer  test.  The  data  were  used 
to  correct  water-level  fluctuations  due  to  barometric  pressure  changes. 
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ANALYSIS  OF  SURFACE  GEOPHYSICS  AND  WELL  SITE  LOCATION 

EM  conductivity  surveys  were  conducted  to  determine  the  most  appropriate  well 
site  locations.  Areas  with  low  specific  conductivity  (around  25  millimhos/meter)  were 
optimal  locations  for  new  wells.  Low  specific  conductivity  values  indicate  clean,  highly 
resistive  strata,  saturated  with  freshwater,  such  as  sand  or  limestone.  These  zones  were 
ideal  for  monitoring  water  quality  changes.  At  each  location  of  EM  conductivity 
measurement,  a  reading  for  the  33,  66,  and  132  foot  intercoil  spacing  for  each  dipole 
configuration  (horizontal  and  vertical)  was  taken  so  that  data  acquisition  ranged  from 
theoretical  depths  of  25  to  197  feet  below  land  surface  (see  Figures  22  through  27). 

The  conductivity  values  changed  depending  on  the  intercoil  spacing  (or  depth)  for  a 
given  location.  However,  these  differences  were  slight  and  it  is  suspected  that  the 
theoretical  depth  of  reading  and  the  actual  depth  may  not  be  the  same,  especially  in  the  salt 
marsh  areas  where  much  of  the  shallow  sediment  is  clay.  This  lithology,  combined  with 
the  high  water  salinity,  may  have  prevented  the  current  from  propagating  to  lower  zones. 
In  the  marsh  areas,  it  is  likely  that  only  the  shallow  zones  were  explored  even  at  the 
maximum  coil  spacing. 

Although  changes  in  conductivity  were  observed  at  a  given  location  for  varying 
intercoil  spacings  or  depth  of  exploration,  the  differences  in  conductivity  were  much  more 
dramatic  for  changes  in  geomorphology,  geology,  and  water  chemistry.  Low  conductivity 
values  were  observed  in  sandy,  freshwater  zones  such  as  along  the  dunes;  high  values 
were  observed  in  saline-saturated  clay  deposits  such  as  in  marsh  areas.  Therefore,  similar 
trends  of  conductivity  occur  as  a  function  of  location  rather  than  varying  depths  at  a  single 
location. 

Readings  taken  around  Dungeness  were  low  for  each  intercoil  spacing  and 
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■  Lin©  of  Equal  Conductivity- 
Interval  varies.  In  millimhos  per  meter 
Data  Point-    Number  Is  conductivity  In  millimhos  per  meter 
Theoretical  depth  of  reading  Is  25  ft  below  land  surface 


Base  from  u  S  Geological  Survey 
Femandina  Beach  Fla  -Ga    1  24  000   1 958  and 
Cumberland  Island  South    124  000    1958 


Figure  22.  EM  conductivity  surface  where  the  coil  spacing  is  33  feet 
and  the  coil  configuration  is  horizontal. 
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—  30  —  Lin©  of  Equal  Conductivity- 
Interval  varies,  in  milllmhos  per  meter 
Data  Point-    Number  Is  conductivity  In  milllmhos  per  meter 

Theoretical  depth  of  reading  Is  50  ft  below  land  surface 


Base  from  U  S  Geological  Survey 
Femandina  Beach,  Fla  -Ga .  1  24,000.  1958  and 
Cumberland  Island  Soulh.  1  24,000.  1958 


Figure 23.  EM  conductivity  surface  where  the  coil  spacing  is  33  feet 
and  the  coil  configuration  is  vertical. 
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Base  Irom  U  S  Geological  Survey 
Fernandlna  Beach.  Fla  -Ga    1  24  000  1958  and 

Cumberland  Island  Souin   1  24,000  1958 


Figure  24.  EM  conductivity  surface  where  the  coil  spacing  is  66  feet 
and  the  coil  configuration  is  horizontal. 


49 


30"4S 


30°43 


+ 

-30 


Line  of  Equal  Conductivity- 
Interval  varies,  In  mllllmhos  per  meter 
Data  Point-    Number  is  conductivity  In  millimhosper  meter 
Theoretical  depth  of  reading  Is  98  ft  below  land  surface 


Base  from  U  S  Geological  Survey 

Femandina  Beach.  Fla.-Ga  .  1  24  000,  1958  and 

Cumberland  Island  South.  1  24.000.  1958 


Figure  25.  EM  conductivity  surface  where  the  coil  spacing  is  66  feet 
and  the  coil  configuration  is  vertical. 
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-30-~  Line  of  Equal  Conductivity- 
Interval  varies,  In  milllmhos  per  meter 
•       Data  Point-    Number  Is  conductivity  in  milllmhos  per  meter 
Theoretical  depth  of  reading  is  98  ft  below  land  surface 


Base  from  u  S  Geological  Survey 

Fernandina  Beach  Ra  -Ga  .  1  24,000.  1958  and 

Cumberland  Island  Souin.  1  24000,  1958 


3000 

=1 


Figure  26.  EM  conductivity  surface  where  the  coil  spacing  is  132  feet 
and  the  coil  configuration  is  horizontal. 
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30°43   - 


Base  from  U  S  Geological  Survey 

Femandlna  Beach.  Fla  -Ga  .  i  24.000.  1958  and 

Cumberland  Island  South.  1  24,000.  1958 


Figure  27.  EM  conductivity  surface  where  the  coil  spacing  is  132  feet 
and  the  coil  configuration  is  vertical. 
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configuration.  The  132  foot  horizontal  and  66  foot  vertical  readings  (89  feet  below  land 
surface)  were  slightly  higher  which  probably  indicates  the  top  of  rock.  Along  the  dunes  the 
conductivity  readings  were  also  low  becoming  slighdy  higher  in  the  interdune  meadows, 
along  the  beach,  and  along  the  west  side  of  the  dunes  adjacent  to  the  marsh.  High 
conductivity  values  were  observed  in  the  low-lying  marsh  areas  (150  to  1000 
millimhos/meter)  which  were  attributable  to  the  presence  of  clay  sediment  saturated  with 
saline  water.  The  only  low  readings  measured  in  the  marsh  were  on  a  hammock  composed 
of  silty  sand.  A  shallow  freshwater  lens  probably  exists  on  the  hammock. 

The  locations  selected  for  the  well  sites  in  the  study  area  were  in  fresh  water  zones. 
Areas  with  low  electric  conductivity  were  selected,  such  as  well  Site  3  at  Dungeness  (see 
Figure  3  for  well  site  locations).  On  the  southern  part  of  the  island  the  conductivity  was 
inconsistent  and  ranged  from  21  to  890  millimhos/meter.  Here  the  geomorphology  varies 
from  low-lying  dunes  to  marsh  areas.  The  lower  values  were  observed  in  the 
topographically  higher  areas  where  freshwater  accumulates  in  silty-sand  deposits,  such  as 
the  areas  selected  for  well  Sites  1  and  2.  Conductivity  increased  at  these  sites  below  depths 
of  50  feet,  but  not  significantly  (see  Figures  22  through  27). 
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HYDROGEOLOGIC  FRAMEWORK 

The  sediments  underlying  Cumberland  Island  comprise  both  local  and  areally 
extensive  aquifers.  For  the  purpose  and  scope  of  this  report,  the  Miocene  sand  aquifer, 
Pliocene-Miocene  aquifer,  and  surficial  aquifer  were  investigated.  The  Floridan  aquifer 
was  beyond  the  scope  of  this  study. 

Miocene  Sand  Aquifer 

In  the  study  area,  aquifers  of  Miocene  age  are  present  from  1 10  to  400  feet  below 
land  surface  and  range  in  thickness  from  10  to  40  feet  (see  Figure  4).  Sand  zones  in  the 
Hawthorn  Formation  have  been  named  the  upper  and  lower  Brunswick  aquifers  (Clark  and 
others,  1990).  Only  the  aquifers  above  these  in  the  Charlton  Formation  were  investigated 
during  this  study,  and  only  the  uppermost  Miocene  sand  aquifer  is  discussed  in  detail. 

The  Miocene  sand  aquifer  is  confined  above  by  alternating  layers  of  argillaceous 
dolomite  and  clay  and  below  by  argillaceous  dolomitic  limestone.  It  is  10  feet  thick,  the  top 
of  which  is  located  136  feet  below  land  surface  at  Site  1  (see  Figure  28).  The  aquifer  is 
recharged  by  precipitation  in  the  Charlton  Formation  outcrop  areas  outside  the  study  area. 
The  areal  extent  of  the  Miocene  sand  aquifer  is  unknown  because  it  was  not  penetrated  by 
other  wells. 

The  horizontal  direction  of  groundwater  flow  in  the  Miocene  sand  aquifer  is  not 
known  because  only  one  well  penetrates  this  zone  in  the  study  area.  The  potential  vertical 
direction  of  groundwater  flow  is  from  the  Miocene  sand  aquifer  to  the  overlying  Pliocene- 
Miocene  aquifer  based  on  water  level  measurements.  At  Site  1 ,  the  head  of  the  Miocene 
sand  aquifer  is  0.13  feet  higher  than  that  of  the  overlying  Pliocene-Miocene  aquifer  based 
on  water-level  measurements  on  June  28,  1990. 
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Clarke  and  others  (1990)  determined  that  the  hydraulic  conductivity  of  the  lower 
Brunswick  aquifer  ranged  from  20  to  57  ft/d  based  on  aquifer  testing.  The  hydraulic 
conductivity  of  the  Miocene  sand  aquifer  as  derived  by  the  method  of  grain  size  distribution 
described  by  Masch  and  Denny  (1966)  is  53  ft/d.  See  Appendix  C  for  the  Masch  and 
Denny  method  of  determining  hydraulic  conductivity  and  calculations. 

The  water  quality  of  the  Miocene  sand  aquifer  is  unsuitable  for  most  uses.  The 
dissolved  solids  content  in  the  aquifer  exceed  State  drinking  standards.  The  chloride 
content  is  6,200  mg/L  having  32  percent  seawater  (Wilson,  1990). 

Pliocene-Miocene  Aquifer 

The  Pliocene-Miocene  aquifer  consists  of  Miocene-age  limestone  and  Pliocene-age 
sand  At  Site  1,  the  Duplin  Formation  is  clayey  sand  resulting  in  low  permeability.  The 
aquifer  consists  of  vuggy,  fractured  dolomite  that  is  enhanced  by  well-developed  secondary 
porosity.  The  Pliocene-Miocene  aquifer  is  present  at  a  depth  of  79  feet  below  land  surface, 
and  it  is  15  feet  thick  (see  Figure  28).  The  aquifer  is  confined  above  by  dense  dolomite  and 
clay  and  below  by  argillaceous  dolomitic  limestone. 

At  Site  2,  the  top  of  the  Pliocene-Miocene  aquifer  is  located  61  feet  below  land 
surface  and  is  33  feet  thick  (see  Figure  28).  The  aquifer  has  low  permeability  due  to  the 
high  clay  content  and  poorly  developed  secondary  porosity.  At  this  site  the  aquifer  consists 
of  dolomitic  limestone  and  is  confined  above  by  sandy  clay  and  clayey  sand  and  below  by 
argillaceous  dolomitic  limestone. 

At  Site  3  and  northward  on  the  island,  the  Pliocene-Miocene  aquifer  consists  of  two 
zones  having  different  lithologies  and  different  hydraulic  properties.  The  upper  zone  of  the 
aquifer  consists  of  sand  of  the  Duplin  Formation  ranging  from  2  to  42  feet  thick.  The 
lower  zone  consists  of  about  10  feet  of  fractured  limestone  of  the  Charlton  Formation. 
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These  two  zones  are  separated  by  8  feet  of  sandy  clay  and  clayey  sand  The  aquifer  is 
confined  above  by  about  30  feet  of  clay  sediments  of  Pleistocene  age  and  below  by 
alternating  layers  of  argillaceous  dolomite  and  clay  (see  Figure  28). 

The  Duplin  Formation  consists  of  sandy  clay,  clayey  sand  lenses,  based  on  core 
and  cutting  samples.  The  formation  increased  in  clay  content  at  Sites  1  and  2  causing  lower 
hydraulic  conductivity  values  at  these  sites.  Variation  in  the  lithology  of  the  Pliocene  age 
deposits  indicates  that  this  zone  is  anisotropic  and  heterogeneous.  The  weathered  and 
fractured  Charlton  limestone  also  was  anisotropic  and  heterogeneous  based  on  aquifer 
testing  results.  The  heterogeneities  were  due  to  thickness  or  facies  changes  in  the 
limestone. 

The  horizontal  groundwater  flow  gradient  is  low  in  the  Pliocene-Miocene  aquifer  as 
shown  by  the  aquifer's  potentiometric  surface  map  of  June  28,  1990  (see  Figure  29).  The 
direction  of  groundwater  flow  was  solved  for  by  applying  the  three  point  method  where  the 
exact  location  of  an  equipotential  line  is  solved  for  by  determining  the  gradient  between  the 
highest  and  lowest  head  values.  Table  2  summarizes  the  ground- water  levels  for  June  28, 
1990  made  during  high  tide  midway  between  spring  and  neap  tides.  Table  3  summarizes 
the  groundwater  densities  at  each  well  (see  Appendix  D  for  density  corrections).  Potential 
horizontal  and  vertical  groundwater  flow  directions  are  based  on  June  28,  1990  data.  For 
wells  tapping  the  Pliocene-Miocene  aquifer  near  Dungeness,  the  Nightingale  well,  and  the 
Sea  Camp  well,  the  water  levels  vary  areally  0.4  feet  per  mile.  The  horizontal  gradient  in 
this  area  is  so  small  that  the  gradient  direction  could  easily  change  if  water-level  and/or  land 
surface  altitude  measurement  errors  are  made.  Because  the  horizontal  gradient  in  this  area 
is  nearly  zero,  horizontal  groundwater  flow  probably  is  influenced  more  by  tidal 
fluctuations  and  storm  surges  than  by  regional  recharge. 

It  seems  logical  that  the  regional  horizontal  groundwater  flow  direction  in  the 
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EXRANATION 


-  2.0 POTENTIOMETRIC  CONTOUR- 

Shows  altitude  at  which  water  would 
have  stood  in  tightly  cased  wells.   Contour 
interval  1  ft.  Datum  is  sea  level. 


DATA  POINT-- 
Number  on  top  is  well  identification 
Number  on  bottom  is  altitude  of  water  level, 
in  feet  above  sea  level. 


30°  43' 


Base  from  U.S.  Geologic*!  Survey 
Fem.mi.rui  Be«:h.  FLA-GA.  1:24.000.  1970 


Figure  29.  Potentiometric  surface  map  of  the  Pliocene-Miocene  aquifer  for  June  28,  1990.  Arrow 
indicates  potential  direction  of  groundwater  flow. 
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Table  2.  Groundwater  level  summary  for  June  28, 1990. 


WeU 
identi- 
fication 

Depth  to 

water 

level,  in 

ft 

Land 

surface 

elevation, 

in  ft 

Measuring 

point,  in 

ft 

Water 

level 

elevation, 

in  ft. 

Corrected* 

water  level 

elevation,  in  ft 

Hydrogeologic 
zone 

KBMPl 
KBMP2 
KBMP3 

KBMP4 

KBMP5 
KBMP6 

KBMP7 

KBMP8 

KBMP9 

KBMPIO 

KBMPll 

KBMP12 

Nightin- 
gale well 

Sea  Camp 
well 

4.4 

4.6 

0.46 

0.66 

1.8 

Miocene  sand 

5.8 

4.6 

1.78 

0.5 

0.5 

Surficial 

5.1 

4.6 

0.65 

0.15 

1.67 

Pliocene- 
Miocene  age 

4.24 

6.86 

1.25 

3.87 

5.52 

Pliocene- 
Miocene  age 

3.37 

6.86 

0.92 

4.41 

4.44 

Surficial 

4.43 

6.86 

1.4 

3.83 

5.32 

Pliocene- 
Miocene  age 

15.01 

15.93 

1.63 

2.55 

2.62 

Pliocene- 
Miocene  age 

11.92 

15.93 

0.92 

4.93 

4.93 

Surficial 

15.06 

15.93 

1.62 

2.48 

2.48 

Pliocene  sand 

14.54 

15.93 

1.19 

2.58 

2.58 

Miocene 

confining 

layer 

16.32 

16.73 

2.14 

2.55 

2.62 

Pliocene- 
Miocene  age 

14.55 

16.42 

0.72 

2.59 

2.59 

Pliocene- 
Miocene  age 

15.15 

17.86 

0.25 

2.96 

2.96 

Pliocene- 
Miocene  age 

16.73 

18.98 

0.47 

2.72 

2.72 

Pliocene- 
Miocene  age 

*  Refer  to  Table  3  for  percent  seawater  in  each  well. 
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Pliocene-Miocene  aquifer  should  be  to  the  southeast,  because  the  aquifer  is  recharged  in 
outcrop  areas  in  the  northwest.  However,  this  direction  of  horizontal  groundwater  flow 
was  not  observed  in  the  field.  The  direction  observed  may  represent  a  local  flow  system  in 
the  aquifer.  The  water  level  at  well  Site  2  is  almost  3  feet  higher  than  at  Sites  1  and  3.  The 
high  water  level  at  Site  2  may  be  due  to  the  high  clay  content  in  the  aquifer  compared  to  the 
lithology  of  the  aquifer  at  Sites  1  and  3.  A  higher  water  level  occurs  in  a  zone  with  low 
hydraulic  conductivity  than  in  a  zone  under  the  same  conditions  with  a  high  hydraulic 
conductivity.  The  high  water  level  at  Site  2  contributes  to  the  westward  direction  of 
horizontal  groundwater  flow  (see  Figure  29).  The  westward  direction  of  groundwater  flow 
indicates  that  recharge  is  occurring  from  offshore.  Adequate  data  is  not  available  to 
determine  if  recharge  is  occurring  from  the  ocean  to  the  Pliocene-Miocene  aquifer  or  if  the 
high  water  level  at  Site  2  is  an  anomaly  and  represents  local  flow  in  the  aquifer. 

The  vertical  groundwater  flow  gradient  for  the  Pliocene-Miocene  aquifer  varies  with 
time  and  space.  At  Site  3,  the  surficial  aquifer  water  level  is  2.45  feet  higher  than  that  of 
the  Pliocene-Miocene  aquifer.  Therefore,  the  vertical  gradient  flow  potential  is  from  the 
surficial  aquifer  downward  to  the  Pliocene-Miocene  aquifer.  At  Sites  1  and  2  the  Pliocene- 
Miocene  aquifer  water  level  is  slightly  higher  than  that  of  the  surficial  aquifer  ranging  from 
0.88  to  1.17  feet.  Here  the  vertical  direction  of  groundwater  flow  is  the  reverse  of  Site  3, 
and  the  vertical  groundwater  flow  potential  is  from  the  Pliocene-Miocene  aquifer  upward  to 
the  surficial  aquifer.  The  water  level  measurements  were  made  June  28,  1990  during  high 
tide.  It  is  possible  that  during  low  tide  the  potential  groundwater  flow  direction  reverses 
and  groundwater  flows  from  the  surficial  aquifer  to  the  Pliocene- Miocene  aquifer. 

The  Pliocene-Miocene  aquifer  receives  recharge  in  the  outcrop  areas  outside  the 
study  area  in  western  Camden  and  eastern  Charlton  counties.  The  aquifer  also  receives 
recharge  by  leakage  from  along  the  St  Marys  River,  from  the  Cumberland  Sound,  and  from 
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the  overlying  surficial  aquifer  and  the  underlying  Miocene  sand  aquifer  depending  on  the 
hydraulic  gradient  between  the  aquifers  (see  Figure  28). 

Regional  recharge  to  the  Pliocene-Miocene  aquifer  is  from  the  Georgia  Coastal 
Plain.  The  area  of  regional  discharge  is  likely  off  the  east  coast  of  Cumberland  Island. 
However,  the  observed  horizontal  flow  direction  is  toward  the  west  and  possibly 
represents  local  groundwater  flow  as  previously  discussed.  The  aquifer  is  intermittendy 
pumped  by  wells  on  the  island,  but  total  withdrawal  is  negligible. 

The  transmissivity  of  the  Pliocene-Miocene  aquifer  ranges  from  235  to  650  square 
feet  per  day  (ft2/d),  and  the  storage  coefficient  from  1.05  x  10~5  to  5.6  x  10"5  based  on 
aquifer  testing  conducted  during  the  investigation.  Only  the  limestone  part  of  the  aquifer 
was  stressed  during  the  aquifer  test;  the  well  screen  did  not  extend  into  the  sand  of  the 
Duplin  Formation.  The  hydraulic  conductivity  for  the  Miocene  limestone  ranges  from  34 
ft/d  to  94  ft/d  based  on  aquifer  testing.  The  hydraulic  conductivity  of  the  Pliocene  sand  is 
100  ft/d  based  on  the  method  of  deriving  this  value  from  grain  size  analysis  as  described  by 
Masch  and  Denny  (1966). 

The  quality  of  water  from  the  Pliocene-Miocene  aquifer  is  suitable  for  drinking. 
However,  near  Dungeness,  salinity  has  increased  in  the  past  2  years  apparendy  due  to 
saltwater  intrusion.  In  December  1988,  KBMP  1 1  began  producing  salty  water  (David 
Hines,  oral  comm.,  1989)  and  now  contains  5.2  percent  seawater.  Well  KBMP  12, 
located  200  feet  away,  penetrates  the  same  aquifer  zone  and  yields  freshwater.  Well 
KBMP  7,  located  247  feet  from  KBMP  11,  yields  water  that  contains  4.2  percent  seawater 
(Wilson,  1990).  For  most  of  the  extent  of  the  island  the  water  meets  all  state  drinking 
water  standards  and  has  low  dissolved  solids  content  (Brown,  1984).  However,  near  the 
coast,  the  water  chemistry  approaches  that  of  seawater  as  was  observed  at  Sites  1  and  2, 
which  produced  water  having  68  and  89  percent  seawater  content,  respectively  (Wilson, 
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1990). 


The  clay  deposits  that  overlie  the  Pliocene-Miocene  aquifer  consist  of  sediments  of 
Pleistocene  age  and  are  approximately  30  feet  thick.  The  zone  varies  from  a  clean  clay  on 
the  south  end  of  the  island  to  a  sandy  clay  on  the  north  end.  This  confining  layer  extends 
the  length  and  width  of  the  island  and  is  present  below  the  mainland,  but  is  not  continuous 
across  the  Kings  Bay  channel  as  a  result  of  dredging,  channel  scour,  or  both. 

Surficial  Aquifer 

The  surficial  aquifer  is  bound  by  the  Cumberland  River  on  the  west,  the  Atlantic 
Ocean  on  the  east,  St  Andrew  Sound  on  the  north,  and  Cumberland  Sound  on  the  south.  It 
is  confined  below  by  5  to  40  feet  of  clay  and  sandy  clay.  The  aquifer  extends  the  length 
and  width  of  the  island  and  ranges  from  15  to  35  feet  in  saturated  thickness. 

The  surficial  aquifer  primarily  consists  of  well-sorted  fine  sand  that  was  deposited 
by  long  shore  wind  and  wave  processes  (Hoyt,  1967).  These  processes  deposit  sediment 
in  a  directional  configuration  that  parallels  the  shore.  This  results  in  an  anisotropic  aquifer 
medium  with  a  higher  hydraulic  conductivity  in  the  lateral  direction  parallel  to  shore.  The 
aquifer  is  homogeneous  based  on  the  lithology. 

A  freshwater  lens  occurs  in  the  topographically  higher  areas  and  generally  beneath 
the  interior  forest  and  back  dunes  based  on  surface  geophysical  surveys.  However,  the 
freshwater  dimensions  are  not  as  extensive  as  the  aquifer  dimensions,  as  indicated  by 
surface  geophysical  data.  Little  freshwater  exists  in  the  marshes  except  at  higher  elevations 
such  as  on  hammocks. 

The  water  table  ranges  from  4  to  20  feet  above  sea  level  on  the  island's  interior  and 
drops  to  near  sea  level  along  the  coast  (Brown,  1984).  Based  on  June  28,  1990  data,  the 
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direction  of  horizontal  groundwater  flow  is  toward  the  southwest  in  the  study  area  (see 
Figure  30).  This  direction  was  solved  for  by  applying  the  3  point  method.  Horizontal 
groundwater  flow  in  the  surficial  aquifer  tends  to  mimic  the  topography,  and  it  is  therefore 
possible  that  the  dunes  along  the  east  coast  of  the  island  act  as  a  groundwater  divide.  If  this 
is  assumed,  the  direction  of  flow  on  the  west  side  of  the  dunes  is  toward  the  west  and  on 
the  east  side  of  the  dunes,  toward  the  east  All  water  levels  in  the  surficial  aquifer  were 
measured  on  the  west  side  of  the  dunes  and  the  direction  of  groundwater  flow  is  toward  the 
southwest.  If  water  level  data  were  available  east  of  the  dunes,  these  data  may  indicate  an 
eastward  direction  of  groundwater  flow. 

The  potential  vertical  groundwater  flow  gradient  at  Site  3  was  from  the  surficial 
aquifer  to  the  Pliocene-Miocene  aquifer  based  on  June  28,  1990  data.  The  surficial  aquifer 
water  level  was  2.45  feet  higher  than  the  head  of  the  Pliocene-Miocene  aquifer.  At  Sites  1 
and  2  the  potential  vertical  gradient  is  reverse;  the  water  level  of  the  Pliocene-Miocene 
aquifer  is  1.17  feet  higher  than  the  surficial  aquifer  at  Site  1  and  0.88  feet  higher  at  Site  2. 

The  vertical  groundwater  flow  direction  of  the  surficial  aquifer  may  change  with  the 
wet  and  dry  seasons.  During  wet  seasons,  the  aquifer  has  a  higher  water  level  and  may  be 
a  source  of  recharge  for  the  Pliocene-Miocene  aquifer.  During  dry  seasons,  the  water  level 
is  lower  and  may  receive  leakage  from  the  Pliocene-Miocene  aquifer.  Also,  on  the  south 
end  of  the  island  where  the  water  level  of  the  Pliocene-Miocene  aquifer  is  approximately 
0.88  feet  higher  than  the  surficial  aquifer,  the  vertical  direction  of  groundwater  flow  may 
reverse  daily  due  to  tides.  Although  the  Pliocene-Miocene  aquifer  fluctuates  from  0.5  to 
0.8  feet  during  tide  cycles,  the  surficial  aquifer  is  not  affected  by  the  tides  because  of  its 
high  storage  coefficient  value. 

The  surficial  aquifer  is  recharged  primarily  by  rainfall  and  occasional  storm 
washover.  The  average  annual  rainfall  for  the  area  is  52  inches.  Approximately  40  inches 
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Figure  30.  Potentiometric  surface  map  of  the  surficial  aquifer  for  June  28,  1990.  Arrow  indicates 
potential  direction  of  groundwater  flow. 
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are  lost  to  evapotranspiration,  6  inches  to  surface  run  off,  and  6  inches  are  available  to 
recharge  the  aquifer  (Brown,  1984).  The  topsoil  does  not  impede  the  downward 
groundwater  flow  of  rainfall  to  the  surficial  aquifer  because  it  is  only  2  to  3  inches  thick, 
relatively  immature,  and  highly  permeable.  Surface  runoff  is  low  as  indicated  by  the  lack 
of  interior  freshwater  streams. 

The  water  level  in  the  surficial  aquifer  changed  little  from  June  1989  to  June  1990. 
The  greatest  water-level  fluctuation  observed  was  a  rise  of  1.73  feet  from  June  1989  to 
December  1989  due  to  recharge  from  precipitation  during  the  fall.  Assuming  a  porosity  of 
30  percent  for  silty  sand  (Freeze  and  Cherry,  1979),  similar  to  that  of  the  surficial  aquifer 
in  the  area,  then  the  6  inches  of  recharge  should  produce  a  1 .67  foot  rise  in  water-level. 
This  value  was  somewhat  lower  than  the  observed  water-level  fluctuation.  The  amount  of 
recharge  was  an  estimated  value  and  may  be  higher  for  the  year.  Also,  the  actual  volume  of 
water  available  as  recharge  may  be  higher  as  a  result  of  an  antecedent  soil-moisture  surplus. 

The  surficial  aquifer  hydraulic  conductivity  values  derived  by  the  method  described 
by  Masch  and  Denny  (1966)  varied  little  between  each  well  site.  The  values  obtained  from 
samples  from  wells  KB  MP  2,  KB  MP  5,  and  KBMP  8  were  22.7,  19.4  and  22  ft/d, 
respectively.  The  relevant  calculations  are  described  in  Appendix  C. 

The  quality  of  water  from  the  surficial  aquifer  is  suitable  for  drinking  and  domestic 
use,  but  could  easily  become  contaminated  because  the  water  table  is  near  land  surface  and 
has  little  overburden  to  filter  the  recharge  (Brown,  1984).  The  total  dissolved  solids  of 
water  from  wells  tapping  the  surficial  aquifer  in  the  area  averages  398  mg/L.  However,  on 
the  south  end  of  the  island  the  total  dissolved  solids  have  been  reported  as  high  as  1,318 
mg/L  (Wilson,  1990). 
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Aquifer  Interaction 

The  Miocene  sand,  Pliocene-Miocene,  and  surficial  aquifers  in  the  study  area 
interact  by  leakage  through  confining  layers.  If  a  difference  in  head  exists  between  two 
aquifers,  groundwater  has  the  potential  to  move  vertically  from  one  aquifer  to  another  in 
the  direction  of  the  aquifer  having  the  lower  potential  head.  Appendix  E  provides  the 
leakage  equation  and  calculations. 

Based  on  water  levels  of  June  28,  1990,  at  Sites  1  and  2  the  groundwater  head  of 
the  Pliocene-Miocene  aquifer  was  0.88  to  1.17  feet  higher  than  that  of  the  surficial  aquifer. 
The  potential  groundwater  flow  direction  is  from  the  Pliocene-Miocene  to  the  surficial 
aquifer.  Because  the  water-levels  in  the  Pliocene-Miocene  aquifer  fluctuate  from  0.5  to  0.8 
feet  during  a  tide  cycle,  the  direction  of  vertical  leakage  may  reverse  during  low  tide  (see 
Figure  31).  Therefore,  during  high  tide,  the  Pliocene-Miocene  aquifer  is  likely  a  potential 
source  for  the  surficial  aquifer,  and  during  low  tide  it  is  a  potential  sink  for  the  surficial 
aquifer.  The  calculated  leakage  rate  values  between  the  surficial  and  Pliocene-Miocene 
aquifer  at  Sites  1  and  2  were  2.09  x  10'5  ft/d  and  3.17  x  10^  ft/d,  respectively. 

At  Site  3  the  head  in  the  surficial  aquifer  was  2.3 1  feet  higher  than  that  of  the 
Pliocene-Miocene  aquifer.  The  aquifers  are  separated  by  about  30  feet  of  clay.  Here  the 
calculated  leakage  rate  from  the  surficial  aquifer  to  the  Pliocene-Miocene  aquifer  was  4.08  x 
10-4  ft/d. 

At  Site  1  the  head  in  the  Miocene  sand  aquifer  was  only  0.13  feet  higher  than  the 
head  in  the  Pliocene-Miocene  aquifer.  These  aquifers  are  separated  by  42  feet  of 
argillaceous  dolomite.  The  leakage  rate  from  the  Miocene  sand  aquifer  to  the  Pliocene- 
Miocene  aquifer  was  3.1  x  10-6  ft/d  (see  Appendix  E  for  calculations). 
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Analytical  Model 

An  analytical  model  was  developed  to  mathematically  describe  the  hydrogeologic 
groundwater  flow  regime  of  the  surficial  and  Pliocene-Miocene  aquifers.  A  dynamic 
groundwater  flow  system  exists  in  the  study  area  due  to  fluctuating  tides  which  cause 
change  in  the  vertical  and  horizontal  groundwater  flow  directions.  The  equations  and 
schematic  diagram  (Figure  32)  represent  the  groundwater  flow  system  at  a  time-average 
mean  where  the  water  levels  are  at  a  constant  altitude. 

For  this  study  the  analytical  model  describes  the  groundwater  flow  system  in  the 
area  of  Dungeness  because  here  the  geology  and  hydrogeology  are  typical  for  a  large  part 
of  the  island.  The  recharge,  discharge,  and  volumetric  leakage  rates  for  a  1  square  mile 
area  were  calculated  as  part  of  the  analytical  model.  The  calculations  for  the  model  are 
provided  in  Appendix  E. 

The  surficial  aquifer  is  in  a  steady-state  condition  with  constant  long-term  recharge, 
equivalent  discharge,  and  no  appreciable  withdrawals  due  to  pumping.  Approximately  6 
inches  are  available  to  recharge  the  aquifer  each  year  (104  million  gallons  of  precipitation 
per  square  mile  per  year).  About  2.85  x  105  gallons  per  square  mile  recharge  the  aquifer 
each  day  and  nearly  the  same  amount  is  discharged  to  a  combination  of  the  low  lying 
wetlands  (1.92  x  105  gal/d),  the  underlying  Pliocene-Miocene  aquifer  (8.5  x  104  gal/d), 
and  along  the  coast  at  the  seepage  face  (7.58  x  103  gal/d).  See  Appendix  E  for 
calculations. 

The  equation  governing  surficial  aquifer  discharge  along  the  seepage  face  is  derived 
from  the  Dupuit  theory  of  free  surface  flow  (Raudkivi  and  Callander,  1976)  and  is 

expressed  as: 

q=[K  (h12_^_h22i]  (1) 

2L 

where 
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q  =  specific  discharge  [ft/d] 

hi  =  water  level  [ft], 

h2  =  water  level  at  the  seepage  face  [ft], 

K  =  hydraulic  conductivity  of  the  surficial  aquifer  [ft/d], 

L  =  length  from  center  of  the  island  to  the  coast  [ft], 

The  assumptions  for  discharge  along  a  seepage  face  are  that,  1)  the  flow  lines  are  horizontal 

and  equipotentials  are  vertical,  and  2)  the  hydraulic  gradient  is  assumed  to  be  equal  to  the 

slope  of  the  phreatic  surface  (Freeze  and  Cherry,  1979).  The  total  discharge  along  the 

seepage  face  was  calculated  by  assuming  a  phreatic  surface  of  zero  at  the  seepage  face  and 

subtracting  this  value  from  the  measured  water  level  at  Site  3.  Leakage  to  or  from  an 

overlying  or  underlying  aquifer  is  expressed  as: 

q=  K'Ah  (2) 

b' 
where 

K'  =  vertical  hydraulic  conductivity  of  confining  layer  [ft/d], 

Ah  =  difference  in  potentiometric  head  between  aquifers  [ft],  and 

b'  =  confining  bed  thickness  [ft]. 

The  surficial  aquifer  is  a  potential  source  of  recharge  for  the  underlying  aquifer,  and 
is  under  steady  state  conditions,  unconfined,  homogeneous,  anisotropic,  and  receives  areal 
recharge  from  precipitation.  The  partial  differential  equation  that  describes  the  three- 
dimensional  flow  in  the  surficial  aquifer  as  modified  for  aquifer  assumptions  and  boundary 
conditions  from  Freeze  and  Cherry  (1979)  is  expressed  as: 

d_  (Kx  (h  -  zb)  3hl  +  d_  (Ky  (h  -  zb)  3h)  +  d  (Kz  (h  -  zb)  3h)  +  Qprecipitation 
dx  dx        dy  dy      dz  dy 

(3) 

-  (K'h    Ah)  -  Qevapotranspiration     ~    [K    (hi2    -   hr2)!    =    0 

b  2L 

where 
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Zb  =  altitude  of  the  base  of  the  surficial  aquifer  [ft], 

Q  =  volumetric  inflow,  outflow  rate  [ft3/d]. 

The  altitude  of  the  bottom  of  the  aquifer  (zb)  was  subtracted  from  the  head  in  the  aquifer  to 

account  for  saturated  thickness  changes  from  the  center  of  the  island  to  the  coast 

The  Pliocene-Miocene  aquifer  consists  of  two  zones  hydraulically  connected  but 
geologically  distinct  At  Site  3,  the  sand  layer  is  separated  from  the  limestone  layer  by 
approximately  8  feet  of  sandy  clay.  A  difference  in  potentiometric  head  of  0.14  feet  exists 
between  the  sand  and  limestone  water  bearing  zones.  However,  if  the  two  water-bearing 
zones  are  considered  as  one  aquifer,  then  an  average  or  effective  transmissivity  for  both 
zones  may  be  assumed.  The  Pliocene-Miocene  aquifer  receives  approximately  8.5  x  104 
gal/d/mi2  recharge  from  the  overlying  surficial  aquifer.  Based  on  water  levels  and  the 
hydraulic  conductivity  of  the  zones,  approximately  6.2  x  104  gal/d/mi2  leak  from  the 
Pliocene  sand  to  the  Miocene  limestone. 

The  Pliocene-Miocene  aquifer  is  a  potential  sink  for  discharge  from  the  overlying 
surficial  aquifer,  under  steady  state  conditions,  confined,  heterogeneous,  and  anisotropic. 
The  equation  governing  the  flow  in  the  Pliocene-Miocene  aquifer  as  modified  for  aquifer 
assumptions  and  boundary  conditions  from  Freeze  and  Cherry  (1979)  is  expressed  as: 

a.(Txahi  +  a_(Tyah)  +  a_(Tzah)  +  ac,  Ah)  =  o  (4) 

dx        dx        By        dy        dz         dz  b 

where 

h  =  water  level  [ft], 

Ka'  =  vertical  hydraulic  conductivity  of  confining  layer  above  the  aquifer  [ft/d], 

b  =  thickness  of  aquifer  [ft], 

T  =  transmissivity  [ft2/d]. 

As  defined  in  the  equation,  the  aquifer  receives  leakage  from  the  overlying  surficial  aquifer.  The 

schematic  diagram  (Figure  32)  portrays  the  Pliocene-Miocene  aquifer  flow  as  described  by 

equation  (4). 
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Aquifer  test 

The  48-hour  aquifer  test  was  conducted  December  9-11,  1989,  at  Site  3  on 
Cumberland  Island.  Well  KBMP  7  was  used  as  the  pumping  well;  and  KB  MP  1 1  and 
KBMP  12  were  used  as  observation  wells  (see  Figure  33  for  location  of  wells).  In 
addition  to  these  wells,  KBMP  8,  KBMP  9,  and  KBMP  10  were  used  as  observation  wells 
in  the  overlying  and  underlying  zones  to  monitor  leakage  during  the  test.  Figure  34 
schematically  shows  the  hydrogeology  at  the  test  site  and  well  construction  of  the  pumping 
and  observation  wells. 

Prior  to  the  aquifer  test,  a  step-drawdown  test  was  performed  to  determine  the  well 
efficiency  and  the  optimum  pumping  rate  for  the  aquifer  test.  Based  on  data  analysis,  the 
optimum  pumping  rate  was  1.47  ft3/min  (11  gal/min).  The  pumping  rate  used  during  the 
test  was  1.54  ft3/min  (1 1.5  gal/min).  (See  Appendix  F  for  the  step  drawdown  procedure 
and  calculations  of  the  well  efficiency  and  optimum  pumping  rate.) 

The  aquifer  test  methods  of  analysis  were  selected  based  on  the  similarity  between 
the  theoretical  assumptions  of  the  methods  and  the  real  system  flow.  If  the  method  chosen 
involved  a  type  curve,  the  plot  of  the  data  matched  the  type  curve  shape  fairly  well.  Often 
the  flow  in  an  aquifer  involves  parameters  that  are  not  assumed  by  the  method,  such  as 
heterogeneities  or  boundaries.  In  these  cases,  the  shape  of  the  plotted  data  curve  indicated 
the  type  of  boundary  present  in  the  system,  such  as  a  recharge  or  impermeable  boundary. 
By  evaluating  the  hydrogeological  system,  the  assumptions  of  the  method,  and  the  curve 
shape  of  the  data,  the  direction  of  the  boundary  or  heterogeneity  (lateral  or  vertical)  was 
determined. 

The  following  section  discusses  the  methods  applied  for  aquifer  test  analysis.  The 
Hantush- Jacob  method  and  the  Streltsova  method  evaluate  for  vertical  heterogeneities  and 
the  Bixel-Van  Poollen  method  evaluates  for  lateral  discontinuities  in  the  aquifer. 
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EXPLANATION 


KB  MP  10 

133 


DATA  POINT-- 

Number  on  top  is  well  identification 
Number  on  bottom  is  total  depth  of  well,  in  ft. 

1"=  50" 


KBMP  11 
95 
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KBMP  7     8 

89    30 

KBMP  10 
133 
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Figure  33.  Location  and  depth  of  pumping  and  observation  wells  for  aquifer  test,  site  3. 
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Interpretation 

Nonsteady  radial  flow  in  a  confined  aquifer  is  expressed  mathematically  as: 

92h  +  J_9h  =  §_dh_  =  I2h  (5) 

ar2    t  dv    t  at     ri  at 

which  is  known  as  the  diffusion  equation,  where 
h  =  head  in  the  aquifer  [ft] 
r  =  radial  distance  from  the  pumping  well  [ft] 
S  =  storage  coefficient  [dimensionless] 
T  =  transmissivity  [ft2/d] 
t  =  time  [d] 

T[  =  diffusivity  coefficient  [ft2/d]  (Lohman,  1972). 

Appendix  G  provides  derivation  of  the  diffusion  equation  based  on  Darcy's  law  and  the 
principle  of  continuity  equation.  The  diffusion  equation  has  a  number  of  solutions  based 
on  the  boundary  conditions  of  the  aquifer.  The  hydraulic  properties  of  the  Pliocene- 
Miocene  aquifer  were  determined  by  methods  which  solve  the  diffusion  equation  by 
applying  boundary  conditions  to  the  equation. 

Theis  (1935)  applied  the  theory  of  heat  flow  to  arrive  at  an  analytical  solution  of  the 
diffusion  equation  that  involves  initial  assumptions  and  boundary  conditions  for  the  flow  in 
the  aquifer.  The  boundary  conditions  to  be  satisfied  are  as  follows  (Theis,  1935): 

h  (r,t) 

h(r,0)  =  h0  for    r>0 

h(«>,t)  =  h0  for  t  >  0 

The  Theis  equation  is  written  in  terms  of  drawdown: 
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s  = 


47tT 


where 

Q  =  discharge  rate  [ft3/d} 
and 


(6) 


u  =  r2_S_ 
4Tt 


(7) 


See  Figure  35  for  assumptions. 

The  exponential  integral  of  equation  (6)  is  known  as  the  well  function  W(u)  for  a 
nonleaky  aquifer  where: 

W(u) 


,.,.[*? 


(8) 


Values  of  W(u)  for  values  of  u  are  provided  in  tables  in  Lohman  (1972)  and  Freeze  and 
Cherry  (1979).  By  substituting  W(u)  in  equation  (6)  drawdown  may  be  expressed  as: 


or 


and 


or 


s  = 

0 

47tT 

W(u) 

logio  s  = 

=  [logio 

t  = 

Q  ] 

47tT 

=  r2_3_ 
4Tu 

+  logio 

W(u) 

logio  t  =  [logio  _r2_SJ  +  logio  1/u 
4T 


(9) 


(10) 


(11) 


(12) 


(Lohman,  1972). 

The  bracketed  parts  of  equations  (10)  and  (12)  are  constant  if  the  discharge  rate  is 
constant  during  the  aquifer  test.  If  values  of  drawdown  versus  time  are  plotted  on 
logarithmic  paper  and  superimposed  to  the  type  curve  (W(u)  versus  1/u),  a  fit  between  the 
data  curve  and  type  curve  are  made.  By  selecting  a  common  match  point  from  the  fit, 
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WW 

\\\  \ 

Direction  of  groundwater  flow 


log  W(u) 


log  1/u 


The  governing  equations  are: 

T=  _Q_      W(u) 
4rcs 


S  =  4Ttu 

r2 


The  Theis  assumptions  are: 

1 .  The  aquifer  is  confined; 

2.  nonsteady  flow  occurs; 

3.  the  aquifer  is  horizontal  having  an  infinite  aerial  extent; 

4.  the  pumped  well  fully  penetrates  the  aquifer  and  is  infinitesimal  in  diameter, 

5.  the  aquifer  is  homogeneous  and  isotropic; 

6.  flow  is  horizontal;  and 

7.  discharge  is  derived  exclusively  from  aquifer  storage. 


Figure  35.  The  Theis  governing  equations  and  assumptions. 
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values  of  W(u),  1/u,  s,  and  t  are  determined.  Solution  of  equation  (9)  is  possible  for  most 
instances  where  the  storage  capacity  of  the  producing  well  is  negligible  and  the  distance  to 
the  observation  well  is  equal  to  or  greater  than  200  times  the  pumping  well  radius 
(Streltsova,  1988). 

Hantush- Jacob  Method 

The  data  curve  from  the  aquifer  test  was  not  similar  to  the  Theis  type  curve  and  a 
match  could  not  be  made.  The  assumptions  of  the  Hantush- Jacob  Leaky  method  (1955) 
match  closely  to  the  real  system  and  the  shape  of  the  type  curve  matched  the  data  curve  (see 
Figure  36  for  assumptions).  Therefore,  this  analytical  method  was  used  to  analyze  the  test 
data. 

The  Hantush- Jacob  solution  of  the  diffusion  equation  is  as  follows: 

22h  +_i2h  -hK'  =  S_2JL  =  I2h  (13) 

dr2       r  3r     T  b'      T  dt         x\  dt 
where 

K'  =  vertical  hydraulic  conductivity  of  the  confining  layer  [ft/d] 

b'  =  confining  layer  thickness  [ft]  (Hantush  and  Jacob,  1955). 

The  governing  boundary  conditions  are: 

h  (r,t) 

h(r,0)  =  h0  for   r  >  0 

h(oo,t)  =  h0  for  t  >  0 

Q  =  0  f  or  t  <  0 

Q  is  constant  at  values  >  0  for  t  >  0. 

limr     9h_=  -  O  fort>0 

r->°   dr        2kT 
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Direction  of  groundwater  flow 


log  L(u,v) 


W(u) 


log  1/u 


The  governing  equations  are: 


T=     O    L(u.v) 
4ns 


S   =  4Ttu 

r2 


K'  =   4T  v2  Jbl 
r2 


The  Hantush-Jacob  leaky  assumptions  are: 

1.  The  aquifer  is  a  leaky  confined  aquifer, 

2.  the  aquifer  is  overlain  or  underlain  by  a  leaky  confining  bed  having  uniform  hydraulic  conductivity  and 
thickness; 

3.  the  pumping  well  maintains  a  constant  discharge; 

4.  the  aquifer  is  horizontal  having  infinite  aerial  extent; 

5.  flow  in  the  confined  aquifer  is  two  dimensional  and  radial  in  the  horizontal  plane;  flow  from  the 
confining  bed  is  vertical; 

6.  water  is  not  released  from  storage  in  the  confining  bed  during  pumping;  and 

7.  the  pumping  well  has  an  infinitesimal  diameter  and  fully  penetrates  the  aquifer 

8.  water  level  in  the  overlying  permeable  zone  remains  constant  (Bedinger  and  Reed,  1988). 


Figure  36.  The  Hantush-Jacob  Leaky  governing  equations  and  assumptions. 
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The  Hantush- Jacob  solution  of  the  diffusion  equation  is  similar  to  the  Theis  solution,  but 
the  second  part  of  the  equation  is  more  complicated.  The  equation  is  expressed  as: 


where 


s=4if[2Ko(2v)~|  yexp(-^->dy] 


L(u,v)  =  [2Ko(2v)-      1  exp  C^l)dy] 


(14) 


(15) 


L(u,v)  =  well  function  for  a  leaky  aquifer 
y  =  the  variable  of  integration 
u  =  r2S/4Tt 


2  Vb'T 


T 

If  the  right  side  of  equation  (14)  is  represented  as  the  leakance  function  of  u  and  v, 

L(u,v),  the  equation  may  be  written  as: 

s  =  _Q_L(u,v)  (16) 

4tcT 
and 

K>4Tv2  (17) 

b'  r2 

(Hantush  and  Jacob,  1955). 

Again,  values  of  drawdown  versus  time  are  plotted  on  logarithmic  paper  and 
superimposed  to  the  type  curve  (L(u,v)  versus  1/u),  and  a  fit  between  the  data  curve  and 
type  curve  is  made.  By  selecting  a  common  match  point  from  the  fit,  values  of  L(u,v),  1/u, 
s,  and  t  are  determined. 

After  about  100  minutes  from  the  starting  time  of  drawdown  and  recovery,  tidal 
fluctuations  had  a  greater  impact  on  the  water  levels  in  each  well  than  the  magnitude  of 
pumping  or  recharge.  The  data  were  corrected  for  the  influencing  tidal  and  barometric 
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pressure  fluctuations  before  they  were  matched  to  the  theoretical  curves.  (See  Appendix  I 
for  well  tidal  efficiency  and  the  method  applied  to  correct  for  tidal  fluctuations) 

Drawdown  and  recovery  data  from  wells  KB  MP  1 1  and  KBMP  12  provided  a 
good  fit  with  the  Hantush-Jacob  theoretical  curve  (see  Figures  37  to  40).  The 
transmissivity  values  calculated  by  the  Hantush-Jacob  method  ranged  from  235  to  504  ft2/d 
(Table  4).  (Aquifer  test  calculations  and  data  tables  are  provided  in  Appendix  H.)  The 
storage  coefficient  may  be  solved  for  by  rearranging  equation  (7).  The  storage  coefficient 
values  ranged  from  1.4  x  10"5  to  5.6  x  10"5.  Once  the  transmissivity  and  storage 
coefficient  are  known  the  diffusivity  may  be  calculated  by: 


•n  =  x 
s 


(18) 


r\  =  diffusivity  of  the  aquifer  [ft2/d] 

(Freeze  and  Cherry,  1979).  The  diffusivity  values  ranged  from  4.5  x  106  to  3.5  x  107  ft2/d 

(see  Table  4). 

The  delayed  drawdown  and  recovery  response  in  KBMP  12  for  the  first  3  minutes 
after  pumping  began  may  be  due  to  poor  lateral  fracture  interconnection  or  facies  changes 
between  KBMP  12  and  the  pumping  well.  In  addition  to  delayed  drawdown  and  recovery, 
the  total  drawdown  in  well  KBMP  12  was  less  than  in  well  KBMP  1 1  which  is  located 
further  from  the  producing  well.  This  indicates  that  less  interconnected  fractures  exist 
between  KBMP  12  and  the  pumping  well. 

The  Hantush-Jacob  method  was  designed  for  hydraulic  analysis  of  leaky  aquifers. 
The  calculated  storage  values  indicate  that  the  aquifer  is  well  confined,  and  no  leakage  was 
observed  from  the  overlying  and  underlying  zone  during  the  test.  The  aquifer  responds  as 
if  it  were  a  leaky  aquifer  because  it  is  composed  of  two  lithologic  zones  with  the  overlying 
zone  having  a  much  greater  transmissivity  than  the  underlying  zone.  As  water  is 
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discharged  from  the  underlying  limestone  unit,  the  overlying  sand  provides  recharge  to  the 
limestone.  This  is  reflected  in  the  plotted  data  as  the  curve  flattens  late  in  the  test. 

Straight  line  method 

Streltsova  Method 

Streltsova  (1988)  presented  a  method  that  numerically  solves  for  the  transmissivity 
in  a  stratified  aquifer  during  transient  groundwater  flow.  The  model  evaluates  for  the  effect 
that  a  permeable  layer  separated  by  an  impervious  zone  has  on  the  underlying  producing 
zone.  (Figure  41  illustrates  the  direction  of  groundwater  flow  in  a  stratified  aquifer  during 
discharge.)  Change  in  drawdown  with  time  for  a  well  producing  at  a  constant  rate  from  a 
stratified  aquifer  results  in  multiple  slopes  on  a  semilog  graph.  Early  drawdown  data  on  a 
semilog  graph  plots  a  straight  line  slope  that  represents  the  transmissivity  of  the  producing 
layer.  The  transmissivity  of  this  zone  is  calcuated  by: 

Ti=    2.3  0  (19) 

4rc  mi 

where 

Ti  =  transmissivity  of  the  producing  zone  [ft2/d] 

mi  =  slope  of  the  initial  segment  [ft/min]  (see  Figure  41  for  assumptions). 

When  the  semipervious  layer  begins  to  influence  system  flow,  the  curve  grades  into  a 

transitional  slope,  the  shape  and  duration  of  which  depend  on  the  storage  ratio  of  the  two 

regions  (S'/S  where  S'  is  the  storage  of  the  overlying  permeable  zone).  When  the  second 

permeable  layer  begins  to  contribute  to  the  groundwater  flow,  the  slope  of  the  final  curve 

(m')  is  dependent  upon  the  transmissivity  of  both  permeable  layers  (Ti  and  T2).  The  final 

segment  slope  is  related  to  the  initial  slope  by  the  relationship: 


> 

- T2 

\t\l\ 

\  Kt\ 

|       ;           T, 

Direction  of  groundwater  flow 


log  time 


The  governing  equations  are: 
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Ti=2JJ} 
4ran 


T->  =  Ti  (mi  -m') 
m' 


The  Streltsova  assumptions  are: 

1.  stratified  aquifer 

2.  vertical  flow  in  the  semipervious  layer 

3.  horizontal  flow  in  the  highly  transmissive  layers 
Each  layer  is: 

4.  homogeneous,  isotropic 

5.  compressible 

6.  of  infinite  radial  extent 


Figure  41.  Streltsova  (1988)  governing  equations  and  assumptions. 
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m"=        miTi (20) 

Ti  +  T2 

where 

m'  =  slope  value  of  the  final  segment  [ft/min] 

T2  =  transmissivity  of  the  permeable  zone  overlying  the  producing  zone  [ft2/d]. 

The  transmissivity  of  each  layer  is  provided  in  Table  4.  As  the  data  indicate,  the 
transmissivity  is  usually  one  and  one-half  times  greater  when  the  overlying  sand  zone 
contributes  to  the  groundwater  flow.  The  transmissivity  values  calculated  by  the  Streltsova 
straight  line  method  range  from  380  to  792  ft2/d.  A  transition  curve  does  not  exist  in  the 
test  data  (Figures  42  through  45).  This  indicates  that  the  storage  of  the  producing  and 
overlying  permeable  zones  are  similar. 

Bixel-Van  Poollen 

The  change  in  drawdown  with  time  for  a  well  producing  at  a  constant  rate  from  a 
formation  with  radial  discontinuity  is  represented  by  multiple  slopes  on  a  semilog  plot.  A 
numerical  evaluation  of  transient  drawdown  for  a  well  centered  in  an  aquifer  with  radial 
discontinuities  was  made  by  Bixel  and  Van  Poollen  (1967).  See  Figure  46  for  assumptions 
and  schematic  diagram  of  an  aquifer  with  radial  discontinuities.  Early  drawdown  data  on  a 
semilog  graph  plots  a  straight  line  of  a  slope  that  represents  the  transmissivity  of  the  aquifer 
region  from  the  producing  well  to  the  radial  discontinuity.  The  transmissivity  of  this  region 
in  the  aquifer  is  calculated  by  equation  (19).  When  the  region  outside  the  discontinuity 
begins  to  influence  the  drawdown,  the  straight  line  grades  into  a  transitional  curve,  the 
shape  and  duration  of  which  depend  on  the  storage  ratio  of  the  two  regions.  Finally,  when 
the  groundwater  is  produced  from  the  region  outside  the  discontinuity,  the  drawdown 
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The  governing  equations  are: 


Ti    =  2.3  0              S    = 

T 

4Jlmi 

11 

T?=Ti  (2mi  -m') 

m' 

T_         n    = 


r2 


2.246  (t0) 


The  Bixel  and  Van  Poollen  assumptions  are: 

1.  Fluids  in  the  aquifer  are  slightly  compressible; 

2.  permeability  and  porosity  are  independent  of  pressure; 

3.  homogeneous,  isotropic  medium; 

4.  constant  aquifer  thickness; 

These  parameters  are  constant  but  may  be  different  on  opposite  sides  of  the  discontinuity 

5.  constant  pumping  rate; 

6.  continuous  flow  across  the  discontinuity; 

7.  the  aquifer  is  of  infinite  areal  extent. 


Figure  46.  Bixel  and  Van  Poollen  (1967)  governing  equations  and  assumptions. 
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curve  exhibits  a  final  straight  line  (m').  The  final  slope  that  represents  the  outside  region, 
m',  is  related  to  the  initial  slope,  mi,  by  the  relationship  of: 

m'  =  mi( IL_-> 

l+(n-l>^ 

Ti  (21) 

where 

T2  =  transmissivity  of  the  aquifer  outside  the  radial  discontinuity  [ft2/d]. 

(n-1)  is  the  number  of  image  wells.  An  image  well  is  equivalent  to  a  fictitious  recharge  or 

discharge  straight  line  boundary  (Bear,  1979).  If  a  linear  recharge  boundary  is  assumed, 

then  (n-1)  =  1  or  n  =  2.  Substituting  (n  =  2)  into  equation  (21)  results  in: 

m'  =  mi<— ^— > 

Ti  (22) 

The  transmissivity  values  calculated  by  the  Bixel-Van  Poollen  straight  line  method 

range  from  380  to  2,021  ft2/d  (Table  5).  Once  the  transmissivity  and  diffusivity  of  the 

aquifer  are  known,  the  storage  coefficient  may  be  determined  by  rearranging  equation  (7). 

The  diffusivity  of  the  aquifer  may  be  calculated  by  an  equation  derived  by  Streltsova  (1988) 

that  is  expressed  as: 

ti  = t2 (23) 

2.246  (to) 

where 

to=  intercept  of  straight  line  at  zero  drawdown  [rnin.].  The  diffusivity  values  calculated  by 

equation  (23)  range  from  5.98  x  103  to  4.31  x  104  ft2/d.     The  storage  values  range  from 

1.1  x  10"5  to  4.4  x  10'5.  Based  on  aquifer  testing,  the  hydraulic  conductivity  for  the 

Miocene  limestone  ranges  from  34  to  101  ft/d.  The  hydraulic  conductivity  of  the  Pliocene 

sand  is  100  ft/d  based  on  the  method  of  deriving  this  value  from  grain  size  analysis  as 

described  by  Masch  and  Denny  (1966). 
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Transmissivity  values  calculated  for  well  KBMP  7  were  slightly  higher  than  those  values 
obtained  from  the  observation  wells  (see  Figures  47  and  48).  These  values  were  calculated 
from  equation  (19).  Data  recorded  after  100  minutes  were  not  included  in  the  straight-line 
fit  because  the  water  level  began  to  fluctuate  with  changing  tides.  The  first  4  minutes  of 
drawdown  were  also  excluded  in  the  analysis  because  of  the  casing  storage  effect. 

Discussion 

The  transmissivity  of  an  aquifer  is  directly  related  to  the  hydraulic  conductivity  and 
the  saturated  thickness  of  the  aquifer  medium;  the  thicker  the  saturated  zone  and  more 
interconnected  the  pore  spaces,  the  greater  the  transmissivity  of  the  aquifer.  The  storage 
coefficient  of  an  aquifer  is  directly  related  to  the  porosity  of  the  medium,  and  the 
compressibility  of  the  rock  matrix  and  fluid.  The  greater  these  values,  the  greater  the 
storage  coefficient. 

For  a  homogeneous,  isotropic  aquifer,  aquifer  test  data  should  plot  a  straight  line  on 
a  semilog  graph.  Breaks  in  the  straight  line  graph  reflect  vertical  heterogeneities  or  lateral 
discontinuities  in  the  aquifer  medium.  A  three  segment  curve,  excluding  early  borehole 
effects,  indicates  variations  of  transmissivity  and  storage  values  in  the  aquifer  (see  Figure 
49).  The  first  segment  represents  the  transmissivity  inside  the  discontinuity  in  an  aquifer 
with  lateral  discontinuities.  In  a  stratified  aquifer,  this  segment  represents  the 
transmissivity  of  the  producing  zone.  The  second  segment  represents  the  transition  of 
transmissivity  in  the  medium  from  the  region  inside  the  discontinuity  to  the  region  outside 
the  discontinuity,  or  in  a  stratified  aquifer,  the  transition  from  the  producing  zone  to  the 
overlying  permeable  zone.  The  shape  and  duration  of  the  second  segment  depends  on  the 
storage  ratio  between  the  two  regions.  If  Si  =  S2  =  S  then  a  transient  curve  will  not  exist 
on  a  semilog  plot.  The  third  segment  represents  the  transmissivity  of  the  region  outside  the 
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If  Si  =  s2 
and  Ti  *  10  x  T2 


T  =  transmissivity  [ft2/d] 
S  =  storage  coefficient 
s  =  drawdown  [ft] 


log  time 


If  Si  *  10  x  S2 
and  Ti  *  10  x  T2 


log  time 


Figure  49.  Configuration  of  aquifer  test  plots  with  varying  transmissivity  and  storage  values. 
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lateral  or  vertical  discontinuity.  In  a  stratified  aquifer,  the  third  segment  represents  the 
transmisssivity  of  a  combination  of  the  overlying  permeable  zone  and  the  producing  zone. 
In  a  stratified  aquifer,  a  three  segment  curve  indicates  that  the  transmissivity  and  storage  of 
the  producing  zone  and  the  overlying  permeable  zone  are  not  equal  (Streltsova,  1988). 

A  single  break  in  the  straight  line  or  a  two  segment  curve  on  a  semilog  plot, 
excluding  early  borehole  effects,  indicates  that  the  aquifer  permeability  and/or  thickness 
vary  across  the  discontinuity.  The  transmissivity  inside  the  discontinuity  (Ti)  and  the 
transmissivity  outside  the  discontinuity  (T2)  are  not  equal  (Figure  49).  It  may  be  assumed 
that  the  storage  value  is  uniform  in  the  aquifer  because  a  transition  curve  does  not  exist  (the 
storage  ratio  equals  one). 

Aquifer  test  data  from  the  Pliocene-Miocene  aquifer  yields  a  3  segment  curve  on  a 
semilog  graph  (see  Figures  42  to  45).  The  first  slope  represents  borehole  and  near 
borehole  effects.  The  second  segment  represents  the  transmissivity  of  the  limestone  part  of 
the  aquifer,  and  the  third  slope  represents  contribution  from  the  sand  unit  of  the  aquifer 
and/or  a  lateral  recharge  boundary.  A  transition  curve  is  not  present  which  indicates  that 
the  storage  coefficient  is  uniform  in  the  aquifer.  Based  on  analysis  of  the  Streltsova  and 
Bixel-Van  Poollen  methods,  this  curve  shape  indicates  differences  of  transmissivity 
between  the  sand  and  limestone  zones  of  the  aquifer  and  across  the  lateral  discontinuity. 

The  first  segment  of  early  time  data  represents  near  well  bore  conditions  that 
includes  factors  such  as  well  bore  storage,  skin  effect,  formation  damage,  and  near  well 
bore  formation  transmissivity.  The  length  of  time  for  the  well  storage  effect  to  subside  was 
about  4.3  minutes  after  initial  pumping  time  (see  Appendix  F  for  calculations). 

The  second  segment  represents  the  transmissivity  of  the  limestone  part  of  the 
aquifer  or  the  producing  zone  inside  the  lateral  discontinuity.  This  is  assumed  because  the 
limestone  has  a  much  lower  transmissivity  than  the  overlying  sand  and  should  therefore 
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have  a  greater  slope  value  than  a  slope  represented  by  the  combination  of  the  two  water- 
bearing zones  and  the  zone  outside  the  lateral  discontinuity.  The  transmissivity  of  this  zone 
was  calculated  from  equation  (19).  See  Table  4  for  values  calculated  by  each  straight  line 
method. 

The  third  segment  represents  vertical  groundwater  contribution  from  the  overlying 
sand  zone  and/or  a  lateral  recharge  boundary.  The  recharge  effects  are  represented  by  the 
low  slope  value  of  the  third  segment.  The  third  segment  indicates  that  the  magnitude  of 
recharge  is  greater  than  the  magnitude  of  drawdown  and  represents  contribution  from  both 
water-bearing  zones.  The  transmissivity  values  calculated  for  the  third  slope  were  2  to  4 
times  greater  than  the  transmissivity  calculated  from  the  slope  of  the  second  segment  (See 
Table  5  for  transmissivity  values.) 

Vertical  heterogeneities  and  lateral  discontinuities  are  present  in  the  Pliocene- 
Miocene  aquifer  as  observed  in  the  aquifer  test  data.  As  a  result  of  these  heterogeneities, 
the  values  of  transmissivity  vary  vertically  and  laterally.  Vertical  heterogeneity  in  the 
aquifer  medium  is  due  to  the  overlying  Pliocene  sand  having  a  much  greater  transmissivity 
value  than  the  Miocene  limestone.  As  the  overlying  sand  zone  contributes  to  the  aquifer 
flow,  the  transmissivity  values  increase  one  and  one-half  times  (see  Table  4).  The  lateral 
discontinuity  could  be  the  result  of  thickness  or  facies  changes  in  the  aquifer  or  possibly 
recharge  from  the  ocean. 

The  Streltsova  (1988)  model  solves  for  the  transmissivity  of  each  permeable  zone 
of  a  vertical  section  in  the  stratified  aquifer.  Application  of  this  method  provided  the 
transmissivity  values  of  the  limestone  (Ti)  and  the  combination  of  the  limestone  and  sand 
units  (T2)  of  the  Pliocene-Miocene  aquifer.  Values  indicate  that  the  combination  of  the  two 
zones  exhibits  a  higher  transmissivity  value  than  the  limestone  zone  alone. 
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The  Bixel-Van  Poollen  model  solves  for  transmissivity  in  an  aquifer  having  radial 
discontinuities.  Application  of  this  method  with  the  assumption  of  one  lateral  recharge 
boundary  provided  the  transmissivity  of  the  zone  inside  the  discontinuity  (Ti)  and  the  zone 
outside  the  discontinuity  (T2).  Values  indicated  that  the  formation  outside  the  discontinuity 
had  a  greater  transmissivity  than  the  formation  inside  the  discontinuity.  Sufficient  data  was 
not  available  to  determine  the  lateral  source  of  recharge  in  the  aquifer. 

The  transmissivity  calculated  by  the  straight  line  methods  ranged  from  380  to  653 
ft2/d  for  the  Ti  values.  The  Ti  values  are  considered  the  valid  values  for  the  aquifer 
because  the  pumping  well  penetrated  only  the  limestone  unit  of  the  aquifer.  Transmissivity 
values  calculated  for  the  Pliocene-Miocene  aquifer  in  Kingsland,  Georgia,  15  miles  west  of 
the  study  area,  were  750  ft2/d  based  on  aquifer  tests  (Brown,  1984). 

Potential  for  Seawater  Encroachment  due  to  Channel  Dredging 

An  evaluation  of  the  volume  of  seawater  infiltration  was  conducted  along  the 

Cumberland  Sound  Channel  because  the  confining  unit  overlying  the  Pliocene-Miocene 

aquifer  is  breached  in  some  areas  below  the  channel.  For  this  discussion,  it  is  assumed  that 

the  confining  layer  overlying  the  Pliocene-Miocene  age  aquifer  was  completely  removed  in 

the  channel,  and  the  aquifer  is  in  direct  contact  with  the  channel  bed.  Following  Darcy's 

Law  (Freeze  and  Cherry,  1979),  specific  discharge  is  expressed  as: 

q  =  -K  Ah_  (24) 

AL 

where 

q  =  specific  discharge  [ft/hr] 

K  =  hydraulic  conductivity  [ft/hr] 

h  =  water  level  [ft] 

L  =  distance  [ft]. 
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The  difference  in  hydraulic  head  of  the  channel  and  the  Pliocene-Miocene  aquifer  fluctuates 

with  time  as  a  result  of  the  fluctuation  of  the  tides  (Figure  31).  Therefore  the  equation  is 

time  dependent  and  can  be  expressed  as  a  function  of  time: 

q(t)=  -KAh_(t)  (25) 

AL 

where 

t  =  time  [hr] 

The  change  in  head  difference  with  respect  to  time  can  be  approximated  by  a  series 

of  straight  lines  as  shown  in  Figure  50.  The  change  in  head  may  be  expressed  as: 


dh  =  m 
dt 


(26) 


The  first  slope  may  be  expressed  by  the  integral: 


j   dh=  J    mdt 

./ho  JO 


By  integrating  over  a  24-hour  period  the  equation  is  expressed: 

Ahj  =  m;t  +  bi 

i  =  1,5 


(27) 


(28) 


m  =  slope  of  the  line  [ft/hr] 

b  =  y-intercept  [ft] 

i  =  number  of  slope  intervals 

During  a  24-hr  period  in  July  1989,  the  head  in  the  aquifer  was  higher  than  channel  stage 

during  two  low-tide  cycles  and  one  high-tide  cycle  (Figure  50).  However,  water  levels 

within  the  channel  were  greater  than  the  aquifer  during  alternating  high  tides.  Therefore, 

during  83  percent  of  the  day  (20  hrs)  the  potential  for  groundwater  flow  was  toward  the 

channel,  and  during  17  percent  of  the  day  (4  hrs),  the  potential  for  flow  was  from  the 

channel  into  the  aquifer. 
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The  volume  of  water  discharged  per  unit  area  can  be  calculated  by  substituting  equation 
(28)  into  equation  (25)  and  integrating.  During  a  24-hour  period  specific  discharge  can  be 
expressed  by  the  integral: 


j     qdt  =  JM 

Jo  alJ0 


24 

(mt  +  b)dt 

(29) 


Integrating  the  equation  yields  the  volume  per  unit  area: 

qt|204=J^<im2+bt>2o4 

AL     2  (30) 

The  volume  of  water  discharged  per  square  foot  during  a  24-hr  period  (1  day)  is  expressed 

as: 

^ = *-  ie£ + b,  $  *  ef- + ^i + eg + brf? 

A     AL        2  2  2 

+  ef  +  b4t)f5  +  ef  +  b4t>l>  +  (Sfi  +  b5t>H]         (3  X) 

When  the  hydraulic  gradient  was  sloped  from  the  Pliocene-Miocene  aquifer  toward 
the  channel  (20  hrs),  the  volume  of  water  discharged  from  the  aquifer  per  unit  area  was 
0.02  ft3/ft2  (0.15  gal/ft2).  When  the  hydraulic  gradient  was  sloped  from  the  channel 
toward  the  Pliocene-Miocene  aquifer  (4  hrs),  the  volume  of  seawater  infiltrating  the  aquifer 
was  0.001  ft3/ft2  (0.007  gal/ft2).  The  net  volume  of  freshwater  discharged  from  the 
Pliocene-Miocene  aquifer  per  unit  area  was  0.019  ft3/ft2  d  (0.142  gal/ft2  d).  These  values 
are  small  and  indicate  that  little  movement  occurs  between  the  aquifer  and  the  channel. 
Seawater  intrusion  from  the  channel  was  insignificant  near  Site  3  because  the  total  volume 
of  freshwater  discharging  from  the  aquifer  was  an  order  of  magnitude  greater  than  the  total 
volume  of  seawater  entering  the  aquifer.  The  seepage  velocity  in  the  aquifer  is  0.064  ft/d 
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which  indicates  groundwater  movement  in  the  aquifer  is  slow.  Any  seawater  that  intrudes 
the  aquifer  will  not  travel  far. 

The  volume  of  water,  inflow  and  outflow,  could  not  be  estimated  at  Sites  1  and  2 
due  to  lack  of  water-level  data.  However,  the  high  chloride  concentrations  reported  in  the 
aquifer  at  the  sites  (ranging  from  13,000  to  17,000  mg/L)  indicate  that  seawater 
encroachment  has  occurred  at  the  southern  end  of  the  island  (Wilson,  1990).  Because  there 
was  no  water-quality  data  for  the  Pliocene-Miocene  aquifer  prior  to  dredging,  it  may  not  be 
possible  to  conclude  that  channel  dredging  promoted  seawater  intrusion. 
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CONCLUSIONS 

The  geology  of  the  study  area  consists  of  alternating  layers  of  sand,  clay, 
limestone,  and  dolomite.  The  upper  Charlton  Formation  consists  of  consolidated  to  semi- 
consolidated  calcareous  sediment.  An  open  marine  environment  existed  during  deposition 
as  indicated  by  the  abundant  marine  fossils  and  burrowing  found  in  the  core.  The 
overlying  Duplin  Formation,  ranging  from  2  to  42  feet  thick,  consists  of  coarse  sand,  shell 
fragments  and  clay  lenses.  The  formation  increases  in  clay  content  at  the  southern  most 
well  sites.  The  Holocene- Pleistocene  undifferentiated  deposits  consist  of  a  basal  clay  layer, 
5  to  40  feet  thick,  overlain  by  silty  sand  ranging  from  22  to  62  feet  thick. 

The  Miocene  sand  aquifer  consists  of  coarse  to  very  coarse  calcareous  sand,  poorly 
sorted  with  lenses  of  argillaceous  dolomite.  The  aquifer  is  recharged  by  precipitation  in 
outcrop  areas  outside  the  study  area.  The  leakage  rate  from  the  Miocene  sand  aquifer  to  the 
Pliocene-Miocene  aquifer  was  3.1  x  10"6  ft/d.  The  hydraulic  conductivity  derived  from 
grain  size  analysis  from  well  KB  MP  1  was  53  ft/d. 

At  sites  1  and  2,  the  Pliocene-Miocene  aquifer  consists  of  fractured  dolomitic 
limestone.  At  site  3  the  aquifer  consists  of  two  hydraulically  connected  zones  that  are 
geologically  distinct:  20  to  30  feet  of  Pliocene-age  sand  underlain  by  about  10  feet  of 
Miocene-age  fractured  limestone.  These  two  zones  are  separated  by  8  feet  of  sandy  clay. 

The  horizontal  gradient  of  the  Pliocene-Miocene  aquifer  on  the  south  end  of 
Cumberland  Island  was  toward  the  west  based  on  June  28,  1990  water  levels.  The 
seepage  velocity  ranged  from  0.049  to  0.085  ft/d  and  changed  with  fluctuating  tides  and 
storm  surges.  The  vertical  gradient  varied  with  time  and  space.  At  well  sites  1  and  2  the 
direction  of  groundwater  flow  was  from  the  Pliocene-Miocene  aquifer  to  the  surficial 
aquifer,  at  well  site  3  the  direction  of  flow  was  reverse  from  the  surficial  aquifer  to  the 
Pliocene-Miocene  aquifer.  The  magnitude  and  direction  of  the  vertical  gradient  changed 
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with  fluctuating  tides  and  storm  surges. 

Based  on  aquifer  testing,  the  transmissivity  of  the  Pliocene- Miocene  aquifer  ranged 
from  235  to  650  ft2/d  and  the  storage  coefficient  from  1.05  x  10"5  to  5.6  x  10"5. 
Transmissivity  values  calculated  for  the  Pliocene-Miocene  aquifer  in  Kingsland,  Georgia, 
15  miles  west  of  the  study  area,  were  750  ft2/d  based  on  aquifer  tests  (Brown,  1984). 
Transmissivity  values,  storage  values,  and  the  straight  line  plots,  indicate  that  the  porosity 
of  the  aquifer  medium  is  uniform  but  the  hydraulic  conductivity  varies  throughout  the 
medium.  Variation  of  the  hydraulic  conductivity  is  the  result  of  lateral  and  vertical 
heterogeneities  in  the  medium.  Lateral  heterogeneities  may  include  fades  changes,  changes 
in  the  aquifer  thickness,  and  possibly  recharge  from  the  ocean.  Vertical  heterogeneities  are 
due  to  the  sand  unit  of  the  aquifer  having  a  much  higher  transmissivity  than  the  limestone 
unit. 

The  surficial  aquifer  is  composed  of  sands  and  clays  of  the  Holocene  and 
Pleistocene  ages.  The  aquifer  ranged  from  15  to  35  feet  of  saturated  thickness  and  was 
confined  below  by  5  to  40  feet  of  clay.  The  horizontal  direction  of  groundwater  flow  was 
toward  the  southwest  based  on  June  28,  1990  water  levels.    The  vertical  leakage  rate  from 
the  surficial  aquifer  to  the  Pliocene-Miocene  aquifer  at  site  3  was  4.08  x  1(H  ft/d.  At  sites 
1  and  2  the  vertical  gradient  between  these  aquifers  was  the  reverse.  The  calculated  leakage 
values  from  the  Pliocene-Miocene  aquifer  to  the  surficial  aquifer  was  2.09  x  10~5  ft/d  at  site 
1,  and  3.17  x  10"4  ft/d  at  site  2.  The  hydraulic  conductivity  of  the  surficial  aquifer  was 
about  21  ft/d  based  on  grain  size  analysis  of  samples  from  three  sites. 

The  surficial  aquifer  was  in  a  steady-state  condition  with  constant  long-term 
recharge,  equivalent  discharge  and  no  appreciable  withdrawals  due  to  pumping. 
Approximately  2.85  x  105  gallons  per  square  mile  recharged  the  aquifer  each  day  and 
nearly  the  same  amount  was  discharged  to  the  low  lying  wetlands  (1.92  x  105  gal/d),  the 
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underlying  Pliocene-Miocene  aquifer  (8.5  x  104  gal/d),  and  along  the  coast  at  the  seepage 
face  (7.58  x  103  gal/d). 

In  the  channel  adjacent  to  site  3,  the  volume  of  water  discharged  from  the  Pliocene- 
Miocene  aquifer  into  the  channel  was  an  order  of  magnitude  greater  than  the  total  volume  of 
seawater  that  entered  the  aquifer  during  a  24-hour  period  in  July  1989.  The  net  volume  of 
freshwater  discharged  from  the  aquifer  per  unit  area  was  0.019  ft3/ft2  d  (0.142  gal/ft2). 
Therefore,  seawater  encroachment  from  the  channel  was  insignificant  near  site  3.  Near 
sites  1  and  2,  however,  high  chloride  concentrations  in  the  aquifer,  ranging  from  13,000  to 
17,000  mg/L,  indicated  that  seawater  encroachment  has  occurred.  Because  water- 
chemistry  data  does  not  exist  for  the  Pliocene-Miocene  aquifer  prior  to  dredging,  it  is  not 
possible  to  conclude  that  channel  dredging  promoted  seawater  encroachment  into  the 
aquifer. 
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FUTURE  WORK 

Much  work  could  be  done  to  enhance  the  hydrogeologic  description  of  Cumberland 
Island.  The  feasible  possibilities  are  discussed  below. 

1.  Install  continuous  water-level  recorders  on  surficial-aquifer  wells  and  on  ponds  on  the 
island's  south  end  to  determine  the  extent  of  surface  and  ground  water  interaction. 
Continuous  water-level  recorder  data  along  with  data  from  weather  stations  (to  be  installed) 
will  yield  information  that  can  be  used  to  determine  an  annual  surficial-aquifer  water 
budget. 

2.  Install  continuous  water-level  recorders  on  the  Pliocene-Miocene  aquifer  wells  and  the 
Miocene  sand  aquifer  well  to  determine  vertical  and  horizontal  gradients  of  the  shallow 
aquifers  during  tidal  fluctuations  and  storm  surges. 

3.  Conduct  aquifer  tests  at  sites  1  and  2  for  the  Pliocene-Miocene  aquifer.  If  pumping  is 
not  feasible,  slug  tests  could  be  preformed.  Slug  tests  could  also  be  conducted  at  each  well 
site  for  the  surficial  aquifer.  Shallow  piezometers  or  well  points  could  be  installed  near  the 
control  well  to  monitor  water-level  changes  during  the  test 

4.  Install  3  or  more  monitoring  wells/piezometers  in  the  Pliocene-age  sand  to  determine  the 
horizontal  groundwater  flow  gradient  and  seepage  velocity  of  this  zone. 
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Table  Al.  Raw  grain  size  analysis  data 


Biufcuiiud 

Sieve  dumeiex  (e>        Beaker  weigtal(g)         unfair  wagti  (g) 

Sample  weight  (g) 

Cosseted  wogtii 

Co 

CummlMive 
paceol 

Sue  1.  13*M6  ft. 

interval 

Percent  error  = 
0.25% 

Sue  2,  20-30  ft. 
interval 

Percent  error  = 
031  * 

Site  3. 10-20  ft. 

interval 

Percent  error  = 
0.05% 

Site  3.  60-70  ft. 

interval 

Pexcem  error  = 
05* 

Sue  1.  10-11  5ft 
iotovil 

Pacax  ujui  = 
Oil* 

47.64 

8331 

3561 

-1 

48.07 

0.43 

121 

121 

0 

3145 

5.81 

1636 

1757 

1 

602 

1256 

3536 

5293 

2 

61.06 

1342 

37.78 

9071 

3 

49.41 

1.77 

458 

95.69 

4 

48.19 

055 

154 

9723 

P" 

48.62 

0.98 

2.76 

9959 

50.1 

85.61 

3551 

•1 

50.14 

0.04 

0.113 

0.113 

0 

5057 

0.47 

133 

1.44 

1 

51.69 

159 

4.49 

5.93 

2 

5759 

7.49 

21.16 

27.09 

3 

68.68 

1858 

5249 

7958 

4 

57.2 

7.1 

20.06 

99.64 

P"> 

5023 

0.13 

0367 

100 

6438 

108  65 

4427 

0 

64.42 

004 

009 

0.09 

1 

64.41 

0.068 

0.068 

0.158 

15 

64.66 

0.633 

0.633 

0.791 

2 

65.75 

3.1 

31 

3.89 

3 

98.69 

7754 

7754 

81.43 

4 

7252 

184 

18.4 

99.83 

P« 

64  46 

0.181 

0.181 

100 

502 

93.89 

43.69 

-2 

5054 

034 

0.78 

0.78 

-15 

52.17 

157 

453 

531 

-1 

53.86 

3.66 

842 

13.73 

0 

6959 

19.79 

4555 

5928 

1 

6037 

10.17 

2341 

8269 

2 

54  68 

4.48 

1031 

93 

3 

5253 

2.73 

628 

9928 

4 

5049 

029 

0.66 

9954 

jm 

5022 

0.02 

0.046 

9959 

50.11 

98.11 

48 

-l 

50.12 

0.01 

0.02 

0.02 

0 

5018 

0.07 

0.147 

0.167 

i 

5021 

0.1 

021 

0377 

u 

5075 

0.64 

134 

1.72 

2 

54.69 

458 

9.62 

1134 

3 

85.44 

3533 

7421 

8555 

4 

5666 

655 

13.76 

9931 

P« 

5044 

033 

0.69 

100 

119 


APPENDIX  B 
Core  Descriptions 


120 


Description  of  split  spoon  samples  and  core  samples  from  well  KB  MP  1 ,  Site  1 . 

Driller:  Corps  of  Engineers,  Savannah,  GA 
Latitude/Longitude:  30°43'11",  81028'13" 
Altitude:  4.59  ft.  above  mean  low  water 
Total  Depth:  155.5  ft. 

(Colors  were  wet  and  matched  to  the  Rock  Color  Chart,  GSA,  1984) 

Slit  spoon  samples 


Depth  in 
ft,  below 
land  surface 

Holocene  to  Pleistocene 
Undifferentiated 


Description 


5.0-6.5 

10-11.5 
15-16.5 

20-21.5 

23.0-23.5 

25-26.5 

30-31.5 


35-36.5 


Sand:  Yellowish  gray  (5Y7/2),  v.f.g.  sand,  subangular, 
well  sorted,  20%  heavy  minerals,  20%  organics. 

Sand:  Same  as  previous  section 

Sand:  Yellowish  gray  (5Y7/2),  v.f.g.  sand,  subangular, 
well  sorted,  20%  heavy  minerals,  20%  organics,  trace  of 
shell  fragments. 

Sand:  Olive  gray  (5Y4/1),  v.f.g.,  subangular,  20%  heavy 
minerals,  10  clay,  1%  mica  flakes. 

Clay:  Dk.  gray  (N3)  fat,  clean,  trace  of  shell  fragments, 
trace  of  mica 

Sandy  Clay:  Med.  gray  (N5),  v.f.g.  sand,  subangular,  alt. 
layers  of  fat  gray  clay  2-3  inch,  in  thickness  (30%),  slightly 
calcareous. 

Sandy  Clay:  Med.  gray  (N5),  v.f.g.  sand,  subangular,  alt. 
layers  of  fat  gray  clay  2-3  inch,  in  thickness  (30%),  10% 
organic  material,  slighdy  calcareous,  trace  of  shell 
fragments. 

Same  as  30-31.5  section 


Pliocene  Series 
40-41.5 


Clayey  Sand:  Brown  gray  (5YR4/1),  v.f.-f.g.,  subangular, 
20%  shell  fragments,  10%  clay,  slightly  calcareous. 
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45-46.5  Clayey  Sand:  Mod.  olive  brown  (5Y4/4),  f.g.-m.g.  sand, 

subangular,  15%  clay,  10%  heavy  minerals,  trace  of  shell 
fragments,  trace  of  mica,  slightly  calcareous. 

50-51.5  Clayey  Sand:  Same  as  previous  section 

60-61.5  Clayey  Sand:  Same  as  45-45.6  section 


Core  Descriptions 

Middle  Miocene 

Upper  Hawthorne  Formation 

66.2-89.15  Dolomite  to  dolomitic  sand,  very  It.  to  It.  gray  (N7),  massive,  high 

sand  content,  vugs,  moderate  porosity,  hard  dense  dolomite 
at  top  of  section  grading  into  softer,  pitted,  fossiliferous 
dolomitic  limestone  at  80  feet.  Contact  based  on  top  of  rock. 

66.2-66.7  Dolomite:  med.  gray  (N5),  7%  f.g.  sand,  hard,  dense,  small 

vugs  present,  phosphate  present,  low  effective  porosity. 

66.7-67.0  Dolomitic  Sand:  very  It.  gray  (N8),  poorly  indurated, 

argillaceous,  30%  calcareous  clay,  5%  shell  fragments, 
0.5%  heavy  minerals,  phosphate  present. 

67.0-67.8  Dolomite:  very  It.  gray  (N8),  arenaceous,  50%  f.g.  sand, 

subrounded,  medium  hard,  low  effective  porosity. 

67.8-68.92  Dolomite:  very  It.  gray  (N8),  arenaceous,  50%  v.f.g.  to  f.g. 

sand,  subrounded,  0.5%  heavy  minerals,  several  small 
vugs,  low  effective  porosity. 

68.92-69.72  Dololmitic  Sand:  yellowish  gray  (5Y8/1),  arenaceous,  50% 

f.g.  to  m.g.  sand,  0.5%  heavy  minerals,  phosphate  present, 
medium  hard,  moderate  effective  porosity. 

69.72-70.6  Dolomite:  It.  gray  (N7),  arenaceous,  45%  f.g.  sand, 

subangular,  7%  calcareous  clay,  0.5%  heavy  minerals,  2% 
fossils,  several  small  vugs,  low  effective  porosity. 

70.6-  75.5  Cavity 

75.5-76.1  Dolomite:  1l  gray  (N7),  hard,  50%  f.g.  sand,  subrounded 

grains,  several  vugs  ( 1  inch  diameter)  not  interconnected, 
poor  effective  porosity. 

76.1-76.9  Dolomite:  It.  gray  (N7),  hard,  50%  f.g.  sand,  subrounded 

grains,  0.5%  heavy  minerals,  several  horizontal  fractures, 
small  vugs,  moderate  effective  porosity. 
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76.9-77.9  Dolomite:  It.  gray  (N7),  50-60%  f.g.  sand,  0.5%  heavy 

minerals,  subrounded  grains,  several  fractures  and  joints, 
small  vugs,  med.  hard,  moderate  effective  porosity. 

77.9-79.6  Dolomite:  It.  gray  (N7),  hard,  45%  f.g.  sand,  1%  heavy 

minerals,  10%  calcareous  clay,  fractured,  jointed,  vuggy, 
good  secondary  porosity. 

79.6-80.1  Core  loss 

80.1-80.7  Limestone:  white  (N9),  soft,  50%  f.g.  sand,  poorly 

indurated,  10%  calcareous  clay  lenses,  semi-plastic, 
fractured,  fair  effective  porosity. 

80.7-81.2  Dolomitic  Limestone:  It.  gray  (N7),  50%  f.g.-m.g.  sand, 

30%  weathered  shells  &  molds,  fossiliferous,  1%  mica, 
0.5%  heavy  minerals,  highly  weathered  and  fractured,  high 
effective  porosity. 

81.2-81.9  Dolomite:  It.  gray  (N7),  hard,  40%  f.g.  sand,  20% 

weathered  shell  fragments  (pelecepods),  0.5%  heavy 
minerals,  phosphate  present,  many  small  vugs  (.25mm), 
weathered,  highly  pitted,  good  effective  porosity. 

81.9-  82.4  Dolomite:  It.  gray  (N7),  medium  hard,  50-60%  f.g.  sand, 

1%  heavy  minerals,  highly  pitted,  many  small  vugs. 

82.4-83.0  Dolomitic  Sand:  It.  gray  (N7),  50%  f.g.  sand,  subrounded 

to  subangular,  poorly  indurated,  fractured,  weathered,  high 
effective  porosity. 

83.0-83.4  Dolomite:  same  as  8 1 .9-82.4  section 

83.4-86.2  Dolomite:  It.  gray  (N7),  hard,  50-60%  f.g.-m.g  sand, 

subrounded  grains,  1%  heavy  minerals,  1%  very  coarse 
smokey  qtz.  grains,  horizontal  fractures,  some  vugs, 
moderate  effective  porosity. 

86.2-88.0  Dolomitic  Sand:  It.  gray  (N7),  soft,  poorly  indurated,  f.g.- 

m.g  subrounded  grains,  1%  heavy  minerals,  1%  mica 
flakes,  fine  sand  seams,  fractured,  weathered,  high  effective 
porosity. 

88.0-88.4  Core  loss 

88.4-89.15  Same  as  86.2-88.0  section. 


89.15-94.0       Dolomite:  lt.gray  (N7),  harder,  tighter  section,  sandy,  very  pitted,  vuggy, 
very  fossiliferous,  weathered,  good  effective  porosity.  Contact  based  on 
denser  rock  and  rapid  increase  in  fossil  content. 
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89.15-92.3  Dolomite:  It.  gray  (N7),  35%  f.g.-m.g.  sand,  subangular  to 

subrounded,  10%  very  coarse  smokey  qtz.  grains, 
subrounded  sand,  20%  weathered  shell  frags. 
(pelecepods),l%  heavy  minerals,  0.5%  mica  flakes, 
fossiliferous,  vuggy,  fractured,  pitted,  good  effective 
porosity. 

92.3-92.6  Dolomite:  It.  gray  (N7),  40%  f.g.-m.g.  sand,  subrounded  to 

subangular,  10%  v.c.  smokey  qtz.  grains,  poorly  indurated, 
surgary,  1%  heavy  minerals,  fair  effective  porosity. 

92.6-94.0  Dolomite:  It.  gray  (N7),  50%  f.-m.g.  sand,  surgary,  20  % 

fossils,  shells,  1%  heavy  minerals,  med.  hard  texture, 
pitted,  vuggy,  fossiliferous,  fractured,  fair  effective 
porosity. 

94.0-102.3      Dolomite:  It.  gray  (N7),  argillaceous,  arenaceous,  highly  bioturbated 

(Ophiomorphia  nodosa),  high  phosphate  content.  Gradational  contact  based 
on  increase  in  clay  content,  bioturbation  and  decrease  in  shell  content. 

94.0-97.3  Dolomite:  It.  gray  (N7-8),  argillaceous,  arenaceous,  35% 

v.f-f.g.  sand,  subrounded  to  subangular,  15%  clay,  5% 
coarse  qtz.  grains  and  phosphatic  calcite,  1  %  heavy 
minerals,  1%  mica  flakes,  phosphate  present,  mottled 
(Callianassa  burrows)  continuous  core,  low  permeability, 
good  confining  unit. 

97.3-98.3  Dolomite:  It.  gray  (N7-8),  argillaceous,  arenaceous,  40% 

v.f.-f.g.  sand,  subrounded  to  subangular,  40%  clay,  5% 
e.g.  qtz.  and  phosphatic  calcite,  2%  heavy  minerals,  plastic, 
continuous,  good  confining  unit. 

98.3-102.3  Dolomite:  It.  gray  (N7-8),  argillaceous,  arenaceous,  50% 

v.f-f.g.  sand,  subrounded  to  subangular,  1%  phosphate, 
many  Callianassa  burrows,  med.  hard,  low  permeability. 
This  section  contains  med.  gray  silty  calcareous  clay  lenses 
(30%),  subangular  silt  grains,  20%  v.c.  qtz  &  phosphatic 
calcite  grains,  low  permeability. 

102.3-102.8  Core  loss 

102.8-105.8     Calcareous  Clay:  very  It.  gray  (N8),  little  to  no  burrowing,  high  clay  and 
sand  content,  soft,  plastic,  continuous.  Contact  based  on  change  in 
lithology  from  dolomite  to  a  sandy  calcareous  clay. 

102.8-105  Calcareous  Clay:  very  It.  gray  (N8),  20%  v.f.-f.g.  sand, 

subangular  to  subrounded,  1%  e.g.  qtz.  &  phosphatic 
calcite,  plastic,  soft,  sparse  shell  fragments,  sulfur  smell, 
low  permeability. 
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105-105.8  Calcareous  Clay:  very  It  gray  (N8),  argillaceous, 

arenaceous,  30%  f.g.  sand,  subangular  to  subrounded,  e.g. 
sand  (5%),  subrounded,  4%  phosphate,  plastic,  soft, 
continuous,  low  permeability. 


105.8-1 13.5     Dolomitic  L.S.:  It.  gray  (N7),  argillaceous,  arenaceous,  soft,  mottled, 
fossiliferous.  Contact  based  on  lithology  change,  and  an  increase  in 
burrowing  and  phosphate  content. 

105.8-109.2  Dolomitic  L.S.:  very  It.  gray  (N8)  with  It.  gray  (N7) 

burrows,  argillaceous,  arenaceous,  30%  f.-m.g.sand, 
subangular  to  subrounded,  20%  clay,  10%  large  qtz.  and 
phosphatized  calcite  grains  (v.c.  to  6  mm  in  diam.),  5% 
phosphate,  25%  shell  molds,  sparse  shell  frags.,  mottled 
(Callianassa  burrows),  poorly  indurated,  plastic, 
fossiliferous,  med.  soft,  low  permeability. 

109.2-112.2  Dolomitic  L.S.:  It.  gray  (N7),  argillaceous,  arenaceous,  40% 

f-m.g.  sand,  subrounded  to  subangular,  20%  clay,  15%  v.c. 
•     qtz  &  phosphatized  calcite  grains,  20%  shell  molds  &  casts, 
3%  phosphate,  fossiliferous,  friable,  low  permeability. 

1 12.2-1 13.5  Dolomitic  L.S.:  It.  gray  (N7),  alt.  sand  &  clay  lenses,  50% 

f.g.  sand,  subangular  to  subrounded,  25%  clay,  10%  v.c. 
qtz.  and  phosphatized  calcite  grains,  3%  phosphate,  .5% 
mica,  sparse  fossils,  burrows  present,  friable,  low 
permeability. 

113.5-120.3     Dolomitic  L.S.:    It.  gray  (N7),  arenaceous,  argillaceous,  highly  bioturbated 
from  1 17.2-1 19.5,  high  phosphate  and  coarse  quartz  grain  content. 
Gradational  contact  based  on  increase  bioturbation  and  phosphate  content. 

1 13.5-1 17.2  Dolomitic  L.S.:  It  gray  (N7)  with  med.  It.  gray  (N6) 

burrows,  40%  v.f.-f.g.  sand,  20%  clay,  20%  v.c.  qtz.  and 
phosphatized  calcite,  3%  particulate  phosphate,  0.5%  mica, 
abundant  shell  molds,  fossiliferous  (mostly  pelecepods), 
some  Callianassa  burrows,  friable,  moderate  hard,  low 
permeability. 

1 17.2-120.3  Dolomitic  L.S.:  It.  gray  (N7)  with  med.  It.  gray  (N6) 

burrows,  35%  v.f.-f.g.  sand,  subangular  to  subrounded, 
20%  clay,  20%  v.c.  qtz.  and  phos.  calcite,  subrounded, 
10%  weathered  fossil  molds  &  casts,  highly  bioturbated, 
poorly  indurated,  plastic,  med.  hard,  continuous,  low 
permeability.  20%  of  this  section  contains  med.  gray  silty 
clay  lenses  w/  20%  v.c.  qtz.  grains. 
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120.3-127.6    Dolomitic  L.S.:    It.  gray  (N7),  arenaceous,  argillaceous,  fossiliferous,  low 
permeability.  Gradational  contact  based  on  decrease  bioturbation, 
phosphate,  and  coarse  quartz  grains. 


120.3-121.5 


121.5-122.2 


Dolomitic  L.S.:  It.  gray  (N7),  35%  v.f.-f.g.  sand, 
subrounded  to  subangular,  20%  clay,  7%  v.c.  qtz  &  phos. 
calcite  grains,  burrows  present,  med.  hard,  continuous,  low 
permeability. 

Core  Loss 


122.2-126.4 


126.4-126.6 


126.6-128.3 


Dolomitic  L.S.:  It.  gray  (N7),  35%  f.g.  sand,  20%  clay, 
30%  shell  molds,  5%  v.c.  qtz.  &  phos.  calcite  grains,  .5% 
mica  flakes,  med.  hard,  fossiliferous  (Chesapecten  cf.,  and 
Mercenaria  prodroma  cf.),  friable,  tiny  vugs  present,  low 
permeability. 

Dolomitic  L.S.:  It.  gray  (N7),  argillaceous,  arenaceous,  40% 
v.f.-f.g.  sand,  15%  clay,  7%  v.c.  qtz.  &  phos.  calcite 
grains,  50%  shell  molds  (clams  &  pelecepods),  semiplastic, 
low  permeability. 

Dolomitic  L.S.:  It.  gray  (N7),  argillaceous,  arenaceous,  40% 
f.-m.g.  sand,  subrounded  to  subangular,  30%  clay,  10% 
v.c.  qtz.  &  phos.  calcite  grains,  20%  shell  molds  (pin  shells, 
pelecepods),  5%  phosphatic  fossiliferous  (brachiopod 
Lingula),  continuous,  bioturbated,  low  permeability. 


128.3-135.9     Dolomitic  L.S.:  It.  gray  (N7)  argillaceous,  arenaceous,  highly  bioturbated, 
increase  in  phosphate,  low  permeability.  Gradational  contact  based  on 
increase  bioturbation,  phosphate. 


128.3-129.8 


129.8-130.9 


130.9-131.6 


135.9-136.7 


Dolomitic  L.S.:  It.  gray  (N7),  argillaceous,  arenaceous, 
45%  f.-m.g.  sand,  subangular  to  subrounded,  45%  clay, 
7%  v.c.  qtz.  &  phos.  calcite,  2%  shell  molds,  bioturbated, 
low  permeability. 

Dolomitic  L.S.:  It.  gray  (N7),  argillaceous,  arenaceous,  40% 
f.-m.g.  sand,  subrounded  to  subangular,  20%  clay,  15% 
v.c.  qtz.  &  phos.  calcite  grains,  15%  shell  molds, 
bioturbated,  fossiliferous,  low  permeability. 

Dolomitic  L.S.:  It.  gray  (N7),  argillaceous,  arenaceous,  45% 
f.-m.g.  sand,  subrounded  to  subangular,  20%  clay,  15% 
v.c.  qtz.  &  phos.  calcite  grains,  10%  shell  molds  &  casts, 
bioturbated,  poorly  indurated,  low  permeability. 

Core  loss 
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136.7-145.7  Calcareous  Sand:  It.  gray  (N7),  poorly  indurated  loose  quartz  sand,  poorly 
sorted,  good  permeability,  fair  effective  porosity.  Contact  based  on  change 
in  lithology. 

136.7-  137.6  Calcareous  Sand:  It.  gray  (N7),  poorly  sorted  sand,  45%  c- 

v.c.  sand,  35%  m.g.  sand,  subangular  to  subrounded,  2% 
clay,  poorly  indurated,  friable,  good  permeability,  fair 
effective  porosity. 

137.6-139.1  Calcareous  Sand:  alternating  layers  of  poorly  sorted 

calcareous  sand  (It.  gray  (N7))  &  calcareous  fat  clay  (very  It. 
gray  (N8)),  the  clay  lenses  (30%)  are  about  1"  in  thickness, 
35%  e.g.  sand,  subrounded,  25%  f.-m.g.  sand,  subangular 
to  subrounded,  10%  v.c.  sand,  subrounded,  friable,  poorly 
indurated,  fair  permeability,  fair  effective  porosity. 

139.1  -139.7  Calcareous  Sand:  It.  gray  (N7),  60%  c.-v.c.  sand, 

subrounded,  30%  f.-m.g.  sand,  subangular  to  subrounded, 
2%  clay,  clean,  friable,  poorly  sorted,  poorly  indurated, 
good  permeability,  fair  effective  porosity. 

139.7-140  Dolomitic  L.S.:  very  It.  gray  (N8),  arenaceous,  argillaceous, 

40%  f.-m.g.  sand,  subrounded,  30%  v.c.g.  sand, 
subrounded,  poorly  sorted  sand,  15%  clay,  med.  hard,  vugs 
present,  low  permeability. 

140-141.1  Calcareous  sand:  It.  gray  (N7),  poorly  sorted  sand,  65% 

e.g.  sand,  subrounded,  20%  f.-m.g.  sand,  subangular  to 
subrounded,  10%  v.c.g.  sand,  subrounded,  2%  clay,  0.5% 
mica,  poorly  indurated  to  loose,  friable,  good  permeability, 
fair  effective  porosity. 

141.1-142  Core  Loss 

142-143.8  Calcareous  Sand:  It.  gray  (N7),  poorly  sorted  sand,  60% 

e.g.  sand,  subrounded,  15%  f.-m.g.  sand,  subangular  to 
subrounded,  5%  v.c.g.  sand,  subrounded,  2%  shell  molds 
&  casts,  5%  thin  tan  gray  fat  clay  lenses  about  1/2  inch 
thick,  friable,  poorly  indurated  to  loose,  good  permeability, 
fair  effective  porosity. 

143.8-144.1  Dolomitic  L.S.:  very  It.  gray  (N8),  alternating  layers  of 

calcareous  poorly  sorted  sand  and  calcareous  tan  gray  clay, 
60%  e.g.  subrounded  sand,  15%  f.-m.g.  subangular  to 
subrounded  sand,  5%  e.g.  subrounded  sand,  10% 
calcareous  clay,  med.  hard,  friable,  bioturbated,  fair 
permeability,  fair  effective  porosity. 

144.1-145.7  Calcareous  Sand:  It.  gray  (N7),  poorly  sorted  sand,  60% 

e.g.  rounded  sand,  15%  m.g.  subangular  sand,  5%  v.c.g. 
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rounded  sand,  5%  calcareous  tan  gray  clay  (1/4"  to  1.5"  in 
thickness),  poorly  indurated  to  loose,  good  permeability,  fair 
effective  porosity. 


145.7-155.5     Dolomitic  L.S.:  It.  gray  (N7),  arenaceous,  argillaceous,  bioturbated,  soft, 
good  confining  unit.  Contact  based  on  lithology  change. 


145.7-148.2 


148.2-148.6 


148.6-150.3 


150.3-153 


153-155.5 


Dolomitic  L.S.:  It.  gray  (N7),  40%  e.g.  subrounded  sand, 
15%  v.f.-m.g.  subrounded,  5%  v.c.  subrounded  sand,  15% 
clay,  3%  shell  molds,  continuous,  friable,  some  vugs,  med. 
hard,  low  permeability 

Dolomitic  L.S.:  It.  gray  (N7),  40%  e.g.  subrounded  sand, 
15%  v.f.-m.g.  subrounded  ,  5%  v.c.  subrounded  sand, 
15%  clay,  3%  shell  molds,  poorly  indurated  to  loose,  low 
permeability. 

Dolomitic  L.S.:  It.  gray  (N7),  45%  f.-m.g.  sand, 
subrounded  to  subangular,  10%  v.c.g.  subrounded  sand, 
25%  clay,  0.5%  mica  flakes,  0.5%  large  qtz.  grains  (8mm  in 
diam.),  2%  brachiopods,  20%  shell  molds,  fossiliferous 
(Chesapecten  cf.,  and  brachiopod  Lingula),  bioturbated, 
continuous,  low  permeability. 

Dolomitic  L.S.:  It.  gray  (N7),  50%  f.-m.g.  sand, 
subrounded  to  subangular,  10%  v.c.g.  subrounded  sand, 
25%  clay,  0.5%  mica  flakes,  0.5%  large  qtz.  grains  (8mm  in 
diam.),  2%  phosphate,  2%  brachiopods,  3%  shell  molds, 
mottled,  continuous,  low  permeability. 

Dolomitic  L.S.:  It.  gray  (N7),  40%  f.-m.g.  sand, 
subrounded  to  subangular,  15%  v.c.g.  sand,  subrounded, 
25%  clay,  2%  phosphate,  1%  mica  flakes,  med.soft, 
continuous,  mottled,  low  permeability,  good  confining  unit. 
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Description  of  split  spoon  samples  and  core  samples  from  well  KB  MP  1.  Site  1. 

Driller:  Corps  of  Engineers,  Savannah,  GA 
Latitude/Longitude:  30o43'ir,  8r28'13" 
Altitude:  4.59  ft.  above  mean  low  water 
Total  Depth:  155.5  ft. 

(Colors  were  wet  and  matched  to  the  Rock  Color  Chart,  GSA,  1984) 

Slit  spoon  samples 


Depth  in 
ft.  below 
land  surface 

Holocene  to  Pleistocene 
Undifferentiated 


Description 


5.0-6.5 

10-11.5 
15-16.5 

20-21.5 

23.0-23.5 

25-26.5 

30-31.5 
35-36.5 


Sand:  Yellowish  gray  (5Y7/2),  v.f.g.  sand,  subangular, 
well  sorted,  20%  heavy  minerals,  20%  organics. 

Sand:  Same  as  previous  section 

Sand:  Yellowish  gray  (5Y7/2),  v.f.g.  sand,  subangular, 
well  sorted,  20%  heavy  minerals,  20%  organics,  trace  of 
shell  fragments. 

Sand:  Olive  gray  (5Y4/1),  v.f.g.,  subangular,  20%  heavy 
minerals,  10  clay,  1%  mica  flakes. 

Clay:  Dk.  gray  (N3)  fat,  clean,  trace  of  shell  fragments, 
trace  of  mica 

Sandy  Clay:  Med.  gray  (N5),  v.f.g.  sand,  subangular,  alt. 
layers  of  fat  gray  clay  2-3  inch,  in  thickness  (30%),  slightiy 
calcareous. 

Sandy  Clay:  Med.  gray  (N5),  v.f.g.  sand,  subangular,  alt. 
layers  of  fat  gray  clay  2-3  inch,  in  thickness  (30%),  10% 
organic  material,  slightly  calcareous,  trace  of  shell 
fragments. 

Same  as  30-31.5  section 


Pliocene  Series 
40-41.5 


Clayey  Sand:  Brown  gray  (5YR4/1),  v.f.-f.g.,  subangular, 
20%  shell  fragments,  10%  clay,  slightly  calcareous. 
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45-46.5  Clayey  Sand:  Mod.  olive  brown  (5Y4/4),  f.g.-m.g.  sand, 

subangular,  15%  clay,  10%  heavy  minerals,  trace  of  shell 
fragments,  trace  of  mica,  slightly  calcareous. 

50-5 1.5  Clayey  Sand:  Same  as  previous  section 

60-61.5  Clayey  Sand:  Same  as  45-45.6  section 


Core  Descriptions 

Middle  Miocene 

Upper  Hawthorne  Formation 

66.2-89.15  Dolomite  to  dolomitic  sand,  very  It.  to  It.  gray  (N7),  massive,  high 

sand  content,  vugs,  moderate  porosity,  hard  dense  dolomite 
at  top  of  section  grading  into  softer,  pitted,  fossiliferous 
dolomitic  limestone  at  80  feet.  Contact  based  on  top  of  rock. 

66.2-66.7  Dolomite:  med.  gray  (N5),  7%  f.g.  sand,  hard,  dense,  small 

vugs  present,  phosphate  present,  low  effective  porosity. 

66.7-67.0  Dolomitic  Sand:  very  It.  gray  (N8),  poorly  indurated, 

argillaceous,  30%  calcareous  clay,  5%  shell  fragments, 
0.5%  heavy  minerals,  phosphate  present. 

67.0-67.8  Dolomite:  very  It.  gray  (N8),  arenaceous,  50%  f.g.  sand, 

subrounded,  medium  hard,  low  effective  porosity. 

67.8-68.92  Dolomite:  very  It.  gray  (N8),  arenaceous,  50%  v.f.g.  to  f.g. 

sand,  subrounded,  0.5%  heavy  minerals,  several  small 
vugs,  low  effective  porosity. 

68.92-69.72  Dololmitic  Sand:  yellowish  gray  (5Y8/1),  arenaceous,  50% 

f.g.  to  m.g.  sand,  0.5%  heavy  minerals,  phosphate  present, 
medium  hard,  moderate  effective  porosity. 

69.72-70.6  Dolomite:  It.  gray  (N7),  arenaceous,  45%  f.g.  sand, 

subangular,  7%  calcareous  clay,  0.5%  heavy  minerals,  2% 
fossils,  several  small  vugs,  low  effective  porosity. 

70.6-  75.5  Cavity 

75.5-76.1  Dolomite:  It.  gray  (N7),  hard,  50%  f.g.  sand,  subrounded 

grains,  several  vugs  ( 1  inch  diameter)  not  interconnected, 
poor  effective  porosity. 

76.1-76.9  Dolomite:  It.  gray  (N7),  hard,  50%  f.g.  sand,  subrounded 

grains,  0.5%  heavy  minerals,  several  horizontal  fractures, 
small  vugs,  moderate  effective  porosity. 
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76.9-77.9  Dolomite:  It.  gray  (N7),  50-60%  f.g.  sand,  0.5%  heavy 

minerals,  subrounded  grains,  several  fractures  and  joints, 
small  vugs,  med.  hard,  moderate  effective  porosity. 

77.9-79.6  Dolomite:  It.  gray  (N7),  hard,  45%  f.g.  sand,  1%  heavy 

minerals,  10%  calcareous  clay,  fractured,  jointed,  vuggy, 
good  secondary  porosity. 

79.6-80.1  Core  loss 

80.1-80.7  Limestone:  white  (N9),  soft,  50%  f.g.  sand,  poorly 

indurated,  10%  calcareous  clay  lenses,  semi-plastic, 
fractured,  fair  effective  porosity. 

80.7-81.2  Dolomitic  Limestone:  It.  gray  (N7),  50%  f.g.-m.g.  sand, 

30%  weathered  shells  &  molds,  fossiliferous,  1%  mica, 
0.5%  heavy  minerals,  highly  weathered  and  fractured,  high 
effective  porosity. 

81.2-81.9  Dolomite:  It.  gray  (N7),  hard,  40%  f.g.  sand,  20% 

weathered  shell  fragments  (pelecepods),  0.5%  heavy 
minerals,  phosphate  present,  many  small  vugs  (.25mm), 
weathered,  highly  pitted,  good  effective  porosity. 

81.9-  82.4  Dolomite:  It.  gray  (N7),  medium  hard,  50-60%  f.g.  sand, 

1%  heavy  minerals,  highly  pitted,  many  small  vugs. 

82.4-83.0  Dolomitic  Sand:  It.  gray  (N7),  50%  f.g.  sand,  subrounded 

to  subangular,  poorly  indurated,  fractured,  weathered,  high 
effective  porosity. 

83.0-83.4  Dolomite:  same  as  81.9-82.4  section 

83.4-86.2  Dolomite:  It.  gray  (N7),  hard,  50-60%  f.g.-m.g  sand, 

subrounded  grains,  1%  heavy  minerals,  1%  very  coarse 
smokey  qtz.  grains,  horizontal  fractures,  some  vugs, 
moderate  effective  porosity. 

86.2-88.0  Dolomitic  Sand:  It.  gray  (N7),  soft,  poorly  indurated,  f.g.- 

m.g  subrounded  grains,  1  %  heavy  minerals,  1  %  mica 
flakes,  fine  sand  seams,  fractured,  weathered,  high  effective 
porosity. 

88.0-88.4  Core  loss 

88.4-89.15  Same  as  86.2-88.0  section. 


89.15-94.0       Dolomite:  lt.gray  (N7),  harder,  tighter  section,  sandy,  very  pitted,  vuggy, 
very  fossiliferous,  weathered,  good  effective  porosity.  Contact  based  on 
denser  rock  and  rapid  increase  in  fossil  content. 
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89.15-92.3  Dolomite:  It.  gray  (N7),  35%  f.g.-m.g.  sand,  subangular  to 

subrounded,  10%  very  coarse  smokey  qtz.  grains, 
subrounded  sand,  20%  weathered  shell  frags. 
(pelecepods),l%  heavy  minerals,  0.5%  mica  flakes, 
fossiliferous,  vuggy,  fractured,  pitted,  good  effective 
porosity. 

92.3-92.6  Dolomite:  It.  gray  (N7),  40%  f.g.-m.g.  sand,  subrounded  to 

subangular,  10%  v.c.  smokey  qtz.  grains,  poorly  indurated, 
surgary,  1%  heavy  minerals,  fair  effective  porosity. 

92.6-94.0  Dolomite:  It.  gray  (N7),  50%  f.-m.g.  sand,  surgary,  20  % 

fossils,  shells,  1%  heavy  minerals,  med.  hard  texture, 
pitted,  vuggy,  fossiliferous,  fractured,  fair  effective 
porosity. 


94.0-102.3      Dolomite:  It.  gray  (N7),  argillaceous,  arenaceous,  highly  bioturbated 

(Ophiomorphia  nodosa),  high  phosphate  content.  Gradational  contact  based 
on  increase  in  clav  content,  hiomrhafion  and  decrease  in  shell  content. 


on  increase  in  clay  content,  bioturbation  and  decrease  in  shell  content. 


94.0-97.3  Dolomite:  It.  gray  (N7-8),  argillaceous,  arenaceous,  35% 

v.f-f.g.  sand,  subrounded  to  subangular,  15%  clay,  5% 
coarse  qtz.  grains  and  phosphatic  calcite,  1  %  heavy 
minerals,  1%  mica  flakes,  phosphate  present,  mottled 
(Callianassa  burrows)  continuous  core,  low  permeability, 
good  confining  unit. 

97.3-98.3  Dolomite:  It.  gray  (N7-8),  argillaceous,  arenaceous,  40% 

v.f.-f.g.  sand,  subrounded  to  subangular,  40%  clay,  5% 
e.g.  qtz.  and  phosphatic  calcite,  2%  heavy  minerals,  plastic, 
continuous,  good  confining  unit. 

98.3-102.3  Dolomite:  It.  gray  (N7-8),  argillaceous,  arenaceous,  50% 

v.f-f.g.  sand,  subrounded  to  subangular,  1%  phosphate, 
many  Callianassa  burrows,  med.  hard,  low  permeability. 
This  section  contains  med.  gray  silty  calcareous  clay  lenses 
(30%),  subangular  silt  grains,  20%  v.c.  qtz  &  phosphatic 
calcite  grains,  low  permeability. 

102.3-102.8  Core  loss 

102.8-105.8     Calcareous  Clay:  very  It.  gray  (N8),  little  to  no  burrowing,  high  clay  and 
sand  content,  soft,  plastic,  continuous.  Contact  based  on  change  in 
lithology  from  dolomite  to  a  sandy  calcareous  clay. 

102.8-105  Calcareous  Clay:  very  It.  gray  (N8),  20%  v.f.-f.g.  sand, 

subangular  to  subrounded,  1%  e.g.  qtz.  &  phosphatic 
calcite,  plastic,  soft,  sparse  shell  fragments,  sulfur  smell, 
low  permeability. 
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105-105.8  Calcareous  Clay:  very  It  gray  (N8),  argillaceous, 

arenaceous,  30%  f.g.  sand,  subangular  to  subrounded,  e.g. 
sand  (5%),  subrounded,  4%  phosphate,  plastic,  soft, 
continuous,  low  permeability. 


105.8-1 13.5     Dolomitic  L.S.:  It.  gray  (N7),  argillaceous,  arenaceous,  soft,  mottled, 
fossiliferous.  Contact  based  on  lithology  change,  and  an  increase  in 
burrowing  and  phosphate  content. 

105.8-109.2  Dolomitic  L.S.:  very  It.  gray  (N8)  with  It.  gray  (N7) 

burrows,  argillaceous,  arenaceous,  30%  f.-m.g.sand, 
subangular  to  subrounded,  20%  clay,  10%  large  qtz.  and 
phosphatized  calcite  grains  (v.c.  to  6  mm  in  diam),  5% 
phosphate,  25%  shell  molds,  sparse  shell  frags.,  mottled 
(Callianassa  burrows),  poorly  indurated,  plastic, 
fossiliferous,  med.  soft,  low  permeability. 

109.2-1 12.2  Dolomitic  L.S.:  It.  gray  (N7),  argillaceous,  arenaceous,  40% 

f-m.g.  sand,  subrounded  to  subangular,  20%  clay,  15%  v.c. 
qtz  &  phosphatized  calcite  grains,  20%  shell  molds  &  casts, 
3%  phosphate,  fossiliferous,  friable,  low  permeability. 

1 12.2-1 13.5  Dolomitic  L.S.:  It.  gray  (N7),  alt.  sand  &  clay  lenses,  50% 

f.g.  sand,  subangular  to  subrounded,  25%  clay,  10%  v.c. 
qtz.  and  phosphatized  calcite  grains,  3%  phosphate,  .5% 
mica,  sparse  fossils,  burrows  present,  friable,  low 
permeability. 

113.5-120.3     Dolomitic  L.S.:    It.  gray  (N7),  arenaceous,  argillaceous,  highly  bioturbated 
from  1 17.2-1 19.5,  high  phosphate  and  coarse  quartz  grain  content. 
Gradational  contact  based  on  increase  bioturbation  and  phosphate  content 

1 13.5-1 17.2  Dolomitic  L.S.:  It  gray  (N7)  with  med.  It.  gray  (N6) 

burrows,  40%  v.f.-f.g.  sand,  20%  clay,  20%  v.c.  qtz.  and 
phosphatized  calcite,  3%  particulate  phosphate,  0.5%  mica, 
abundant  shell  molds,  fossiliferous  (mostly  pelecepods), 
some  Callianassa  burrows,  friable,  moderate  hard,  low 
permeability. 

1 17.2-120.3  Dolomitic  L.S.:  It  gray  (N7)  with  med.  It  gray  (N6) 

burrows,  35%  v.f.-f.g.  sand,  subangular  to  subrounded, 
20%  clay,  20%  v.c.  qtz.  and  phos.  calcite,  subrounded, 
10%  weathered  fossil  molds  &  casts,  highly  bioturbated, 
poorly  indurated,  plastic,  med.  hard,  continuous,  low 
permeability.  20%  of  this  section  contains  med.  gray  silty 
clay  lenses  w/  20%  v.c.  qtz.  grains. 
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120.3-127.6    Dolomitic  L.S.:    It.  gray  (N7),  arenaceous,  argillaceous,  fossiliferous,  low 
permeability.  Gradational  contact  based  on  decrease  bioturbation, 
phosphate,  and  coarse  quartz  grains. 

120.3-121.5  Dolomitic  L.S.:  It.  gray  (N7),  35%  v.f.-f.g.  sand, 

subrounded  to  subangular,  20%  clay,  7%  v.c.  qtz  &  phos. 
calcite  grains,  burrows  present,  med.  hard,  continuous,  low 
permeability. 

121.5-122.2  Core  Loss 

122.2-126.4  Dolomitic  L.S.:  It.  gray  (N7),  35%  f.g.  sand,  20%  clay, 

30%  shell  molds,  5%  v.c.  qtz.  &  phos.  calcite  grains,  .5% 
mica  flakes,  med.  hard,  fossiliferous  (Chesapecten  cf.,  and 
Mercenaria  prodroma  cf.),  friable,  tiny  vugs  present,  low 
permeability. 

126.4-126.6  Dolomitic  L.S.:  It.  gray  (N7),  argillaceous,  arenaceous,  40% 

v.f.-f.g.  sand,  15%  clay,  7%  v.c.  qtz.  &  phos.  calcite 
grains,  50%  shell  molds  (clams  &  pelecepods),  semiplastic, 
low  permeability. 

126.6-128.3  Dolomitic  L.S.:  It.  gray  (N7),  argillaceous,  arenaceous,  40% 

f.-m.g.  sand,  subrounded  to  subangular,  30%  clay,  10% 
v.c.  qtz.  &  phos.  calcite  grains,  20%  shell  molds  (pin  shells, 
pelecepods),  5%  phosphatic  fossiliferous  (brachiopod 
Lingula),  continuous,  bioturbated,  low  permeability. 

128.3-135.9     Dolomitic  L.S.:  It.  gray  (N7)  argillaceous,  arenaceous,  highly  bioturbated, 
increase  in  phosphate,  low  permeability.  Gradational  contact  based  on 
increase  bioturbation,  phosphate. 

128.3-129.8  Dolomitic  L.S.:  It.  gray  (N7),  argillaceous,  arenaceous, 

45%  f.-m.g.  sand,  subangular  to  subrounded,  45%  clay, 
7%  v.c.  qtz.  &  phos.  calcite,  2%  shell  molds,  bioturbated, 
low  permeability. 

129.8-130.9  Dolomitic  L.S.:  It.  gray  (N7),  argillaceous,  arenaceous,  40% 

f.-m.g.  sand,  subrounded  to  subangular,  20%  clay,  15% 
v.c.  qtz.  &  phos.  calcite  grains,  15%  shell  molds, 
bioturbated,  fossiliferous,  low  permeability. 

130.9-131.6  Dolomitic  L.S.:  It.  gray  (N7),  argillaceous,  arenaceous,  45% 

f.-m.g.  sand,  subrounded  to  subangular,  20%  clay,  15% 
v.c.  qtz.  &  phos.  calcite  grains,  10%  shell  molds  &  casts, 
bioturbated,  poorly  indurated,  low  permeability. 

135.9-136.7  Core  loss 
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136.7-145.7  Calcareous  Sand:  It.  gray  (N7),  poorly  indurated  loose  quartz  sand,  poorly 
sorted,  good  permeability,  fair  effective  porosity.  Contact  based  on  change 
in  lithology. 

136.7-  137.6  Calcareous  Sand:  It.  gray  (N7),  poorly  sorted  sand,  45%  c- 

v.c.  sand,  35%  m.g.  sand,  subangular  to  subrounded,  2% 
clay,  poorly  indurated,  friable,  good  permeability,  fair 
effective  porosity. 

1 37.6- 1 39. 1  Calcareous  Sand:  alternating  layers  of  poorly  sorted 

calcareous  sand  (It.  gray  (N7))  &  calcareous  fat  clay  (very  It. 
gray  (N8)),  the  clay  lenses  (30%)  are  about  1"  in  thickness, 
35%  e.g.  sand,  subrounded,  25%  f.-m.g.  sand,  subangular 
to  subrounded,  10%  v.c.  sand,  subrounded,  friable,  poorly 
indurated,  fair  permeability,  fair  effective  porosity. 

139.1  -139.7  Calcareous  Sand:  It.  gray  (N7),  60%  c.-v.c.  sand, 

subrounded,  30%  f.-m.g.  sand,  subangular  to  subrounded, 
2%  clay,  clean,  friable,  poorly  sorted,  poorly  indurated, 
good  permeability,  fair  effective  porosity. 

139.7-140  Dolomitic  L.S.:  very  It.  gray  (N8),  arenaceous,  argillaceous, 

40%  f.-m.g.  sand,  subrounded,  30%  v.c.g.  sand, 
subrounded,  poorly  sorted  sand,  15%  clay,  med.  hard,  vugs 
present,  low  permeability. 

140-141.1  Calcareous  sand:  It.  gray  (N7),  poorly  sorted  sand,  65% 

e.g.  sand,  subrounded,  20%  f.-m.g.  sand,  subangular  to 
subrounded,  10%  v.c.g.  sand,  subrounded,  2%  clay,  0.5% 
mica,  poorly  indurated  to  loose,  friable,  good  permeability, 
fair  effective  porosity. 

141.1-142  Core  Loss 

142-143.8  Calcareous  Sand:  It.  gray  (N7),  poorly  sorted  sand,  60% 

e.g.  sand,  subrounded,  15%  f.-m.g.  sand,  subangular  to 
subrounded,  5%  v.c.g.  sand,  subrounded,  2%  shell  molds 
&  casts,  5%  thin  tan  gray  fat  clay  lenses  about  1/2  inch 
thick,  friable,  poorly  indurated  to  loose,  good  permeability, 
fair  effective  porosity. 

143.8-144.1  Dolomitic  L.S.:  very  It.  gray  (N8),  alternating  layers  of 

calcareous  poorly  sorted  sand  and  calcareous  tan  gray  clay, 
60%  e.g.  subrounded  sand,  15%  f.-m.g.  subangular  to 
subrounded  sand,  5%  e.g.  subrounded  sand,  10% 
calcareous  clay,  med.  hard,  friable,  bioturbated,  fair 
permeability,  fair  effective  porosity. 

144.1-145.7  Calcareous  Sand:  It.  gray  (N7),  poorly  sorted  sand,  60% 

e.g.  rounded  sand,  15%  m.g.  subangular  sand,  5%  v.c.g. 
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rounded  sand,  5%  calcareous  tan  gray  clay  (1/4"  to  1.5"  in 
thickness),  poorly  indurated  to  loose,  good  permeability,  fair 
effective  porosity. 


145.7-155.5     Dolomitic  L.S.:  It.  gray  (N7),  arenaceous,  argillaceous,  bioturbated,  soft, 
good  confining  unit.  Contact  based  on  lithology  change. 


145.7-148.2 


148.2-148.6 


148.6-150.3 


150.3-153 


153-155.5 


Dolomitic  L.S.:  It.  gray  (N7),  40%  e.g.  subrounded  sand, 
15%  v.f.-m.g.  subrounded,  5%  v.c.  subrounded  sand,  15% 
clay,  3%  shell  molds,  continuous,  friable,  some  vugs,  med. 
hard,  low  permeability 

Dolomitic  L.S.:  It.  gray  (N7),  40%  e.g.  subrounded  sand, 
15%  v.f.-m.g.  subrounded  ,  5%  v.c.  subrounded  sand, 
15%  clay,  3%  shell  molds,  poorly  indurated  to  loose,  low 
permeability. 

Dolomitic  L.S.:  It.  gray  (N7),  45%  f.-m.g.  sand, 
subrounded  to  subangular,  10%  v.c.g.  subrounded  sand, 
25%  clay,  0.5%  mica  flakes,  0.5%  large  qtz.  grains  (8mm  in 
diam.),  2%  brachiopods,  20%  shell  molds,  fossiliferous 
(Chesapecten  cf.,  and  brachiopod  Lingula),  bioturbated, 
continuous,  low  permeability. 

Dolomitic  L.S.:  It.  gray  (N7),  50%  f.-m.g.  sand, 
subrounded  to  subangular,  10%  v.c.g.  subrounded  sand, 
25%  clay,  0.5%  mica  flakes,  0.5%  large  qtz.  grains  (8mm  in 
diam.),  2%  phosphate,  2%  brachiopods,  3%  shell  molds, 
mottled,  continuous,  low  permeability. 

Dolomitic  L.S.:  It.  gray  (N7),  40%  f.-m.g.  sand, 
subrounded  to  subangular,  15%  v.c.g.  sand,  subrounded, 
25%  clay,  2%  phosphate,  1%  mica  flakes,  med.soft, 
continuous,  mottled,  low  permeability,  good  confining  unit. 
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Description  of  cuttings  from  Well  KBMP  6.  Site  2 

Driller:  Corps  of  Engineers,  Savannah,  GA 
Latitude/Longitude:  30°43*10",  81027,26*' 
Altitude:  6.86  ft.  above  mean  low  water 
Total  Depth:  170  ft 


(Colors  were  dry  samples  matched  to  the  Rock  Color  Chart,  GSA,  1984) 

Description 


Depth  in 
ft  below 
land  surface 


Holocene  to  Pleistocene 
Undifferentiated 

0-10 


10-20 
20-30 

30-40 

45-50 

Pliocene  Series 
50-61.5 


Sand:  Lt.  olive  gray  (5Y6/1),  v.f.-f.g.,  subangular,  well 
sorted,  5%  heavy  minerals,  trace  of  mica 

Sand:  same  as  previous  section 

Sand:  Lt.  olive  gray  (5Y6/1),  v.f.-f.g.,  subangular,  well 
sorted,  7%  shell  fragments,  5%  heavy  minerals,  trace  of 
mica. 

Sand:  Lt.  gray  (N7),  v.f.-f.g.,  subangular,  well  sorted,  3% 
heavy  minerals,  3%  mica. 

Sandy  Clay:  Med.  lt.  gray  (N6),  35%  v.f.-f.g.  sand, 
subangular,  45%  clay,  3%  heavy  minerals. 


Clayey  Sand:  Lt.  gray  (N7),  v.f.-f.g.,  subangular,  10% 
shell  fragments,  10%  clay,  3%  heavy  minerals. 


Middle  Miocene  Series 
Upper  Charlton  Formation 

61.5-170  Dolomitic  Limestone:  olive  gray  (5Y6/1),  argillaceous,  arenaceous, 

phosphatic,  20%  granule  size  sand.  Contact  based  on  top  of  rock. 


61.5-65 

65-70 
70-80 


Dolomitic  L.S.:  Lt.  olive  gray  (5Y6/1),  45%  v.f.-f.g.  sand, 
subangular,  (wash  down),  10%  broken  shells,  5%  clay,  2% 
phosphate. 

Dolomitic  L.S.:  Same  as  previous  section 

Dolomitic  L.S.:  Lt.  olive  gray  (5Y6/1),  45%  f.g.  sand, 
subangular,  15%  clay,  10%  broken  shell  fragments,  5% 
heavy  minerals,  2%  phosphate. 
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80-90  Dolomitic  L.S.:  Yellow  gray  (5Y8/1),  40%  f.-m.g.  sand, 

subangular,  20%  clay,  5%  shell  fragments,  5%  heavy 
minerals,  2%  phosphate,  1%  mica,  argillaceous. 

90-100  Dolomitic  L.S.:  Lt.  olive  gray  (5Y6/1),  40%  f.-m.g.  sand, 

subangular,  25%  clay,  15%  v.c.  sand  to  fine  gravel, 
subrounded,  5%  heavy  minerals,  5%  phosphate,  1%  mica, 
argillaceous. 

100-1 10  Dolomitic  L.S.:  Lt.  olive  gray  (5Y6/1),  40%  f.-m.g.  sand, 

subangular,  25%  clay,  20%  v.c.  sand  to  fine  gravel, 
subrounded,  7%  phosphate,  3%  heavy  minerals,  1%  mica, 
argillaceous. 

1 10-120  Dolomitic  L.S.:  Lt.  olive  gray  (5Y6/1),  40%  f.-m.g.  sand, 

subangular,  30%  calcareous  clay,  25%  v.c.  sand  to  fine 
gravel,  subrounded,  5%  shells,  5%  phosphate,  argillaceous. 

120-130  Dolomitic  L.S.:  Lt.  olive  gray  (5Y6/1),  40%  f.-m.g.  sand, 

subangular,  30%  clay,  20%  c.  sand  to  fine  gravel, 
subrounded,  7%  phosphate,  3%  heavy  minerals,  1%  mica, 
argillaceous. 

130-140  Dolomitic  L.S.:  Lt.  olive  gray  (5Y6/1),  40%  clay,  35%  f.- 

m.g.  sand,  subangular,  poorly  sorted,  20%  c.  sand  to  fine 
gravel,  subrounded,  3%  phosphate,  3%  heavy  minerals,  1% 
mica,  argillaceous. 

140-145  Dolomitic  L.S.:  Lt.  olive  gray  (5Y6/1),  45%  calcareous 

clay,  30%  f.-m.g.  sand,  subangular,  20%  v.c.  sand  to  fine 
gravel,  subrounded,  5%  phosphate,  3%  heavy  minerals, 
argillaceous. 

145-155  Dolomitic  L.S.:  Lt.  olive  gray  (5Y6/1),  30%  f.-m.g.  sand, 

subangular,  45%  calcareous  clay,  20%  v.c.  sand  to  fine 
gravel,  subrounded,  3%  phosphate,  argillaceous. 

155-160  Dolomitic  L.S.:  Lt.  olive  gray  (5Y6/1),  30%  f.-m.g.  sand, 

subangular,  30%  v.c.  sand  to  fine  gravel,  subrounded,  25% 
clay,  3%  phosphate,  argillaceous. 

160-165  Dolomitic  L.S.:  Lt.  olive  gray  (5Y6/1),  30%  f.-m.g.  sand, 

subangular,  45%  calcareous  clay,  20%  v.c.  sand  to  fine 
gravel,  subrounded,  3%  phosphate,  2%  phosphate, 
argillaceous. 

165-170  Dolomitic  L.S.:  Lt.  olive  gray  (5Y6/1),  30%  f.-m.g.  sand, 

subangular,  50%  clay,  15%  v.c.  sand  to  fine  gravel, 
subrounded,  3%  phosphate,  2%  phosphate,  argillaceous. 
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Description  of  cuttings  from  Well  KB  MP  10.  Site  3 

Driller:  Corps  of  Engineers,  Savannah,  GA 
Latitude/Longitude:  30°44'51",  8r28'0" 
Altitude:  15.93  ft.  above  mean  low  water 
Total  Depth:  133  ft 

(Colors  were  dry  samples  matched  to  the  Rock  Color  Chart,  GSA,  1984) 

Depth  in 

ft  below  Description 

land  surface 


Pleistocene  Series 
5-10 

10-20 
20-30 
30-35 

35-40 

40-45 

45-50 
50-55 
55-60 

60-62 

Pliocene  Series 
Duplin  Formation 
62-80 


Sand:  Grayish  orange  (10YR7/4),  v.f.-f.g.,  well  sorted, 
angular  to  subangular,  10%  shell  fragments,  3%  heavy 
minerals,  l%mica. 

Sand:  same  as  previous  section 

Sand:  same  as  0-10  section 

Clayey  Sand:  Yellowish  gray  (5Y8/1),  50%  v.f.g.  silty 
sand,  subangular,  40%  gray  clay  chunks  approximately 
5mm  in  diameter,  3%  shell  fragments. 

Clayey  Sand:  Yellowish  gray  (5Y8/1),  45%  v.f.g.  silty 
sand,  subangular,  45%  gray  clay,  10%  shell  fragments, 
10%  heavy  minerals. 

Clayey  Sand:  Yellowish  gray  (5Y8/1),  50%  gray  clay,  40% 
v.f.g.  silty  sand,  subangular,  10%  shell  fragments,  10% 
heavy  minerals. 

Clay:  Lt.  gray  (N7),  fat,  clean,  10%  oyster  shells  and  other. 

Clay:  same  as  previous  section. 

Sandy  Clay:  Lt  gray  (N7),  50%  fat  gray  clay,  40%  v.f.-f.g. 
sand,  subangular,  15%  broken  shell  fragments,  3%  heavy 
minerals,  2%mica. 

Sandy  Clay:  same  as  previous  section 


Sand:  lt.  gray  (N7),  poorly  sorted,  subangular  to  subrounded, 
grading  into  clayey  sand  and  clay  towards  the  base  of  this  section, 
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good  permeability  from  65-72'.  Contact  based  on  change  in 
lithology. 

62-65  Clayey  Sand:  Lt.  gray  (N7),  55%  f.-m.g.  sand,  subangular, 

20%  clay,  20%  shell  fragments  often  2mm  in  diameter. 

65-72  Sand:  Lt.  gray  (N7),  f.-c.g.  sand,  subangular  to 

subrounded,  poorly  sorted,  25%  of  sand  is  granule  in  size, 
25%  coarse  shell  fragments,  3%  clay,  good  permeability. 

72-75  Clayey  Sand:  Med.  lt.  gray  (N6),  f.-c.g.  sand,  subangular 

to  subrounded,  poorly  sorted,  25%  coarse  shell  fragments, 
20%  clay,  slightly  calcareous. 

75-80  Sandy  Clay:  Yellowish  gray  (5Y8/1),  40%  clay,  30%  f.-c.g. 

sand,  subangular,  poorly  sorted,  15%  limestone  fragments, 
15%  shell  fragments. 

Middle  Miocene  Series 

Upper  Charlton  Formation 

80- 1 33  Limestone,  fractured  at  the  top  of  rock  becoming  a  very  dense 

dolomite  at  89'.  From  100-132'  sandy  clay  with  dolomite 
fragments,  phosphatic.  A  very  dense  dolomite  is  present  at  the  base 
of  this  section.  Contact  based  on  top  of  rock. 

80-88  Limestone:  Very  lt.  gray  (N8),  15%  f.-c.g.  sand, 

subangular,  15%  broken  shell  fragments,  3%  clay, 
fossiliferous. 

89-95  Dolomite:  Med.  lt.  gray  (N6),  very  hard,  slow  drilling  time. 

95-100  Dolomite:  Lt.  gray  (N7),  30%  clay,  20%  oyster,  pelecepod 

&  other  shell  fragments,  10%  f.-c.g.  sand,  subangular  to 
subrounded,  3%  phosphate,  argillaceous. 

100-1 10  Sandy  Clay:  Lt.  gray  (N7),  40%  gray  clay,  30%  v.f.-c. 

sand,  subangular  to  subrounded,  20%  dolomite  fragments, 
10%  shell  fragments,  5%  phosphate  (right  on  with  gamma), 
slighdy  calcareous. 

1 10-1 15  Sandy  Clay:  Lt.  gray  (N7),  40%  gray  clay,  25%  m.-c.  sand, 

subangular  to  subrounded,  poorly  sorted,  3%  heavy 
minerals,  2%  phosphate,  slightly  calcareous. 

1 15- 120  Sandy  Clay:  Lt.  gray  (N7),  40%  gray  clay,  40%  v.f.-c. 

sand,  subangular  to  subrounded,  poorly  sorted,  3%  shell 
fragments,  3%  heavy  minerals,  1%  phosphate,  slightly 
calcareous. 
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120-132  Clayey  Sand:  Lt.  gray  (N7),  40%  gray,  fat,  clay,  35%  v.f.- 

c.  sand,  subangular  to  subrounded,  poorly  sorted,  10%  fine 
gravel,  subrounded,  7%  oyster  &  other  shell  fragments,  5% 
phosphate,  1%  heavy  minerals,  slighdy  calcareous. 

133  Dolomite:  Ll  gray  (N7),  30%  clay,  20%  f.-c.g.  sand, 

subangular,  15%  gravel,  subrounded,  5%  shell  fragments, 
3%  phosphate,  argillaceous,  hard,  slow  drill  time. 
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APPENDIX  C 


Estimation  of  Hydraulic  Conductivity  Based  on  Grain  Size  Analysis 
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Masch  and  Denny  (1966)  presented  a  method  in  which  the  standard  deviation  of  the 
sample  is  incorporated  in  the  process  of  determining  the  hydraulic  conductivity.  Estimation 
of  hydraulic  conductivity  is  more  accurate  when  the  standard  deviation  of  the  gradation 
curve  is  taken  into  account.  First,  the  gradation  curve  is  plotted  using  Krumbein's  0  units 
(see  figure  CI).  Next,  the  standard  deviation  of  the  sample  is  calculated  by: 


s  =  _dl6^L_d84  +  A5^A95  (CI) 

4  6.6 


Knowing  the  median  grain  size  (dso)  derived  from  the  grain  size  analysis  (figure  CI)  and 
the  standard  deviation  of  the  sample,  the  hydraulic  conductivity  may  be  estimated  from 
graph  C2.  The  curves  in  figure  C2  were  developed  experimentally  by  testing  a  group  of 
unconsolidated  sand  samples  (Masch  and  Denny,  1966). 

The  method  employed  to  approximate  hydraulic  conductivity  and  does  not  take  into 
account  factors  such  as  natural  particle  arrangement,  stratification,  or  the  anisotropy  of  the 
zone.  Also,  the  samples  tested  are  not  representative  of  the  grain-size  characteristics  of  the 
entire  zone.  Although  these  factors  are  not  considered,  the  approximation  for  hydraulic 
conductivity  based  on  grain  size  provides  a  suitable  range  of  values  for  field  samples. 


135 


m 

-C 

CO 

&, 

H 

C 

Z 

c« 

s 

<U 

►H 

C 

<N   S 

o 

o< 

z 

c 
o 

'.s 

W 

3 

J 

•c 

""*   CO 

oc 

Z 

•o 

5 

04 

CO 

O 

C 

o 

In 

o 

_• 

U 

<u 

i~ 

I 

S 

E 

oooooooooo 


QSNIVIH^  INHDHHd 


136 


8000 


5000  — 


to 

z 
o 


< 


< 


1000  — 


CD 

< 
UJ 

2 

c 

UJ 

o. 


UJ 

o 


UJ 

o 
o 


300  — 


-1.0 


0.0  10  2.0  3.0 

MD50     DIAMETER   «J>) 


4.0 


C2.  Estimation  of  hydraulic  conductivity  once  dso  and  the  standard  deviation  are 
known  (Masch  and  Denny,  1966). 
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Calculation  of  the  standard  deviation  of  each  sample. 

Site  3,  Pliocene  sand  of  the  Pliocene-Miocene  aquifer,  KBMP  9 

1.15 -(-0.9)     +     2.3  -    (-1.5)     =     1.08 
4  6.6 


Site  1,  Miocene  sand,  KBMP  1 

1.8  -(-0.1)     +     2.85  -    (-0.78)     =     1.025 
4  6.6 


Site  2,  Surficial  aquifer,  KBMP  5 

3.21 -(1.45)     +     3.77  -    (0.7)     =     0.905 
4  6.6 


Site  1,  Surficial  Aquifer,  KBMP  2 

2.99  -  (2.07)     +     3.65  -    (1.65)     =     0.533 
4  6.6 


Site  3,  Surficial  Aquifer,  KBMP  8 

3.12 -(2.15)     +     3.7  -    (2.0)     =     0.5 
4  6.6 
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APPENDIX  D 


Density  Corrections 


139 


Density  Corrections 


Well  number 

Percent  seawater 

p;  fgm/cc) 

&  (lb/ft 

KBMP01 

32 

1.008 

62.93 

KBMP02 

0.2 

1.000 

62.43 

KBMP03 

68 

1.017 

63.49 

KBMP04 

73 

1.018 

63.55 

KBMP05 

3.0 

1.001 

62.49 

KBMP06 

89 

1.022 

63.80 

KBMP07 

4.2 

1.001 

62.49 

KBMP08 

0.7 

1.000 

62.43 

KBMP09 

0.2 

1.000 

62.43 

KBMP10 

0.2 

1.000 

62.43 

KBMP11 

5.2 

1.017 

62.49 

The  following  equation  was  used  to  convert  percent  seawater  to  density: 

p  =  (percent  seawater)(0.025  gm/cc)  +  l.Ogm/cc 

The  freshwater  head  in  a  well  containing  saline  water  is  calculated  by  the  following 

equation: 


Hif  =  Qi  Hip  -  z  (pj  -  pf) 
Pf 
(Lusczynski,  1961)  where: 

i  =  any  point  in  ground  water  of  varible  density 

Hif  =  fresh  water  head  at  i 

Hip  =  point  water  head  at  i 

z  =  elevation  of  i,  measured  positively  upward 

Pi  =  water  density  at  i 

Pf  =  fresh  water  density 


(Dl) 
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KBMP01  (146  ft  deep)  32%  seawater 

Hif  =  (62.93  lb/ft3)  (0.66  ft)  -  (-141.4  ft)  (62.93  lb/ft3  -  62.43  lb/ft3) 

62.43  lb/ft3 

Hif  =  1.80  ft 


KBMP  03  (94  ft  deep)  68%  seawater 

Hif  =  (63.49  lb/ft3)  (0.15  fO  -  (-89.4  ft)  (63.49  lb/ft3  -  62.43  lb/ft3) 

62.43  lb/ft3 

Hif  =  1.67  ft 


KBMP  04  (95  ft  deep)  73%  seawater 

Hif  =  (63.55  lb/ft3)  (3.87  ft)  -  (-88.14  ft)  (63.55  lb/ft3  -  62.43  lb/ft3) 

62.43  lb/ft3 

Hif  =  5.52  ft 


KBMP  05  (44  ft  deep)  3%  seawater 

Hif  =  (62.49  lb/ft3)  (4.41ft)  -  (-37.14  ft)  (62.49  lb/ft3  -  62.43  lb/ft3) 

62.43  lb/ft3 

Hif  =  4.44  ft 


KBMP  06  (71  ft  deep)  89%  seawater 

Hif  =  (63.80  lb/ft3)  (3.83  ft)  -  (-64.14  ft)  (63.80  lb/ft3  -  62.43  lb/ft3) 

62.43  lb/ft3 

Hif  =  5.32  ft 
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KBMP  07  (89  ft  deep)  4.2%  seawater 

Hjf  =  (62.49  lb/ft3)  (2.55  ft)  -  (-73.07  ft)  (62.49  lb/ft3  -  62.43  lb/ft3) 

62.43  lb/ft3 

Hif  =  2.62  ft 


KBMP  1 1  (95  ft  deep)  5.2%  seawater 

Hif  =  (62.49  lb/ft3)  (2.55  ft)  -  (-79.07  ft)  (62.49  lb/ft3  -  62.43  lb/ft3) 

62.43  lb/ft3 

Hif  =  2.62  ft 
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APPENDIX  E 


Leakage  and  Analytical  Model  Calculations 
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Steady  state  vertical  leakage  between  the  shallow  aquifers. 

The  horizontal  hydraulic  conductivity  for  a  clean  clay  based  on  table  2.2  from 
Freeze  and  Cherry  (1979)  is  1.0  x  lO2  ft/d.  The  vertical  hydraulic  conductivity  is 
assumed  to  be  a  half  an  order  of  magnitude  smaller  than  the  horizontal  hydraulic 
conductivity.  Therefore  the  vertical  hydraulic  conductivity  for  the  Pleistocene  aquiclude  is 
assumed  to  be  5.0  x  103  ft/d.  The  vertical  hydraulic  conductivty  for  the  Hawthorn  and 
Charlton  Formations  is  1.0  x  10"3  ft/d  (Brown,  1984).  The  vertical  hydraulic  conductivity 
of  the  confining  layer  between  the  Pliocene  sand  and  the  Miocene  limestone  is  1.7  x  10"2 
ft/d  based  on  aquifer  test  results.  All  water  levels  were  measured  6/28/90. 

The  leakage  equation  is  expressed  as: 

q  =  K'  0„v,  -  0 
b' 

q  =  the  leakage  rate  or  the  volume  of  water  flowing  per  unit  time  through  a  unit  cross- 
sectional  area  [ft/d] 

K'  =  vertical  hydraulic  conductivity  of  the  confining  layer  [ft/d] 
0ext  =  potentiometric  head  of  the  overlying  or  underlying  zone  [ft] 
0  =  potentiometric  head  in  the  aquifer  [ft] 
b1  =  confining  layer  thickness  [ft]  (Bear,  1979). 

Where  the  leakage  rate,  q,  is  negative,  groundwater  is  leaving  the  aquifer.  Where  the 

leakage  rate  is  positive,  groundwater  is  entering  the  aquifer. 

Leakage  from  the  surficial  aquifer  to  the  Pliocene-Miocene  aquifer  at  site  3 
Water  levels  above  mean  sea  level :  KB  MP  8  =  4.93  ft 

KBMP  9  =  2.48  ft 
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q  =  (5.0  x  1Q-3  ft/d)  (2.45  ft) 
30  ft 


q  =  4.08  x  10^  ft/d 

Volumetric  leakage  for  one  square  mile: 

1  mi2  =  2.79  x  107  ft2 

(4.08  x  10-4  ft/d)  (2.79  x  107  ft2)  =  1.14  x  104  ft3/d 

LHxK^ftVd  x  7.481  gal/ft3  =  8.52  x  104 gal/d 

Leakage  from  the  Pliocene-Miocene  aquifer  to  the  surficial  aquifer  at  site  2. 
Water  levels  above  mean  sea  level  :  KBMP  4  =  5.52  ft 

KBMP  5  =  4.44  ft 


q  =  (5.0x  IP3 ft/d)  (-1.08  ft) 
17  ft 

q  =  -3.17X10-4  ft/d 


Leakage  from  the  Pliocene-Miocene  aquifer  to  the  surficial  aquifer  site  1. 
Water  levels  above  mean  sea  level  :  KBMP  3  =  1.67  ft 

KBMP  2  =  0.5  ft 


q  =  (1.0x  10-3  ft/d)  (-1.17  ft) 
56  ft 

q  =  -2.09  x  10-5  ft/d 


Leakage  from  the  Miocene  sand  aquifer  to  the  Pliocene-Miocene  aquifer  at  site  1. 
Water  levels  above  mean  sea  level  :  KBMP  1  =  1.8  ft 


145 


KBMP3  =  1.67ft 


q  =  (1.0  xlO'3  ft/d)  (0.13  m 
42  ft 

q  =  3.1xl0-6ft/d 


Leakage  from  the  Charlton  limestone  (KBMP  7)  to  the  Duplin  sand  (KB MP  9)  at  site  3. 
K'  was  calculated  from  leakance  values  obtained  through  aquifer  testing. 
Water  levels  above  mean  sea  level  :  KBMP  9  =  2.48  ft 

KBMP  7  =  2.62  ft 

q  =  (1.7  x!Q-2  ft/d)  (0.14  ft) 
8  ft 

q  =  2.97  x  10"*  ft/d 

Volumetric  leakage  from  the  Charlton  limestone  to  the  Duplin  sand  for  a  one  square 

mile  area. 

1  mi2  =  2.79  x  107  ft2 

(2.97  x  10"4  ft/d)  (2.79  x  107  ft2)  =  8.29  x  103  ft3/d 

8.29  x  103  ft3/d  x  7.481  gal/ft3  =  6.2  x  10*gal/d 

The  amount  of  recharge  to  the  island  is  an  estimated  value  based  on  results  of 
Brown  (1984).  An  annual  water  budget  for  the  surficial  aquifer  was  not  determined  during 
this  study.  Therefore  the  calculated  recharge  volume  represents  an  approximate  value 
which  is  applied  to  quantitatively  define  the  system  flow. 
Total  recharge  to  the  surficial  aquifer  by  precipitation  for  1  year: 
6  inches  of  recharge  x  1  mi2  =  1.39  x  107  ft3 
1.39  x  107  ft3  x  7.481  gal/ft3  =    104  x  106  gal/mi2 
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Seepage  discharge  from  the  surficial  aquifer: 

q  =  [K  (hi2^_h22)] 
2L 

(RaudMvi  and  Callander,  1976). 

K  =  hydraulic  conductivity  of  the  surficial  aquifer  =  21  ft/d 

hi  =  water  level  at  site  3  =  4.93  ft 

h2  =  water  level  at  the  seepage  face  =  0  ft 

L  =  length  from  center  of  the  island  =  2640  ft 

q  =  21  ft/d  [  (4.93  ft)2  -  (0  ft)2  1 
2  (2640  ft) 

q  =  0.096  ft2/d 

Total  discharge  for  the  distance  of  1  mile: 

1  mile  =  5280  ft 

0.096  ft2/d  x  5280  ft  =  506.88  ft3/d 

506.88  ft3/d  x  7.481  gal/ft3  =  3.79  x  103  gal/d 

The  calculated  specific  discharge  along  the  seepage  face  near  site  3  is  valid  given 

that  a  time-average  mean  is  assumed  where  the  water  levels  are  at  a  constant  altitude. 

Analytical  Model  - 

Recharge  and  discharge  values  were  calculated  for  a  1  square  mile  area. 
Total  recharge  from  precipitation  per  year  =  104  x  106  gal/mi2 
Total  recharge  from  precipitation  per  day  =  2.85  x  105  gal/mi2 
Total  discharge  along  the  seepage  face  3.79  x  103  gal/d  x  2  =  7.58  x  103  gal/d 
Total  discharge  to  the  underlying  Pliocene-Miocene  age  aquifer  =  8.5  x  104  gal/d 
Discharge  to  the  lakes  and  lowlying  wetiands  = 
2.85  x  105  gaVmi2  -  7.58  x  103  gal/d  -  8.5  x  104  gal/d  =  1.92  x  105  gal/d 
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APPENDIX  F 


Step  Drawdown 
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Step  drawdown  test 

The  step-drawdown  test  was  conducted  by  selecting  five  pumping  rates  or  steps 
each  lasting  1.5  hours  and  measuring  drawdown  in  the  pumped  well.  The  data  provided 
from  the  test  were  used  to  determine  the  specific  capacity  of  the  well  and  the  optimum 
pumping  rate. 

The  rates  chosen  for  the  step-drawdown  test  were  based  on  the  thickness  of  the 
aquifer  and  the  estimated  hydraulic  conductivity.  Estimated  hydraulic  conductivity  values 
were  obtained  from  Table  2.2  from  Freeze  and  Cherry  (1979).  The  five  rates  chosen  were: 
1.34,  2.67,  4.01,  5.35  and  6.69  ft3/min;  the  actual  rates  pumped  were  1.47,  3.08,  3.61 
and  4.68  ft3/min.  The  first  three  rates  were  conducted  for  1.5  hours  each.  Thirty  minutes 
into  the  fourth  step  (4.68  ft3/min.)  the  water  level  declined  below  the  intake  hose  causing 
pumping  to  cease  and  the  step-drawdown  test  to  be  discontinued. 

The  data  provided  by  the  three  successful  steps  were  used  to  determine  the  optimum 
pumping  rate  by  comparing  the  specific  capacities  of  each  step.  The  specific  capacity 
decreases  as  the  degree  of  turbulent  head  loss  increases.  Based  on  the  specific  capacities 
calculated  for  the  three  rates,  the  optimum  pumping  rate  for  the  aquifer  test  was  1.47 
ft3/min.  The  rate  used  for  the  aquifer  test  was  1.54  ft3/min  (see  figure  Fib). 

C.E.  Jacob  (1947)  presented  a  method  to  evaluate  the  head  loss  resulting  from 
turbulent  flow  in  an  artesian  well  by  means  of  a  step-drawdown  test.  As  described  by 
Theis  (1935),  drawdown  in  a  confined  aquifer  under  laminar  flow  conditions  is  expressed: 


s=  _Q_   W(u)  (Fl) 

4tcT 
where 


s  =  drawdown  [ft] 
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B  =  3.25  min/ft2  (intercept) 
C  =  0.25  min2 /ft5  (slope) 


0      0.5      1 


1.5      2      2.5 
Q  (ft3  /min) 


3.5 


B. 


Q  (ft^/vi.n) 

s(ft) 

s  (total) 

s/Q  (min/ft2; 

n 

Q/s  (ft$fo,/i) 

1.47 

5.4 

5.4 

3.67 

89% 

0.272 

3.08 

6.78 

12.18 

3.95 

80% 

0.253 

3.61 

2.9 

15.08 

4.17 

78% 

0.239 

Figure  Fl.     A.  Plot  of  s/Q  vs.  Q  to  determine  B  and  C  values.  B.  Table  of  discharge, 
drawdown,  percent  of  head  loss  due  to  laminar  flow,  and  specific  capacity. 
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Q  =  pumping  [fttymin] 

T  =  transmissivity  [ft2/min] 
W(u)  =  well  function 

Equation  (Fl)  can  be  shortened  to 

s  =  BQ  (F2) 


where  B  =  W(u) 
4kT 


Total  drawdown  in  the  discharging  well  includes  both  laminar  flow  and  turbulent  flow 
(Jacob,  1947)  and  is  expressed  as: 

s  =  BQ  +  CQ2  (F3) 

where 

B  =  the  head  loss  in  the  formation  per  unit  discharge  (aquifer  constant)  [min/ft2]; 
C  =  the  head  loss  in  the  well  due  to  turbulence  per  unit  discharge  squared  (well  loss 

constant)  [min2/ft5];  and 
Q  =  the  discharge  rate  [ft3/min].  By  dividing  equation  (F3)  by  Q,  the  equation  becomes 

s/Q  =  CQ  +  B,  (F4) 

which  is  a  linear  equation.  If  s/Q  and  Q  are  plotted  on  arithmetic  graph  paper,  the  slope  of 
the  line  formed  equals  C  and  the  intercept  equals  B  (see  figure  Fla  for  B  and  C  values). 
Once  B  and  C  are  known,  the  percentage  of  total  head  loss  due  to  laminar  flow  can  be 
calculated  by  dividing  the  laminar  head  loss  by  the  total  head  loss 

BO  x    100  (F5) 

BQ  +  CQ2 
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where  n  equals  the  efficiency  of  the  discharging  well  (Driscoll,  1986). 

For  a  pumping  rate  of  1.47  ft3/min.,  89  percent  of  the  head  loss  is  attributable  to 
laminar  flow  (see  figure  Fib  for  n  values  for  each  step).  But  the  percentages  calculated,  n, 
do  not  represent  true  well  efficiency.  To  derive  the  true  well  efficiency,  the  aquifer 
transmissivity  must  be  determined  by  performing  a  constant-rate  test.  Then,  the  well 
efficiency  is  calculated  by  dividing  the  actual  specific  capacity,  (Q/s)Actuai>  by  the  theoretical 
specific  capacity.  The  theoretical  specific  capacity  is  expressed  as: 

(Q/s)TheoreticaI  =  T/2000  (F6) 

(Driscoll,  1986).  The  theoretical  specific  capacity  is  an  empirical  equation  derived  by 
assuming  an  average  well  diameter  (0.5  ft.),  an  average  duration  of  pumping  (24  hrs.),  and 
typical  storage  coefficient  value  for  a  confined  aquifer  (1.0  x  10~3)  (Driscoll,  1986).  The 
well  efficiency  calculated  for  KB  MP  7  at  a  pumping  rate  of  1.54  ft3/min.  is  82  percent 
(Calculations  follow). 
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Step  Drawdown  Data 

Head  loss  in  well  due  to  laminar  flow  (n).  From  Driscoll  (1986),  (n)  is  calculated  from: 


n  =  BO  x    100 

BQ  +  CQ 

From  the  graph  of  Q  vs.  s/Q: 

B  =  3.25  min/ft2 
C  =  0.25  min2/ft5 

ForQ=  1.47ft3/min 

n  = (3.25  min/ft2)  (1.47  ft3/min) =  89% 

(3.25  min/ft2)  (1.47  ft3/min)  +  (0.25  min2/ft5)  (1.47  ft3/min)2 


For  Q  =  3.08  ft3/min 

n  = f3.25  min/ft2-)  (3.08  ft3/min) =  80% 

(3.25  min/ft2)  (3.08  ft3/min)  +  (0.25  min2/ft5)  (3.08  ft3/min)2 


For  Q  =  3.61  ft3/min 

n  = (3.25  min/ft2-)  (3.61  ft3/min) =  78% 

(3.25  min/ft2)  (3.61  ft3/min)  +  (0.25  min2/ft5)  (3.61  ft3/min)2 
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Casing  storage  effects  from  fDriscoll.  1986) 

tc  =  0.6  fcU  -423 
Q/s 

dc  =  diameter  of  casing  =  4  inches 

dp  =  diameter  of  pump/discharge  pipe  =  5/8  inch 

tc  =  0.6  \(4)2  -  (5/8)2l        =    4.1  minutes 
11.5/5.0 

at  4. 1  minutes  s  =  4.68  ft 

l  =  0.6r(4)2-(5/8)2l        =    3.81  minutes 
11.5/4.68 

at  3.81  minutes  s  =  4.62  ft 

tc  =  0.6  TC4)2  -  (5/8)21        =    3.76  minutes 
11.5/4.62 

at  3.76  minutes  s  =  4.62  ft 

At  3.76  minutes  the  casing  storage  effects  become  negligible 

OR 

V  =  7tr2h  =  (3.14)  (0.166  ft)2  (75.83  ft)  =  6.56  ft3   =  49.07  gal 

49.07  gal  =  4.26  minutes 
11.5  gpm 

At  4.26  minutes  the  casing  storage  effects  become  negligible 

Pumping  well  efficiency 

T  =  650  ft2/d  =  4863  gpd/ft 

Total  drawdown  after  371  minutes  =  5.7  ft 
(after  this  time  water  levels  fluctuate  with  the  tides) 

(Q/s)actuai=  11.5  gpm/ 5.7  ft  =  2.0gpm/ft 

(Q/s)theoreticai  =  T/2000  =  4863  gpd/ft  =  2.43  gpm/ft 

2000 

E  =       actual  specific  capacity        =     2.0     =  82%  efficient 
theoretical  specific  capacity         2.43 
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APPENDIX  G 


Derivation  of  the  Diffusion  Equation 
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The  diffusion  equation  is  a  combination  of  the  Darcy's  law  and  the  principle  of 
continuity.  Darcy's  law  is  expressed  as: 

Q  =  -KAAh_  (Gl) 

AL 
where: 

Q  =  discharge  rate  [ft3/d] 
K  =  hydraulic  conductivity  [ft/d] 
A  =  area  [ft2] 
h  =  head  [ft] 
L  =  length  [ft] 

Dupuit  (1848)  made  the  assumption  that  within  the  cone  of  depression  of  a 
pumping  well  the  head  throughout  a  vertical  line  in  the  aquifer  is  constant  and  may  be 
represented  by  the  water  level  at  a  given  point  in  the  aquifer.  If  the  aquifer  is  homogeneous 
and  horizontal,  then  by  the  law  of  continuity  (Qi  -  Q2  =  0)  virtually  equal  quantities  of 
water  are  discharged  from  a  pumping  well  as  are  discharged  through  any  two  concentric 
cylinders  within  the  cone  of  depression.  Under  these  conditions  Darcy's  law  may  be 
expressed  as  a  first  order  differential  equation  where: 


and 


Qi  =  -  2  %  n  K  b  (dh)       [ft3/d]  (G2) 

dri 


Q2  =  -27tr2Kb(dh)       [ft3/d]  (G3) 

dT2 

where 

ri  =  distance  from  pumping  well  to  the  first  concentric  circle  [ft] 

Qi  =  discharge  at  ri  [ft3/d] 

r2  =  distance  from  pumping  well  to  the  second  concentric  circle  [ft] 

Q2  =  discharge  at  r2  [ft3/d] 

b  =  saturated  thickness  of  aquifer  [ft] 

For  steady  state  flow: 

Qi    -   02=  0  (G4) 

For  transient  flow: 

Qi    -    Q2=  SAah  (G5) 

3t 
S  =  Storage  coefficient  [dimensionless] 
combining  equations  (G2)  and  (G3)  yields: 

Q1-Q2  =  -2  7tKb(rAi-r2(^)2) 

dr  dr  (G6) 
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because  Kb  =  T  equation  (G6)  may  be  expressed: 


Q1-Q2  =  -2,tT((r|)2-(r^)1) 


(G7) 


rearranging  the  equation  in  general  form  becomes: 


Qi    -   02=   2tcT   d(rdh/dr)  Ar  (G8) 

9r 

by  substituting  equation  (G5)  into  equation  (G8)  results  in: 

S27tr8h=    2juT  (r92h  +  ah}  (G9) 

9t  dr2     dr 

dividing  both  sides  of  the  equation  by  27iTr  results  in  the  partial  differential  equation  for 
nonsteady  flow  or  the  diffusion  equation: 

l9h_+a2h  =  S9h  (G10) 

r9r     9r2     T3t 
(Lohman,  1972). 

To  solve  the  diffusion  equation  for  nonsteady  radial  flow  boundary  conditions  must 
be  satisfied  and  are  as  follows: 

h  (r,t) 

h(r,0)  =  h0  for   r>0 

h(«»,t)  =  h0  for  t  >  0 

If  the  rate  of  pumping  is  constant  and  the  boundary  conditions  are  applied,  the  equation  (in 
SI  units)  may  be  expressed  as: 


s  =  _Q_f  e^du 
4tcTJ     u 


(Gil) 


where 

u  =  i2S 
4Tt 

expanding  equation  (Gil) 


l-u+ui.ul 

s  =  4^r|         2!u3'       d" 


(G12) 


dividing  by  u  gives 
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s  = 


d»~ 


— '    -     l+£  -^-...du 


4tcT|    u  2!      3! 

by  integrating,  the  equation  becomes 

s=z!^iogeu-u+— -— •••tf 

47cT  2-2!     3-3! 

by  applying  the  limits  and  solving  equation  (G14)  results  in 

O  2  3 

s  =  -r-—  [-  Euler  constant  -  (logeu  -  u  +  — —  .  .  .  ] 

4llT  2-2!     3-3! 


(G13) 
(G14) 


(G15) 


or 

s 


(G16) 


=  ^L[_.5772l6-logeu  +  u--^  +  -^-...] 
4ttT  2-2!      3-3! 

for  values  of  u  <  0.01  equation  (G16)  may  be  terminated  to  the  second  term: 

s  =  _Q_  [-0.577216  -  loge  r2S  ]  (G17) 

47tT  4Tt 

or 

s  =  _Q_  [-loge     4       +  loge  4H]  (G 1 8) 

47tT  2.25  r2S 

by  converting  to  common  log  the  equation  becomes: 

T=  2.30O     login  2.25Tt  (G19) 

47ts  r2S 

in  equation  (G19)  T,  Q,  S,  and  r  are  constants  and  the  equation  may  be  written  : 

T  =  2.30O     [login  2.25Tr   +  login  t]  (G20) 

47ts  S 


by  differential  calculus  the  equation  becomes: 

T  =       2.30O  (G21) 

47tds/dlogio  t 

By  changing  from  derivatives  to  finite  values  the  equation  becomes: 

T  =       2.30O  (G22) 

47tAs/Alogio  t 
(Lohman,  1972)  where 

As/Alogio  t  =  the  slope  of  the  line  (m)  on  semilogarithmic  paper.  Equation  (G22)  may  be 
expressed  as  (Bixel  and  Van  Poollen,  1967): 

T  =       2.30O  (G23) 

47tm 
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APPENDIX  H 


Aquifer  Test  Data  and  Calculations  for  Transmissivity  and  Storage  Values 
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Static  water  level  =  2.76  ft  above  sea  level 


Time,  in       Drawdown,     "Drawdown,                       Time,  in       Drawdown,      "Drawdown, 
minutes           in  feet              in  feet                            minutes           in  feet              in  feet 

0.5 

1.34 

1.34 

191 

5.27 

5.54 

1 

3.65 

3.65 

221 

5.29 

5.58 

1.5 

4.14 

4.14 

251 

5.32 

5.62 

2 

4.35 

4.35 

281 

5.37 

5.66 

3 

4.55 

4.55 

311 

5.42 

5.69 

4 

4.68 

4.68 

371 

5.51 

5.7 

5 

4.8 

4.8 

431 

6.12 

6.16 

6 

4.86 

4.86 

491 

6.23 

6.13 

7 

4.93 

4.93 

551 

6.29 

6.07 

8 

4.99 

4.99 

611 

6.32 

6.06 

9 

5.01 

5.01 

671 

6.3 

6.12 

10 

5.06 

5.06 

731 

6.19 

6.15 

15 

5.14 

5.14 

791 

6.08 

6.18 

20 

5.24 

5.24 

851 

5.89 

6.13 

25 

5.25 

5.25 

911 

5.71 

6.06 

30 

5.28 

5.28 

971 

5.57 

6 

35 

5.3 

5.3 

1031 

5.53 

5.99 

40 

5.31 

5.31 

1151 

5.76 

6.08 

45 

5.31 

5.31 

1271 

5.92 

6.11 

50 

5.32 

5.32 

1391 

6.15 

6.06 

55 

5.31 

5.31 

1514 

6.09 

6.03 

60 

5.33 

5.33 

1574 

6.02 

6.1 

65 

5.34 

5.34 

1601 

5.68 

5.82 

70 

5.33 

5.33 

1721 

5.54 

5.83 

75 

5.34 

5.34 

1841 

5.66 

5.93 

80 

5.35 

5.35 

1961 

5.92 

5.89 

85 

5.33 

5.33 

2301 

6.14 

6.07 

90 

5.36 

5.36 

2391 

5.86 

6.15 

95 

5.34 

5.34 

2514 

5.61 

6.04 

100 

5.32 

5.46 

3634 

5.84 

6.26 

130 

5.3 

5.49 

2754 

6.15 

6.14 

160 

5.3 

5.54 

2814 

6.19 

6.04 

*  Drawdown  data  corrected  for  tidal  and  barometric  fluctuations 


Table  HI .  Drawdown  data  for  KBMP  7,  December  8,  1989. 


Water  level  at  maximum  drawdown  =  3.28  ft  below  sea  level 


Time,  in        Recovery,      "Recovery,                             Time,  in         Recovery,      "Recovery, 
minutes            in  feet             in  feet                                 minutes            in  feet             in  feet 

1 

4.16 

4.16 

80 

5.61 

5.61 

1.5 

4.39 

4.39 

85 

5.63 

5.63 

2 

4.52 

4.52 

90 

5.64 

5.64 

2.5 

4.61 

4.61 

95 

5.65 

5.65 

3 

4.69 

4.69 

100 

5.67 

5.67 

3.5 

4.75 

4.75 

130 

5.74 

5.74 

4 

4.8 

4.8 

150 

5.79 

5.79 

4.5 

4.85 

4.85 

180 

5.87 

5.79 

5 

4.91 

4.91 

210 

5.94 

5.79 

6 

4.97 

4.97 

240 

6.03 

5.79 

7 

5.03 

5.03 

270 

6.11 

5.82 

8 

5.07 

5.07 

300 

6.17 

5.85 

9 

5.11 

5.11 

360 

6.2 

5.87 

10 

5.15 

5.15 

420 

6.13 

5.82 

11 

5.17 

5.17 

480 

6 

5.8 

12 

5.2 

5.2 

540 

5.91 

5.89 

14 

5.23 

5.23 

600 

5.84 

5.99 

15 

5.25 

5.25 

660 

5.79 

6.08 

20 

5.34 

5.34 

720 

5.77 

6.11 

25 

5.39 

5.39 

780 

5.78 

6.04 

30 

5.42 

5.42 

840 

5.87 

5.98 

35 

5.46 

5.46 

900 

6.03 

5.97 

40 

5.49 

5.49 

960 

6.24 

6.01 

45 

5.5 

5.5 

1020 

6.44 

6.07 

50 

5.52 

5.52 

1080 

6.57 

6.13 

55 

5.55 

5.55 

1140 

6.55 

6.1 

60 

5.56 

5.56 

1260 

6.24 

6.02 

65 

5.58 

5.58 

1380 

6.03 

6.16 

70 

5.59 

5.59 

1500 

5.94 

6.23 

75 

5.61 

5.61 

1700 

5.66 

5.73 

*  Recovery  data  corrected  for  tidal  and  barometric  fluctuations 


Table  H2.  Recovery  data  for  KBMP  7,  December  10,  1989. 
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Static  water  level  =  3.05  ft  above  sea  level 


Time,  in   Drawdown,  'Drawdown,         Time,  in   Drawdown,  'Drawdown, 
minutes    in  feet      in  feet           minutes    in  feet     in  feet 

1 

0.01 

0.01 

220 

0.97 

1.22 

2 

0.18 

0.18 

250 

0.98 

1.23 

3 

0.31 

0.31 

280 

1.01 

1.25 

4 

0.4 

0.4 

340 

1.05 

1.21 

5 

0.5 

0.5 

400 

1.17 

1.25 

6 

0.56 

0.56 

460 

1.24 

1.19 

7 

0.6 

0.6 

520 

1.35 

1.18 

8 

0.64 

0.64 

580 

1.41 

1.19 

9 

0.7 

0.7 

640 

1.42 

1.19 

10 

0.73 

0.73 

700 

1.37 

1.24 

15 

0.81 

0.81 

760 

1.29 

1.28 

20 

0.88 

0.88 

820 

1.19 

1.31 

25 

0.91 

0.91 

880 

0.99 

1.23 

30 

0.95 

0.95 

940 

0.82 

1.16 

35 

0.96 

0.96 

1000 

0.73 

1.11 

40 

0.99 

0.99 

1120 

0.79 

1.13 

45 

0.97 

0.97 

1240 

1.07 

1.17 

50 

1.01 

1.01 

1360 

1.22 

1.1 

55 

1.02 

1.02 

1480 

1.27 

1.2 

60 

1.02 

1.02 

1600 

1.16 

1.28 

65 

1.02 

1.02 

1720 

1.01 

1.27 

70 

1.02 

1.02 

1840 

1.1 

1.31 

75 

1.02 

1.02 

1960 

1.31 

1.29 

80 

1.02 

1.02 

2080 

1.45 

1.2 

85 

1 

1.12 

2200 

1.49 

1.29 

90 

0.99 

1.13 

2320 

1.27 

1.31 

95 

0.99 

1.12 

2440 

0.95 

1.27 

100 

0.99 

1.14 

2560 

0.86 

1.24 

130 

0.99 

1.2 

2680 

1.1 

1.28 

160 

0.98 

1.21 

2800 

1.3 

1.24 

190 

0.96 

1.2 

*  Drawdown  data  corrected  for  tidal  and  barometric  fluctuations 


Table  H3.  Drawdown  data  for  KBMP  11,  December  8, 1989. 


Water  level  at  maximum  drawdown  =  1.81  ft  above  sea  level 


Time,  in       Recovery,       'Recovery,                        Time,  in       Recovery,        "Recovery, 
minutes         in  feet              in  feet                           minutes          in  feet               in  feet 

1 

0.14 

0.14 

95 

1.18 

1.18 

2 

0.28 

0.28 

100 

1.19 

1.26 

3 

0.4 

0.4 

120 

1.22 

1.27 

4 

0.48 

0.48 

150 

1.28 

1.26 

5 

0.54 

0.54 

180 

1.34 

1.27 

6 

0.59 

0.59 

210 

1.4 

1.28 

7 

0.64 

0.64 

240 

1.48 

1.3 

8 

0.69 

0.69 

270 

1.59 

1.37 

9 

0.72 

0.72 

360 

1.69 

1.42 

10 

0.74 

0.74 

420 

1.61 

1.34 

15 

0.85 

0.85 

480 

1.51 

1.34 

20 

0.91 

0.91 

540 

1.44 

1.42 

25 

0.95 

0.95 

600 

1.34 

1.47 

30 

0.99 

0.99 

660 

1.32 

1.57 

35 

1 

1 

720 

1.3 

1.6 

40 

1.01 

1.01 

780 

1.31 

1.55 

45 

1.01 

1.01 

840 

1.39 

1.49 

55 

1.08 

1.08 

900 

1.46 

1.44 

60 

1.09 

1.09 

960 

1.86 

1.53 

65 

1.1 

1.1 

1020 

1.99 

1.59 

70 

1.11 

1.11 

1080 

1.96 

1.57 

75 

10.1 

10.1 

1140 

1.73 

1.54 

80 

1.13 

1.13 

1260 

1.54 

1.62 

85 

1.13 

1.13 

1334 

0.99 

90 

1.15 

1.15 

Recovery  data  corrected  for  tidal  and  barometric  fluctuations 


Table  H4.  Recovery  data  for  KBMP  11,  December  10,  1989. 
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Static  water  level  =  2.67  ft  above  sea  level 


Time,  in      Drawdown,      "Drawdown,                       Time,  in       Drawdown,      "Drawdown, 
minutes          in  feet              in  feet                           minutes           in  feet              in  feet 

1 

0.2 

0.2 

281 

0.95 

1.19 

2 

0.21 

0.21 

341 

1.09 

1.25 

3 

0.25 

0.25 

401 

1.11 

1.17 

4 

0.33 

0.33 

461 

1.2 

1.15 

4.5 

0.39 

0.39 

521 

1.31 

1.14 

5 

0.46 

0.46 

581 

1.37 

1.15 

5.5 

0.5 

0.5 

641 

1.37 

1.14 

6 

0.54 

0.54 

701 

1.32 

1.19 

7 

0.57 

0.57 

761 

1.21 

1.2 

7.5 

0.6 

0.6 

821 

1.1 

1.22 

8 

0.65 

0.65 

881 

0.92 

1.16 

9 

0.67 

0.67 

941 

0.76 

1.1 

10 

0.68 

0.68 

1001 

0.66 

1.04 

15 

0.81 

0.81 

1121 

0.73 

1.07 

20 

0.88 

0.88 

1241 

1.05 

1.15 

25 

0.91 

0.91 

1361 

1.2 

1.08 

30 

0.93 

0.93 

1481 

1.23 

1.16 

35 

0.95 

0.95 

1601 

1.1 

1.22 

40 

0.96 

0.96 

1721 

0.97 

1.23 

45 

0.97 

0.97 

1841 

1.07 

1.28 

50 

0.97 

0.97 

1961 

1.26 

1.26 

55 

0.98 

0.98 

2081 

1.18 

1.18 

86 

0.97 

0.97 

2201 

1.23 

1.23 

101 

0.91 

1.06 

2321 

1.17 

1.17 

131 

0.93 

1.13 

2441 

1.23 

1.23 

161 

0.9 

1.13 

2561 

1.18 

1.18 

191 

0.86 

1.11 

2681 

1.23 

1.23 

221 

0.86 

1.12 

2801 

1.16 

1.16 

251 

0.9 

1.15 

*  Drawdown  data  corrected  for  tidal  and  barometric  fluctuations 


Table  H5.  Drawdown  data  for  KBMP  12,  December  8,  1989. 
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Water  level  at  maximum  drawdown  =  1.51  ft  above  sea  level 


Time,  in      Recovery,     'Recovery,                       Time,  in       Recovery,      'Recovery, 
minutes          in  feet           in  feet                           minutes          in  feet            in  feet 

1 

0.09 

0.09 

74 

1.17 

1.17 

2 

0.18 

0.18 

79 

1.19 

1.19 

3 

0.25 

0.25 

84 

1.19 

1.19 

4 

0.31 

0.31 

89 

1.21 

1.21 

4.5 

0.38 

0.38 

94 

1.23 

1.23 

5 

0.42 

0.42 

99 

1.23 

1.28 

5.5 

0.48 

0.48 

120 

1.28 

1.32 

6 

0.54 

0.54 

150 

1.34 

1.33 

6.5 

0.56 

0.56 

180 

1.41 

1.34 

7 

0.6 

0.6 

210 

1.49 

1.36 

7.5 

0.62 

0.62 

240 

1.57 

1.39 

8 

0.67 

0.67 

270 

1.64 

1.41 

8.5 

0.67 

0.67 

300 

1.67 

1.41 

9 

0.72 

0.72 

360 

1.75 

1.45 

10 

0.75 

0.75 

420 

1.7 

1.43 

10.5 

0.75 

0.75 

480 

1.59 

1.42 

11 

0.8 

0.8 

540 

1.49 

1.47 

14 

0.86 

0.86 

600 

1.42 

1.55 

19 

0.92 

0.92 

660 

1.38 

1.63 

24 

0.98 

0.98 

720 

1.36 

1.65 

29 

1.01 

1.01 

780 

1.35 

1.65 

34 

1.03 

1.03 

840 

1.41 

1.51 

39 

1.07 

1.07 

900 

1.54 

1.49 

44 

1.08 

1.08 

960 

1.92 

1.59 

49 

1.1 

1.1 

1020 

2.08 

1.68 

54 

1.12 

1.12 

1080 

2.07 

1.67 

59 

1.14 

1.14 

1140 

1.83 

1.64 

64 

1.14 

1.14 

1260 

1.61 

1.52 

69 

1.16 

1.16 

1330 

1.54 

1.54 

Recovery  data  corrected  for  tidal  and  barometric  fluctuations 


Table  H6.  Recovery  data  for  KBMP  12,  December  10, 1989. 
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Calculations  for  transmissivity  and  storage  values  of  the  Pliocene-Miocene  aquifer  using 
the  Hantush-Jacob  leaky  method  (1955). 


T  =SL-  W(u) 

4tis 


S  =  4Tt&2 
1/u 

K'/b"  =  4Tv2/r2 


KBMP  1 1  Drawdown 

T  =  1.54ft3/min       (1.0) 
4  7t(0.51  ft) 

T  =  0.240  ft2/min 

T  =  346ft2/d 

S  =  4  (346  ft2/d)  (9.38x10^  dV(247  ft)2 

1.0/1.0 

S  =  2.12xl0-5 

K'/b'  =  4  (346  ft2/d)  (0.2)2 

(247  ft)2 

K'/b'  =  9.07X10-4  ft/d/ft 


KBMP  1 1  Recovery 

T  =  1.54ft3/min       (1.0) 
4  %  (0.35  ft) 

T  =  0.350  ft2/min 

T  =  504ft2/d 

S  =  4(504ft2/d)  (4.38x10^  dV(247  ft)2 

1.0/1.0 

S  =   1.44xl0-5 

K'/b'  =  4  (504  ft2/d)  (0.1)2 

(247  ft)2 
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K'/b'  =  3.3X10-4  ft/d/ft 


KBMP  12  Drawdown 

T  =  1.54ft3/min       (1.0) 
4  k  (0.75  ft) 

T  =  0.163  ft2/min 

T  =  235  ft2/d 

S  =  4  (235  ft2/d)  (1.58xl0-3  dV(168  ft)2 

1.0/1.0 

S  =  5.27xlO-5 

K'/b'  =  4  (235  ft2/d)  (0.35)2 

(168  ft)2 

K'/b'  =  4.08xlO-3  ft/d/ft 


KBMP  12  Recovery 

T  =  1.54ft3/min       (1.0) 
4  jc  (0.65  ft) 

T  =  0.189ft2/min 

T  =  271  ft2/d 

S  =  4  (271  ft2/d)  (1.46xl0-3  dVf  1 68  ft)2 

1.0/1.0 

S  =  5.6xl0"5 

K'/b'  =  4  (271  ft2/d)  (0.24)2 

(168  ft)2 

K'/b'  =  2.21x10-3  ft/d/ft 
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Straight  line  calcualtions: 

Transmissivity  calculated  from  second  straight  line  (Bixel  and  Van  Poollen,  1967) 
Ti  =  2.3  O 
4ns 

Transmissivity  calculated  from  third  straight  line 

T2  =  Xi  QiSki  - m'^ 
m 


Diffusivity  (Streltsova,  1988) 

r2_1 

2.246  (to) 


rj  = r2 


Storage  (Freeze  and  Cherry,  1979) 
S  =  _T 

KB  MP  11  drawdown 
2nd  curve 


Tl  =  2.3  x  1.54  ft2/min     =     0.353  ft2/min    =  507  ft2/d 
4  %  0.8  ft 

K  =  J_  =  507  ft2/d  =  72.4  ft/d 
b  7ft 

Vs  =  Ki  =  72.4  ft/d  (5.52  ft  -  2.59  ft)/  10.660  ft    =    0.066  ft/d 
T|  0.30 

ri  (247  ft)2  =    21730.7  ft2/min 

2.246  (1.25  min) 


S  = 507  ft2/d =  1.62xl0-5 

21730.7  ft2/min  x  1440  min/d 

3rd  curve 

T2  =  507  ft2/d  (7(0.8^-0.33)     =    1951.2  ft2/d 
0.33 
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KBMP  1 1  Recovery 

2nd  curve 

Tj  =  2.3  x  1.54ft2/min     =     0.454  ft2/min    =  653  ft2/d 
4tc  0.62  ft 


K  =    653  ft2/d  =  93.2  ft/d 
7  ft 

Vs  =  Ki  =  93.2  ft/d  (5.52  ft  -  2.59  ft)/  10.660  ft    =    0.085  ft/d 
ti  0.30 

ti  (247  ft)2  =    43116.5  ft2/min 

2.246  (0.63  min) 

S  = 653  ft2/d =  1.05xl0-5 

431 16.5  ft2/min  x  1440  min/d 

3rd  curve 

T2  =  653ft2/d(2(0.62)-0.38)     =    1477.8  ft2/d 
0.38 


KBMP  12  drawdown 
2nd  curve 


Ti  =  2.3  x  1.54ft2/min     =     0.303  ft2/min    =  436  ft2/d 
4  ti  0.93  ft 


K  =  436  ft2/d  =  62.3  ft/d 
7  ft 

Vs  =  Ki  =  62.7  ft/d  (5.52  ft  -  2.59  ft)/  10.660  ft    =    0.057  ft/d 
ti  0.30 

Ti  (168  ft)2  =    7180.8  ft2/min 

2.246  (1.75  min) 

S  = 436  ft2/d =  4.2  x  105 

7180.8  ft2/min  x  1440  min/d 
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3rd  curve 

T2  =  436  ft2/d  (2  (0.93)  -  0.33)     =    2021  ft2/d 
0.33 

KBMP  12  recovery 

2nd  curve 

Ti  =  2.3  x  1.54ft2/min     =     0.264  ft2/min    =  379.5  ft2/d 
4  7C  1.07  ft 

K  =  379.5  ft2/d  =  54.2  ft/d 
7ft 

Vs  =  Ki  =  54.2  ft/d  (5.52  ft  -  2.59  ft)/  10.660  ft    =    0.049  ft/d 
ti  0.30 

Ti  (168  ft)2  =     5983.9  ft2/min 

2.246  (2.1  min) 

S  = 379.5  ft2/d =  4.4  x  10-5 

5983.9  ft2/min  x  1440  min/d 

3rd  curve 

T2  =  379.5  ft2/d  (2  (1.07)  -0.45)     =    1425  ft2/d 
0.45 


KBMP  7  drawdown 


Ti  =  2.3  x  1.54  ft2/min     =    644  ft2/d 
4  k  0.63  ft 


KBMP  7  recovery 


Ti  =  2.3  x  1.54ft2/min     =    654  ft2/d 
47t  0.62  ft 


Streltsova  (1988)  straight  line  method 
Transmissivity  calculated  from  second  straight  line 
Ti  =  2.3  O 

4tc  s 
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Transmissivity  calculated  from  third  straight  line 
T2  =  Ti  (2rm   -m') 
m' 


KBMP  1 1  Drawdown 


Ti  =  2.3  x  1.54ft2/min     =     0.353  ft2/min    =  507  ft2/d 
4  n  0.8  ft 

T2  =  507  ft2/d  (0.8  -  0.33)  =  722  ft2/d 
0.33 


KBMP  1 1  Recovery 

Ti  =  2.3  x  1.54  ft2/min     =     0.454  ft2/min    =  653  ft2/d 
4tc  0.62  ft 

T2  =  653  ft2/d  (0.62  -  0.38)  =  412  ft2/d 
0.38 


KBMP  12  Drawdown 

Ti  =  2.3  x  1.54ft2/min     =    0.303  ft2/min    =  436  ft2/d 
4  31  0.93  ft 

T2  =  436  ft2/d  (0.93  -  0.33)  =  792  ft2/d 
0.33 


KBMP  12  Recovery 

Ti  =  2.3  x  1.54ft2/min     =     0.264  ft2/min    =  379.5  ft2/d 
4  7C  1.07  ft 

T2  =  380  ft^/d  (1.07 -0.45)  =  522  ft2/d 
0.45 
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APPENDIX  I 


Tidal  Effects 
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Tidal  Efficiency 

Water  levels  in  wells  penetrating  confined  aquifers  respond  to  nearby  fluctuating 
tides  due  to  the  compression  and  expansion  of  the  overlying  confining  layers  as  the  load  of 
the  tidal  water  changes.  Well  tidal  efficiency  is  the  effectiveness  of  a  well's  response  to 
fluctuating  tides  measured  in  percentage.  The  tidal  efficiency  decreases  with  the  depth  of  a 
well  and  the  distance  from  the  tidal  body.  Well  construction  and  the  hydrogeology  of  the 
screened  interval  also  are  factors  in  the  degree  of  tidal  efficiency  (Gregg,  1966).  A  poorly 
constructed  well  in  a  low  permeable  zone  will  have  a  low  tidal  efficiency. 

The  tides  on  Cumberland  Island  are  nearly  semidiurnal  and  range  from  6  to  8  feet 
They  are  measured  by  a  continuous  tide-stage  recorder  operated  by  the  USGS.  The 
recorder  is  located  at  the  Dungeness  dock  on  the  west  side  of  Cumberland  Island. 

Along  the  west  coast,  the  island  is  affected  by  channel  tides  and  on  the  east  coast  by 
ocean  tides.  These  tides  probably  do  not  have  the  same  magnitude  at  the  same  time  which 
causes  asymmetrical  water-level  changes  in  the  confined  aquifers.  This  difference  in  tide 
stage  could  not  be  accounted  for  in  calculations  since  tidal  data  are  not  available  on  the 
ocean  side  of  the  island. 

The  tidal  efficiency  equation  takes  into  account  the  barometric  pressure  change  and 
is  expressed  as: 

T.E.  =  barometric  change  -  water  level  change  (II) 

barometric  change  -  tidal  change 

(Gregg,  1966). 

The  tidal  efficiency  was  calculated  for  one-half  of  a  tide  cycle  for  each  of  the  wells. 

The  change  in  tide  was  corrected  to  its  freshwater  equivalent  by  multiplying  the  tidal  change 

by  the  seawater  density  (1.025  gm/cm3)  and  dividing  by  the  freshwater  density  (1.0 

gm/cm3)  as  shown  below: 
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tidal  change  =  uncorrected  tidal  change  (ft)  x  1.025  ft.  seawater/unit  weight       (12) 

1.000  ft.  freshwater/unit  weight 

(Gregg,  1966). 

The  drawdown  and  recovery  data  for  each  of  the  monitoring  wells  does  not  form  a 
smooth  curve  or  straight  line  after  about  100  min.  from  the  initial  pumping  or  recovery  time 
due  to  the  influence  of  tidal  fluctuations.  During  the  first  100  min.  of  drawdown  and 
recovery,  the  data  were  not  corrected  for  tidal  fluctuation  because  the  magnitude  of  the 
drawdown  and  recovery  had  a  greater  impact  on  water  levels  than  did  the  tidal  effects.  The 
data  were  corrected  for  tidal  fluctuations  by  calculating  the  tidal  efficiency;  multiplying  this 
value  by  the  tide  level;  then  subtracting  or  adding  the  new  value  to  the  water  level, 
depending  on  high  or  low  tidal  fluctuations,  respectively.  This  technique  minimized  the 
tidal  fluctuations  in  the  data,  but  did  not  entirely  alleviate  them.  This  is  partially  due  to  the 
variation  of  lag  time  from  one-half  of  a  tide  cycle  to  the  next  half  of  a  tide  cycle  (high  to  low 
tide  or  vice  versa)  ranging  from  0.5  to  1.5  hours.  The  time  lag  was  "ignored"  during  tidal 
efficiency  calculations  by  matching  the  high  and  low  ocean  tides  with  the  high  and  low 
water  levels  by  a  shift  in  time.  Variations  in  lag  time  made  it  difficult  to  choose 
corresponding  tide  levels  between  half  tides.  Also,  because  the  tidal  fluctuation  data  were 
recorded  in  15-minute  intervals,  the  corresponding  tide  level  could  not  be  matched  exactly 
to  the  water  level  for  a  given  time. 

The  tidal  efficiencies  calculated  for  KBMP  7  (9.4  percent),  KBMP  1 1  (7.9 
percent),  and  KBMP  12  (8.2  percent)  were  in  close  agreement,  which  is  logical  since  these 
wells  are  in  the  same  zone  and  are  located  relatively  close.  The  tidal  efficiency  for  KBMP  9 
(7.2  percent)  is  located  close  to  those  wells  tapping  the  limestone  zone.  Although  this  well 
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penetrates  a  different  type  of  lithology,  the  well  construction  is  similar  to  KB  MP  7  and  is 
located  close  to  the  previously  mentioned  wells. 

Well  KBMP  10  has  a  tidal  efficiency  of  1.4  percent.  Although  it  is  located  in  the 
same  area  as  the  other  wells  and  the  well  construction  is  similar,  the  hydrogeology  is  much 
different.  The  well  taps  a  zone  that  has  a  low  hydraulic  conductivity.  Therefore, 
propagation  of  tidal  fluctuations  is  not  as  effective  in  this  zone. 
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The  tidal  efficiency  is  calculated  by: 

T.E.  =  barometric  change  -  water  level  change 
(13) 
barometric  change  -  tidal  change 

The  tidal  efficiency  was  calculated  for  one-half  of  a  tide  cycle  for  each  of  the  wells. 
The  change  in  tide  was  corrected  to  its  freshwater  equivalent  by  multiplying  the  tidal  change 
by  the  seawater  density  (1.025  g/cm3)  and  dividing  by  the  freshwater  density  (1.0  g/cm3) 
as  shown  below: 

tidal  change  =  uncorrected  tidal  change  (ft)  x  1.025  ft.  seawater/unit  weight       (14) 

1.000  ft.  freshwater/unit  weight 

(Gregg,  1966). 

The  tide  intensity  and  lag  varied  from  each  tide  cycle  which  caused  the  calculated 
tidal  efficiency  to  vary  for  a  well  for  each  tide  cycle.  The  tidal  efficiency  was  calculated  for 
each  well  for  3  or  4  half  tide  cycles.  The  average  value  was  assumed  for  the  well  tidal 
efficiency. 

KBMP  7  tidal  efficiency: 
8/12/89-9/12/90 

change  in  barometric  pressure  =  0  ft. 

change  in  water  level  =  (6.32  ft.  -  5.53  ft.)  =  0.79  ft. 

change  in  tide  level  =  (-2.75  ft.  +  4.92  ft.)  =  7.67  ft.  x  1.025  ft.  seawater/unit 

weight  =  7.86 

0.79  ft.  =   10.1% 
7.86  ft. 
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11/12/90 

change  in  barometric  pressure  =  0  ft. 

change  in  water  level  =  (19.54  ft.  -  20.35  ft.)  =  0.81  ft. 

change  in  tide  level  =  (-3.72  ft.  +  4.86  ft.)  =  8.58  ft.  x  1.025  ft.  seawater/unit 

weight  =  8.79 

0.81  ft.  =  9.2% 
8.79  ft. 

10/12/90-11/12/90 

change  in  barometric  pressure  =  0  ft. 

change  in  water  level  =  (5.53  ft.  -  6.15  ft.)  =  0.62  ft. 

change  in  tide  level  =  (-1.76  ft.  +  4.92  ft.)  =  6.68  ft.  x  1.025  ft.  seawater/unit 

weight  =  6.85 

0.62  ft.  =  9.1% 
6.85  ft. 

The  average  tidal  efficiency  for  KB  MP  7  is  9.4%. 


KBMP  11  tidal  efficiency: 
8/12/90-9/12/90 

change  in  barometric  pressure  =  0.01  ft. 

change  in  water  level  =  0.69  ft. 

change  in  tide  level  =  7.67  ft.  x  1.025  ft.  seawater/unit  weight  =  7.86 

0.01  -  0.69  ft.  =  8.7% 
0.01  -  7.86  ft. 
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9/12/90-10/12/90 

change  in  barometric  pressure  =  0.01  ft 

change  in  water  level  =  0.63  ft. 

change  in  tide  level  =  8.35  ft  x  1.025  ft.  seawater/unit  weight  =  8.56 

0.01  -  0.63  ft.  =  7.3% 
0.01  -  8.56  ft. 

9/12/90 

change  in  barometric  pressure  =  0  ft. 

change  in  water  level  =  0.51  ft 

change  in  tide  level  =  6.68  ft  x  1.025  ft  seawater/unit  weight  =  6.85 

0.51  ft.  =  7.5% 
6.85  ft 

11/12/90 

change  in  barometric  pressure  =  0  ft. 

change  in  water  level  =  0.69  ft 

change  in  tide  level  =  8.58  ft  x  1.025  ft.  seawater/unit  weight  =  8.79 

0.69  ft  =  7.8% 
8.79  ft. 

The  average  tide  efficiency  for  KBMP  1 1  is  7.8%. 

KBMP  12  tide  efficiency  8/12/90-9/12/90: 

change  in  barometric  pressure  =  0.01  ft 

change  in  water  level  =  0.7 1  ft 

change  in  tide  level  =  7.67  ft  x  1.025  ft.  seawater/unit  weight  =  7.86 

0.01  -0.71  ft  =  8.9% 
0.01  -  7.86  ft. 
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9/12/90 

change  in  barometric  pressure  =  0  ft. 

change  in  water  level  =  0.57  ft. 

change  in  tide  level  =  6.68  ft.  x  1.025  ft.  seawater/unit  weight  =  6.8  ft. 

0.57  ft.   =  8.4% 
6.8  ft. 

11/12/90 

change  in  barometric  pressure  =  0  ft. 

change  in  water  level  =  0.73  ft. 

change  in  tide  level  =  8.58  ft.  x  1.025  ft.  seawater/unit  weight  =  8.79  ft. 

0.73  ft.  =  8.3% 
8.79  ft. 

The  average  tide  efficiency  for  KBMP12  is  8.5% 

The  tide  efficiency  for  KBMP  9: 
8/12/90-9/12/90 

change  in  barometric  pressure  =  0.01  ft 

change  in  water  level  =  0.52  ft 

change  in  tide  level  =  7.67  ft.  x  1.025  ft.  seawater/unit  weight  =  7.86  ft. 

0.01  -  0.52  ft.  =  6.5% 
0.01  -  7.86  ft. 

9/12/90 

change  in  barometric  pressure  =  0  ft. 

change  in  water  level  =  0.53  ft. 

change  in  tide  level  =  6.68  ft.  x  1.025  ft.  seawater/unit  weight  =  6.8  ft. 

0.53  ft.   =  7.8% 
6.8  ft. 


179 


9/12/90-10/12/90 

change  in  barometric  pressure  =  0.01  ft. 

change  in  water  level  =  0.67  ft. 

change  in  tide  level  =  8.35  ft.  x  1.025  ft  seawater/unit  weight  =  8.56 

0.01  -  0.67  ft.  =  7.7% 
0.01  -  8.56  ft. 

11/12/90 

change  in  barometric  pressure  =  0  ft. 

change  in  water  level  =  0.60  ft. 

change  in  tide  level  =  8.58  ft.  x  1.025  ft.  seawater/unit  weight  =  8.79  ft. 

0.60  ft.  =  6.8% 
8.79  ft. 

The  average  tide  efficiency  for  KBMP  9  is  7.2%. 

The  tide  efficiency  for  KBMP  10  8/12/90-9/12/90: 

change  in  barometric  pressure  =  0.01  ft 

change  in  water  level  =  0.19  ft. 

change  in  tide  level  =  7.67  ft.  x  1.025  ft.  seawater/unit  weight  =  7.86  ft. 

0.01  -0.19  ft.  =  2.3% 
0.01  -  7.86  ft. 

9/12/90 

change  in  barometric  pressure  =  0  ft. 

change  in  water  level  =  0.06  ft. 

change  in  tide  level  =  6.68  ft.  x  1.025  ft.  seawater/unit  weight  =  6.8  ft. 

0.06  ft.  =  0.9% 
6.8  ft. 
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9/12/90-10/12/90 

change  in  barometric  pressure  =  0.01  ft 

change  in  water  level  =  0.09  ft. 

change  in  tide  level  =  6.99  ft.  seawater/unit  weight  =  7.16  ft. 

0.01  -  0.09  ft.  =   1.1% 
0.01 -7.16  ft. 

11/12/90 

change  in  barometric  pressure  =  0  ft. 

change  in  water  level  =  0.09  ft 

change  in  tide  level  =  8.58  ft  x  1.025  ft.  seawater/unit  weight  =  8.79 

0.09  ft.  =   1.0% 
8.79  ft. 

The  average  tide  efficiency  for  KB  MP  9  is  1.3% 
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APPENDIX  J 


Calculations  for  Potential  Seawater  Intrusion 
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The  following  equation  is  used  to  determine  the  net  volume  of  water  exiting  the  Pliocene- 
Miocene  aquifer  in  the  channel. 

Y  =  JL  [(Sl£  +  blt)3„  +  (5*£  +  bag  +  (52?  +  b3t)J5 
A     AL        2  2  2 

+  (^  +  b4t)f5  +  &f  +  b4t&  +  <Sj£  +  b5t$]  (J1) 

The  total  volume  of  water  exiting  the  Pliocene-Miocene  aquifer  into  the  channel  is 

calculated  by: 

x = i  [(aaf + blt)30  +  <s*£ + b2t)9  +  (^  +  b3t)J5  +  <a£ + b4t)f5  ] 

A     AL        2  2  2  2  (J2) 

The  total  volume  of  seawater  entering  the  Pliocene-Miocene  aquifer  from  the  channel  is 
calculated  by: 


V_  K   [vm4t2 


Well  site  3 

K  =  2.3  ft/hr 
AL  =  5405.7  ft. 


[<^+h«ia+<s^+i*tgs] 

A     AL       2  2  (J3) 


V  =    2.3  ft/hr  [[l.66  ft/hr  Ghr)2  +  (0ft.)(3hr)  -0]o  + 

A        5405.7  ft.  2 


[-0.8  ft/hr  C9hr)2  +    7.3ft  (9hr)  -    -0.8  ft/hr  (3hr)2  +  7.3ft  (3hr)]l        + 
2  2 


[p.76  ft/hr  (15hr)2  +    -6.7ft  (15hr)  -    0.76  ft/hr  (9hf)2  +  -6.7ft  (9hr)J95         + 

2  2 


[-0.96  ft/hr  (20hr)2  +    19.2ft  (20hr)  -    -0.96  ft/hr  (15hrt2  +  19.2ft  (lShr)]??] 

2  2 


+ 
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V  =  2.3  ft/hr      48.99ft.hr. 
A       5405.7  ft 


V  =  0.02  feet  per  20  hours 
A 


V  =  2.3  ft/hr[[-0.96  ft/hr  (21hr)2  +  19.2ft  (21hr)  -  -0.96  ft/hr  (20hr)2  +  19.2ft  (20hr)]20+ 
A      5405.7  2  2 

[0.33  ft/hr  (24hr)2  +    -8.0ft  (24hr)  -   0.33  ft/hr  (21hr^  +  -8.0ft  (21hr)]2i] 
2  2 

V  =  2i2_ftZhr  2.21  ft.hr. 
A         5405.7  ft 

V  =  0.001  feet  per  4  hours 
A 

The  net  volume  of  flushing  over  a  24  hour  period  is  0.019  ft3/ft2/d 


As  the  nation's  principal  conservation  agency,  the  Department  of 
the  Interior  has  responsibility  for  most  of  our  nationally  owned 
public  lands  and  natural  and  cultural  resources.  This  includes 
fostering  wise  use  of  our  land  and  water  resources,  protecting  our 
fish  and  wildlife,  preserving  the  environmental  and  cultural  values 
of  our  national  parks  and  historical  places,  and  providing  for 
enjoyment  of  life  through  outdoor  recreation.  The  department 
assesses  our  energy  and  mineral  resources  and  works  to  ensure  that 
their  development  is  in  the  best  interests  of  all  our  people.  The 
department  also  promotes  the  goals  of  the  Take  Pride  in  America 
campaign  by  encouraging  stewardship  and  citizen  responsibility  for 
American  Indian  reservation  communities  and  for  the  people  who  1  ive 
in  island  territories  under  U.S.  administration. 
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